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ABSTRACT 



In the Minimal Supersymmetric Standard Model (MSSM), there are numerous sources of 
flavour-violation in addition to the usual Kobayashi-Maskawa mixing matrix of the Stan- 
dard Model. We reexamine the renormalisation group equations (RGEs) with a view to 
investigating flavour effects in a supersymmetric theory with an arbitrary flavour structure 
at some high scale. To incorporate (two-loop sized) threshold effects in the one-loop RGEs, 
we calculate the /3-functions using a sequence of (non-supersymmetric) effective theories 
with heavy particles decoupled at the scale of their mass, keeping track of the fact that 
many couplings (such as gauge and gaugino couplings) which are equal in an exact su- 
persymmetric theory may no longer be equal once the supersymmetry (SUSY) is broken. 
We find that this splitting, which is ignored in the literature, may be larger than two-loop 
terms that are included. In addition, gaugino couplings develop flavour structure that is 
absent without including decoupling effects. A program (to be incorporated into ISAJET) 
has been developed, which includes flavour-violating couplings of superparticles and solves 
the two-loop threshold RGEs subject to specified high scale inputs. The weak scale flavour 
structure derived in this way can be applied to the study of flavour-changing decays of SUSY 
particles. As an illustration, we revisit the branching ratio of the flavour-violating decay 
of the top squark. We find that, in the minimal supergravity (mSUGRA) class of models, 
previous estimates for the width of this decay have been too large by a factor 10 — 25. 
However, this decay rate is sensitive to the flavour structure of the high scale boundary 
conditions. We analyse the consequences of introducing non-universality in the high scale 
soft SUSY-breaking mass matrices and find that under these conditions the partial width 
can be altered by a large amount. 
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Chapter 1 



Introduction 



On the 10th September, 2008, the Conseil Europeen pour la Recherche Nucleaire 
(CERN) began circulating a beam around the 27 km diameter Large Hadron Collider (LHC), 
175 m below the border between Switzerland and France on the outskirts of Geneva [1]. 
Once fully commissioned, the LHC will collide beams of protons travelling at close to the 
speed of light (with 7 ~ 7400) and analyse the products of these collisions. At a total cost of 
around €5 billion the experiments at the LHC will address some of the biggest questions of 
contemporary particle physics. Although a highly successful theory for almost 35 years, the 
Standard Model of particle physics (SM) (see, for example, [2, 3, 4]), which is the current 
state of the art, is generally recognised to be incomplete. 

One potential discovery at the LHC is the Higgs boson, the elementary particle 
which has been postulated to give all other elementary particles their mass. Close study 
of the properties of the Higgs boson will give us some indication of the solutions to the 
inconsistencies of the SM, and it is the only particle in the SM that is yet to be observed 
experimentally. Additionally, the LHC will be looking for hints about the nature of dark 
matter — the unknown material that makes up over 80% of the mass of the universe. 
Contrary to suggestions by the writers of Star Trek, dark matter is not as obvious as a big 
purple cloud and has until now not been seen using standard optical techniques. We know, 
however, from studies of the dynamics of galactic rotation, that typical galaxies contain a 
substantial amount of dark matter out to large distances from the centre [5] . The detectors 
at the LHC may be able to infer the presence of dark matter produced as a result of the 
particle collisions. 

Returning to the problems suffered by the SM, one extremely popular solution 
involves introducing a novel new symmetry, or supersymmetry (SUSY) [6, 7, 8], which links 
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Tabic 1.1: The elementary particles of the Standard Model, including matter fermions and gauge 
bosons. Also indicated is the charge of the matter fermions, Q. 

bosons and fermions, the fundamental groups to which all particles can be assigned. Since 
SUSY is a connection between these two groups, it involves introducing a bosonic partner 
for every SM fermion and vice-versa, thereby approximately doubling the total number 
of elementary particles. In addition, if we insist that SUSY is respected by all particle 
interactions, we find that the additional particles must have precisely the same mass, a 
fact not supported by experiment. Supersymmetry must therefore be broken in a way that 
we can introduce additional mass contributions for the SUSY partners without ruining the 
mechanisms which fix our SM blues (described later). When we break the symmetry in this 
manner, we say we have broken it softly. 

In the SM, the particle content of which is listed in Table 1.1, the interactions 
available to the various matter particles are dictated by choosing certain symmetries that 
these interactions must obey. There are three fundamental symmetries (or gauge groups) 
observed by the SM: colour (or SU(3)), weak isospin (SU(2)) and hypercharge (U(l)) [4]. 
The elementary particles can "feel" a force associated with each symmetry by interacting 
with the "gauge bosons": gluons (which carry colour), charged and neutral weak bosons 
and the photon. Each matter particle is introduced into the theory by specifying how it 
transforms under the various symmetries together with its couplings to the Higgs sector. 
The matter fermions are: leptons, which do not carry colour, and quarks, which interact 
with all three force carriers. 
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The various matter fermions (i.e., quarks and leptons) of the SM are collected 
into three "generations" corresponding to the three rows in the left-hand section of Ta- 
ble 1.1. Within each generation of quarks there are two different "flavours" with differing 
charge. We find that fundamental (or tree-level) gauge interactions between quarks of 
different flavour are extremely restricted in the SM, with so-called "flavour-changing" inter- 
actions only taking place between quarks of different charge. Although this charged current 
flavour-changing interaction mainly takes place within generations, there is also a small 
coupling between quarks of different charge and generation. When we introduce SUSY, we 
find that there is a potential proliferation of flavour-changing interactions, and unless the 
coupling strengths are strongly suppressed such interactions would breach bounds for SM 
flavour- violation obtained by experiment [9]. There must therefore be some mechanism by 
which SUSY contributions to SM flavour-violating processes are suppressed. Our aim is to 
investigate the flavour characteristics of SUSY theories in a general manner. 

Interest in SUSY phenomenology was sparked once it was understood that its 
properties would allow meaningful predictions from theories at scales far in excess of those 
currently available (or likely to ever be available) in terrestrial experiments [10]. SUSY 
can therefore provide us with a framework within which it is possible to obtain the low 
energy consequences of high scale inputs, or ansatze. An early pointer to the possibility for 
a simplified high scale description of the rich structure of weak scale SUSY is the apparent 
unification of the three gauge couplings in the SM. Although they have very different values 
at low scales (~ 100 GeV), they are seen to come together [11] at approximately 10 16 GeV 
if we introduce a SUSY theory ^ at around 10 3 GeV. 

Scale dependence appears in the couplings as a result of our renormalisation pro- 
cedure. The tree-level interaction is augmented by higher-order corrections, coming from 
particles circulating in loops, and the parameters are "renormalised" by defining new pa- 
rameters which are scale dependent. If we try to carry out perturbative calculations at the 
weak scale using parameters renormalised at a hierarchically larger scale, large logarithms 
will render the perturbative expansion invalid. Instead, we sum these large contributions 
using renormalisation group methods [4]. The couplings are said to 'run' with scale and it 
is in this sense that we say the gauge couplings come to a point. 

^ Since the masses of the SUSY particles are essentially unknown, the point at which SUSY influences the 
scale dependence of the couplings is not fixed. However, we have good reason to expect that SUSY particles 
reside in the region of 10 3 GeV, as we shall see in the next chapter. 
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We can derive a system of coupled differential equations called the renormalisation 
group equations (RGEs) which predict how each Lagrangian parameter will depend on the 
scale. RGEs therefore play an important role in the study of the phenomenology of theories 
valid at high energy scales [12] . Although softly broken weak scale SUSY contains a large 
number of parameters in addition to those of the SM, this is a result of our ignorance of 
the mechanism of SUSY breaking. Until this is understood, we resort to models of unified 
theories, where specific physics assumptions lead to models with more manageable numbers 
of independent parameters. The predictions these unified theories make for weak scale 
phenomenology are found by solving the RGEs to obtain the low energy interactions. 

We are therefore in a position to consider various means of introducing flavour- 
violation at some specified high scale, using the RGEs to derive the weak scale parameters 
associated with our favourite theory. Once we have found the weak scale predictions of our 
high scale structure we will be able to investigate the implications for flavour-violation in 
the various decays and interactions of the theory. We must, however, remember that our 
boundary conditions will include not only the values of various unified SUSY couplings, but 
also the size of the weak scale SM interactions which reproduce (to a first approximation) 
the currently observed levels of flavour-violation. 

In many models, the masses of SUSY particles are spread over a wide range. To 
account for the effect that this has on the form of the RGEs, we re-derive the RGEs for 
both the SM [13, 14] and the Minimal Supersymmetric Standard Model (MSSM) [15, 16], 
the SUSY extension which introduces the fewest additional particles and interactions. We 
construct a collection of effective theories with varying particle content depending on the 
scale at which they are valid. The high energy theory contains the same particle content as 
the MSSM, and we remove, or "decouple" , particles from the effective theory at a series of 
thresholds until we obtain the SM as the effective theory valid when all SUSY particles are 
decoupled [17, 18]. 

We find that once we decouple at least one SUSY particle, the effective theory 
'knows' about the SUSY breaking so that couplings that were related in the supersymmetric 
limit are now distinct and have their own RGEs. As a result, the full system of RGEs is 
much larger than usually written in the literature. Specifically, the couplings of the gauginos 
(superpartners to the gauge bosons) to the fermion-sfermion system evolve differently from 
the gauge couplings, while the couplings of the higgsinos (superpartners to the Higgs bosons) 
evolve differently from the usual fermion "Yukawa" coupling. Indeed, gaugino couplings that 
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were flavour-diagonal in the SUSY limit can now develop flavour- violating terms, which must 
be accounted for in the implementation of our decoupling procedure. 

On completion of the program for a complete two-loop calculation of weak scale 
flavour structure, we will apply our results to a sample process — the decay of the lightest 
stop (t\), a superpartner to the heaviest quark. Although the t\ has interactions that would 
permit standard tree-level two-body decays, i.e., t\ — > tZj, where Z; L is a neutralino (one 
of the additional SUSY fermions with zero charge), and t\ — > bW{, where Wi is a chargino 
(charged cousin of the neutralino), the stop may be lighter than the sum of the products, 
in which case such decays would be forbidden. However, if we allow for mixing between 
squarks of different flavours, the t\ may be able to decay via its flavour-changing coupling 
to a charm quark and a neutralino. In this case the decay could compete with, or even 
exceed, the rate of other three-body [19] (or four-body [20]) decays, which are suppressed 
by phase space factors. 

Outline 

The next chapter outlines the various arguments for physics beyond the Standard 
Model. We construct the basic SUSY model, the MSSM, introducing its particle content 
and interactions. After discussing the breaking of both supersymmetry and the electroweak 
symmetry, we briefly mention efforts to detect SUSY particles in colliders and elsewhere. 

In Chapter 3 we write out our notation for the mixing of the quarks and their 
superpartners the squarks. This mixing is responsible for flavour- violating interactions 
and (for SM particles) experimentally well constrained. We also show the interactions 
responsible for squark mixing in the up-type sector. These relations will be necessary for 
our analysis of stop decay. 

Chapter 4 considers the general equations we will use to derive the RGEs of the 
MSSM with full thresholds. In this chapter we discuss our choice of renormalisation scheme, 
our way of writing the various mass eigenstates, and the consequences of allowing for broken 
SUSY in the RGEs. 

The general relations from Chapter 4 are applied to the MSSM in Chapter 5. We 
first write the interaction terms available in softly broken SUSY, and use them to derive 
samples of the full threshold RGEs, contained in full in Appendix A and Appendix B. 

Chapter 6 outlines our method for iteratively solving the RGEs, which was carried 
out by a FORTRAN code, RGEFLAV. Along with a sample input file in Appendix E, we write 
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the general layout of the code, the specifics of the boundary conditions and the iterative 
procedure. Finally, we detail the expected output. 

Once we have described the code used, Chapter 7 presents some sample solutions 
for the various SM and SUSY couplings. We consider examples of all important SUSY 
couplings, and present some surprising consequences of a consistent treatment of flavour 
structure in the RGEs. 

In Chapter 8, the flavour-violating decay alluded to above is considered. We 
compare our result with previous literature and show some scenarios in which there are 
significant changes to the partial widths and branching ratios. This will serve as an example 
of the application of our work. 

Finally, we summarise in Chapter 9. 
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Chapter 2 

Supersymmetry at the TeV Scale 



The Standard Model of particle physics (see, for example, Ref. [2]) is one of the 
most successful physical models to date in that its predictions have withstood experimental 
scrutiny over a wide range of energy scales. Right up to the present day — over 30 years 
since it appeared in its current form — it continues to withstand the assault of sustained 
investigation, being used as the yard-stick for experimental data. 

Having said this, there are a number of cracks in its armour, which perhaps point 
to models with wider scope that include the SM as a low energy approximation. If we are to 
search for a model to replace the SM, we must ensure that the successes of the last 30 years 
are not ignored. A natural path to follow would therefore be to search for a new theory 
which contains the SM in the limit of validity of current experimental results. If we take 
this approach, there are a number of issues with the SM which may point us in the right 
direction: 

1. Perhaps the most obvious limitation is the lack of an inclusion of gravitational ef- 
fects. Since these effects become important when gravitational corrections begin to 
compete with other forces, conventional wisdom states that gravity need only be in- 
cluded at scales in the region of the reduced Planck mass, approximately 10 18 GeV. 
Gravitational effects therefore continue to be of subordinate importance.''' 

2. The Standard Model does not contain a candidate for the source of dark matter - 
the non-luminous matter responsible for a variety of astronomical effects including 
the anomalous rotation curves of most galaxies [5] . 

^Many models have been developed [21] which suggest this scale may be a gross overestimate, and that 
gravitational effects could become important at scales soon to be probed. In this respect, however, we will 
trust conventional wisdom. 
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3. We have no guiding principle to inform our choice of gauge group and particle content, 
particle masses and mixings. All these 'ingredients' of the SM must be chosen on an 
empirical basis. 

4. The theory is not stable when we include scalar fields (which are part of the SM), 
and it is desirable to find some mechanism to prevent scalar masses from suffering 
corrections of the order of the Planck scale. This is the problem of fine-tuning within 
the SM. 

It is due to these issues that we can be certain of the existence of a new theory at energy 
scales we are yet to probe with experiment. This chapter puts forward some arguments for 
weak scale supersymmetry (SUSY) [6, 7, 8] as an interesting candidate for a new theoretical 
framework that contains the SM at low energies. Sec. 2.1 contains a more detailed descrip- 
tion of the fine-tuning problem and considers the scale at which SUSY may exist. This is 
followed by a discussion of why SUSY is considered by some to be a good candidate for the 
new physics. 

2.1 What is the Scale of the New Physics? 

The key question that needs to be addressed when considering the presence of 
physics beyond the SM is at what scale the new physics becomes important. We know that 
no clear signature of new physics has appeared at currently probed scales (up to around a 
few hundred GeV), but the only other clear scale available is the Planck scale. Is there any 
indication of an intermediate scale? To answer this, we must take a closer look at the Higgs 
boson and the fine-tuning problem. 

The Higgs boson, the scalar particle which is introduced to give mass to the other 
fundamental particles, is the final piece in the puzzle for the SM. Current experimental data 
[22] suggests that the SM Higgs boson mass must lie somewhere between 100 — 150 GeV, 
and it is widely believed to be within easy reach of the LHC [1], which is scheduled to 
resume operation at CERN early next year. 

When we consider the radiative corrections to the Higgs boson propagator which 
arise from higher order contributions in perturbation theory, we find a number of quadrati- 
cally divergent pieces coming from both boson and fermion loops. The computation of such 
contributions involves integrating the momentum in the loop up to a cutoff value, A, which 
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can be interpreted as the scale at which the SM is no longer valid and where new physics 
is expected to intervene. The physical Higgs boson mass can therefore be expressed as the 
sum of the mass parameter in the Lagrangian, to# , and the loop contributions, which enter 
with a coefficient, c ~ 1, 

Since the loop contributions are quadratically divergent, the corrections to the squared 
mass of the Higgs boson depend on the square of the cutoff, and therefore, without any new 
physics below the Planck scale, these corrections would be roughly 30 orders of magnitude 
larger than the squared mass itself. 

Many theorists are concerned about the extremely fine cancellation required be- 
tween the separate terms in the sum over contributing processes. In order to allay these 
concerns and make the mass stable, we would like the corrections to the mass of the Higgs 
boson to be of the same order as the mass itself, in other words A ~ TeV. This is an ex- 
tremely exciting suggestion, since we would therefore expect the LHC to not only find the 
Higgs boson, but also begin probing the scale at which we hope to find deviations from the 
predictions of the SM. 

2.2 Why SUSY is a Good Candidate 

We can choose to continue in a number of directions, explaining away the fine- 
tuning problem with various arguments. Some possibilities are: 

1. The Planck scale is not the scale at which gravity becomes important. 

If there exist extra spatially compact dimensions in which gravitational effects can 
propagate [21], we would not observe this propagation at large distances, and gravity 
would appear weaker than it really is. Assuming the effects of the extra dimensions 
start becoming apparent at the TeV scale, this would provide the cut-off we are looking 
for, and we would need to incorporate gravitational effects into our theory. A scary 
prospect. 

2. There is no Higgs particle. The Higgs particle is still to be discovered experi- 
mentally, so perhaps there is another mechanism by which the fundamental particles 
acquire mass. An example of a theory where this occurs is technicolour [23], in which 
the job of the Higgs particle is carried out by composites of a new set of fermions. 
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Unfortunately, technicolour runs into problems since it appears to require more and 
more complicated structures to avoid issues with unreasonably large flavour-changing 
processes. The models become increasingly cumbersome and often conflict with pre- 
cision electroweak data. At any rate, since the scale of the technicolour composites is 
500 — 1000 GeV, these new technicolour particles should be within reach of the LHC. 

3. The world is finely tuned. The cancellations in the contributions to the Higgs 
boson mass are exact to an extremely high precision to all orders in perturbation 
theory. 

If we decide we do not like any of these suggestions, we could try instead to find a way 
to make the quadratic divergences cancel for a symmetry reason, and this is the approach 
which leads to a discussion of supersymmetry. 

When we separate out the divergent one-loop contributions to the Higgs boson 
mass and look at each term individually, we find that the bosonic loop contributions always 
appear with the opposite sign to the fermionic contributions. This is a general feature 
and the reason why SUSY is so interesting. If we arrange for our Lagrangian to obey a 
new symmetry, where the contributions from fermions and bosons are related, we can take 
advantage of the opposite signs to ensure that the quadratic divergences cancel to all orders 
in perturbation theory [10]. This cancellation of divergences is similar to the mechanism 
observed in the massless limit of Quantum Electrodynamics (QED), where large radiative 
corrections to the fermion mass are forbidden by chiral symmetry. 

The essence of supersymmetry, then, is a connection between bosons and fermions. 
We introduce a transformation which connects the two states, and a new operator, Q, which 
generates the transformations, such that 

Q jboson) = | fermion) and Q | fermion) = | boson) . (2-2) 

The operator Q is fermionic and therefore carries half-integer spin. We continue by build- 
ing our theory out of supermultiplets, irreducible representations of the supersymmetry 
algebra. These supermultiplets contain both bosons and fermions, which are referred to as 
superpartners. When constructing a supersymmetric theory with the SM as a low energy 
effective theory, we immediately encounter the problem that there are no observed super- 
partners to ordinary matter. This leads us to roughly double the number of particles so 
that each SM particle gains a superpartner and as we build our theory we must be mindful 
of the fact that we need to arrange for these superpartners to be unobserved to date. 
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In addition to providing an elegant solution to the fine-tuning problem, there are 
a number of other issues which may be eased in a supersymmetric universe. An almost 
accidental consequence of the introduction of SUSY at approximately 1 TeV is to point 
towards some kind of high scale unification of the SM gauge groups. In some SUSY models 
various SUSY parameters are also unified at this scale, which may be an indication of the 
type of physics beyond our simple description of TeV scale SUSY. 

2.3 The Minimal Supersymmetric Model 

Once we decide to pursue supersymmetry as the solution to some of our Standard 
Model problems, we may legitimately ask what a supersymmetrized SM would look like. To 
construct such a model we must promote our fields to superfields (indicated by a caret, as in 
S), promote our particle multiplets to particle supermultiplets, and introduce soft SUSY- 
breaking terms allowed by the symmetry. The resulting theory, known as the Minimal 
Supersymmetric Model (MSSM) [6] , represents the TeV scale limit of our supersymmetric 
universe, and is minimal in the sense that it contains the smallest possible number of 
additional particle states and interactions while keeping the SM as the low energy limit. 

2.3.1 Particle Content 

In much the same way as we choose the gauge groups of the SM, to construct 
a supersymmetric model we need to choose a group for our supermultiplets. So that we 
can step directly from the SM to a supersymmetric equivalent, this is chosen to be the 
same as the SM, namely SU(3)c x SU(2)l x U(1)y- This will fix the gauge supermultiplet 
to contain both the expected gauge bosons from the Standard Model in addition to their 
superpartners, the gauginos. 

We then promote the SM fermion fields to superfields — one superfield for each 
chirality of SM fermion — and use the charge conjugates of the right-handed fermion fields 
so that all superfields are left-chiral. We now have spin 'sfermion' partners for each 
chirality state of each fermion, so we give them a subscript 'L' or 'R' to indicate which 
fermion they are partnered with. It is important to remember that since the sfermions are 
scalar particles they do not have chirality, and these subscripts are just labels to specify the 
fermion with which they are connected. 



11 



The final superfields we need to introduce are for the Higgs sector of the theory. We 
include a doublet Higgs superfield, H u , which contains the usual Higgs doublet and a spin 
1/2 doublet of higgsinos. We assign it weak hypercharge Y = 1 and stipulate that it must 
transform as a 2 under 577(2) l- In the SM, the Higgs doublet would be used to give mass via 
Yukawa interactions to the up-type quarks and its complex conjugate would have Yukawa 
interactions which would give mass to the down-type quarks and the charged leptons. In 
supersymmetry, such Yukawa interactions are specified by the so-called superpotential, /, 
which contains the interactions of various chiral superfields. When we choose a form for 
the superpotential, supersymmetry prevents us from using both a field and its complex 
conjugate. For this reason, the SM Higgs mechanism cannot be precisely copied. Instead, 
we must introduce an additional Higgs doublet, called which transforms as a 2* under 
SU(2)l and is assigned weak hypercharge Y = — 1. We can now see that the labels on each 
Higgs superfield have been chosen to remind us of their role in the superpotential. 

The final choices for the matter and Higgs content of the MSSM are shown in 
Table 2.1. 



2.3.2 Interactions 



Taking the interaction of the particles with gauge bosons to be given by the min- 
imal coupling prescription, our general SUSY Lagrangian is [6]: 
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where is the appropriate gauge covariant derivative, Si and tpi are, respectively, the scalar 
and fermion components of the left-chiral superfield Si. In addition, for each gauge group 
a with corresponding gauge index A, F^ uaA is the gauge field, \ aA are the superpartners 
of the gauge bosons, g a is the coupling strength, t aA are the group generators, and ^ aA are 
constants which are non-zero only for £7(1) gauge groups. The coupling strengths are: g' 
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Superfield 



Content 
Spin 1/2 Spin 



SU(3) C SU(2) L U(1) Y 



Q 



u c 



D c 



L 



E c 



Hi, 



UL 



UR 



en 



d L 



u 



R 



d R 



&L 



'R 



K 
/»?. 



-i 



Table 2.1: The chiral superfields of the MSSM, with particle labels and gauge transformation proper- 
ties. Of course, for each quark and lcpton there arc three generations, only one of which is displayed. 
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for £7(1); g for 51/(2); g s for 577(3). In essence, (2.3) is a master formula, representing the 
complete Lagrangian for renormalisable, supersymmetric gauge theories. 

The first line of (2.3) lists the expected kinetic energies with appropriate covariant 
derivatives. The next line contains the interaction of gauginos with the scalar and fermion 
components of the superfields. It describes how gauginos couple matter fermions to their 
superpartners, and Higgs scalars to higgsinos. The terms on the third line describe scalar 
quartic interactions, of which there are a large number due to the square occurring outside 
the sum on the first term. They will be important for the Higgs potential as well as 
describing the interactions between squarks and sleptons. The first set of terms are called 
D-terms, and the second set, which arise from the superpotential, are called F-terms. The 
final line of (2.3) contains the interactions of a scalar with two fermions. These give rise to 
the Yukawa terms of the SM quarks, and additional higgsino-sfermion-fermion "Yukawa" 
terms. 

The next step is to choose a superpotential for the MSSM, which describes the 
interactions between the chiral superfields and, via the last line of (2.3), will contain the 
Yukawa terms for quarks. We must ensure that the operators we choose result in renor- 
malisable interactions and obey all the gauge symmetries of the SM in addition to super- 
symmetry. Before writing down the most general superpotential possible, we will add an 
additional constraint, not motivated by supersymmetry, which is that all operators must 
conserve baryon number (B) and lepton number (L). Although this is achieved for free in 
the SM, there are valid renormalisable operators which conserve supersymmetry and would 
be included in the most general superpotential possible, but could give rise to, for example, 
a proton decay rate significantly higher than current experimental bounds. 

One method for justifying the absence of these operators from the superpotential 
would be to introduce a new parity that forbids them. This parity can be defined (using s 
for the spin of the field) as 

R = (_l)3(B-L)+2 a ? ( 2 4 ) 

and is commonly referred to as i?-parity [24]. Under this symmetry, all SM particles are 
.R-parity even, and all SUSY partners are -R-parity odd. Note that although requiring im- 
parity conservation means that all allowed renormalisable operators will conserve B and L, 
it is possible to construct .R-parity conserving non-renormalisable operators which do not 
conserve B and L. 
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The introduction of imparity conservation has profound implications for SUSY 
phenomenology. Namely, that SUSY particles must always be produced from SM collisions 
in pairs, that the decay state of a SUSY particle always contains at least one SUSY particle, 
and that the lightest SUSY particle is stable. This final implication, that there is an i?-odd 
particle which is stable with a mass at the weak scale, may turn out particularly useful, as 
we shall see in Sec. 2.4. 

The superpotential for the MSSM with R-parity conservation is: 

/ = niftHaa + (fuhe ab Q?H b u UC + {h)i 3 QtH d aD c 3 + (feh^H^ , (2.5) 

which follows the notation of [6], using the superfields listed in Table 2.1. The couplings, 
(fu,<z,e)«?) are (3 x 3) Yukawa coupling matrices (where i,j label the matter fermion or 
sfermion flavour) and jj, is the Higgs self coupling. The indices a, b are SU(2) doublet labels 
and the matrix e a & is the completely antisymmetric tensor with ei2 = 1. This tensor is 
required in order that we combine the doublets Q and H u , which each transform as a 2 
under SU(2), in an SU(2) invariant manner, and will introduce certain minus signs into 
the resulting Lagrangian. 

Now that we have specified our particle content, supersymmetric interaction La- 
grangian and superpotential, we can turn to the question of why SUSY particles have not 
yet been detected. 

2.3.3 Supersymmetry Breaking 

In a world where supersymmetry is an exact symmetry of nature, the superpoten- 
tial given in (2.5), combined with the general Lagrangian for SUSY, (2.3), would be all we 
need to find the masses and interactions of the various particles in the theory. We would 
be able to calculate all tree-level processes and decays. Unfortunately, the Lagrangian we 
have constructed, being completely supersymmetric, describes SUSY particles which have 
precisely the same mass as their SM counterparts.''' If we lived in this exact SUSY universe, 
we would know about it already from experiment. 

We are thus led to the conclusion that our exact SUSY Lagrangian must be aug- 
mented by terms which are not invariant under the SUSY transformations. In adding these 
terms, we must be careful not to lose the cancellation of divergences which was our initial 

^In this instance the matter fermions and their superpartners would be massless, but it is possible 
to construct extensions to the MSSM with exact SUSY in which they have non-zero masses. For more 
information, see the exercise on p. 160 of Ref. [6]. 
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motivation. Luckily, when the full analysis in carried out, we find that there are certain 
terms which result in mass differences between SM particles and their partners without in- 
troducing new quadratic divergences. These terms only introduce logarithmic divergences 
that are therefore not sufficiently large to resurrect the fine-tuning problem. If we add 
terms to the master formula which break the supersymmetry without destroying the can- 
cellation of quadratic divergences, we are said to have broken SUSY softly. In the MSSM, 
we construct a list of all such soft terms (soft masses as well as additional soft interactions) 
consistent with the symmetries of the underlying theory and add these to the rest of the La- 
grangian [25]. Using, again, the fields listed in Table 2.1, the soft supersymmetry-breaking 
(SSB) terms are: 

£ssb = - {u{ k (m 2 Q ) kl u Ll + d{ k (m 2 Q ) kl d Ll + u Rk (mlr) k iu m + d Rk (m? D ) kl d m 
+ Ak^L)hWu + e[ k (m 2 L ) k ie L i + e Rk {m 2 E ) H e m 
+m 2 Hu h°J h °u + m 2 H h^K + ™ 2 Hd h°Jh d + m 2 Hd K?hj) 

+ {u Rk {aI)kiu L ih° u - u Rk {sZ) k id Ll K + d Rk ( a d)kid L ih° d 

+d Rk ( a d )klULihj + e Rk (a^) M u L ih d + e Rk (a£) kl e u h° d + h.c.} 

- X - [M1A0A0 + M 2 X P X P + M 3 § A g A ] 

- % - [M[ A075A0 + M' 2 \ p ^\ p + M^ Al5 g A ] 

+ {b(h+h^ + h u h° d ) + h.c.}, 

where, in this instance, the roman letters k, I are the flavour indices, and p and A are gauge 
group indices. The Hermitian matrices nig, m^, m^, m| and m 2 E are labelled according 
to the scalar in the interaction and are called soft sfermion mass terms. The soft Higgs 
boson mass terms, and rn 2 ld , and the six gaugino mass parameters M a , M a , are all 

real while the bilinear term, b, is complex. The trilinear terms, a^, aj, aj, are also matrices 
in flavour space. Their elements are frequently written as A. kl f kl which is especially useful 
in the context of models of the SUSY breaking mechanism. We have ignored additional 
trilinears that couple via the conjugate Higgs field and are usually given a coupling c k i 
{e.g., d Rk {c^) k iduh*u)- Although there is no principle forbidding the inclusion of these 
terms, they are negligible in many models. 
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The addition of a collection of SSB terms to the Lagrangian is a useful way of 
parametrizing the means by which the SUSY breaking is communicated to the superpart- 
ners, which is as yet unknown. There is a range of models that attempt to shed light on this 
mechanism [26, 27, 28, 29], that typically involve the transmission of the breaking from a 
so-called "hidden sector" . Such a model leads to a characteristic prediction for the pattern 
of SSB parameters in the Lagrangian renormalised at some high scale, Q high- 
One such model, which we will refer to regularly throughout this dissertation, is 
the minimal supergravity model (mSUGRA) [26] . Within the framework of mSUGRA it 
is possible to obtain a simple set of 'boundary conditions' at a so-called grand unification 
(GUT) scale, Mqut ~ 10 16 GeV, defined to be the scale at which the running gauge 
couplings unify We make the assumption that there is no new physics other than the 
MSSM between the weak scale and Mgtjt, and furthermore that the GUT scale boundary 
conditions consist of a small set of unified parameters. These include: a single value for 
the three gauge couplings, a unified gaugino mass m^, and unified soft SUSY breaking 
scalar masses uiq and trilinear parameters Aq. In addition, we require two parameters from 
the Higgs sector, and we shall see shortly that when we break the electroweak symmetry 
using the MSSM version of the Higgs mechanism we obtain relations between the various 
terms. We will be able to write b in terms of tan (3, the ratio of the two Higgs field vacuum 
expectation values, v u /vd, and will also be able to fix /j? so as to obtain the observed value 
of Mz- Since fixing /j, 2 doesn't fix the overall sign of fi, mSUGRA is therefore completely 
specified by the parameters 

mo, mi/2, A , tan/3, sign(/z) . (2.7) 

Notice that we have said nothing about the flavour structure of many of the cou- 
plings. Since we have a number of scalars with the same gauge quantum numbers, SUSY 
potentially contains many new sources of flavour-violation [10, 30] which can exceed exper- 
imental bounds. All the sfermion mass matrices and a-parameters are matrices in flavour 
space and therefore it is expedient to find a procedure which will enable an estimation of 
their off-diagonal entries given our general high scale inputs from the various models. In 
mSUGRA, all off-diagonal entries in the soft masses are set to zero^ at MquT; but this 
does not imply that they are also zero at the weak scale. The procedure we use to obtain 
phenomenological predictions at the weak scale, is discussed in more detail in Chapter 4. 
^ Explicitly writing the flavour structure of the mSUGRA conditions we therefore have rrig >v D L E — moll 

and SL u ,d,e = Aof u ,d,e- 
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Once we have used this procedure to obtain the weak scale Lagrangian, we will be 
in a position to find the states with definite mass. As in the SM, these physical states are not 
necessarily the same as the current eigenstates. We will describe the fermionic states formed 
by mixing between the charged winos and higgsinos as "charginos", Wi, where i = 1 is the 
lighter and i = 2 the heavier state. Similarly, we will describe the states formed by mixing 
between the bino, the neutral wino and neutral higgsinos as "neutralinos" , Zj, where the 
states are labelled i = 1, ... 4. Squarks and sleptons that share the same quantum numbers 
will also mix, forming two eigenstates of each type and flavour. If we allow for an arbitrary 
flavour structure in all MSSM parameters, these flavour eigenstates will be a mixture of 
both /x, and f R 'chirality' states, and all three current flavour eigenstates, leading to the 
possibility for mixing between a total of six matter sfermions. For example, the lighter 
of the top squark mass eigenstates (labelled with a 1 in analogy with the charginos and 
neutralinos) will be formed from the current eigenstates in the following way: 

ti = a x u L + a 2 c L + a 3 t L + a±u R + a 5 c R + a 6 t R . (2.8) 

The method for constructing this eigenstate from the Lagrangian will be discussed in more 
detail later. 

2.3.4 Electroweak Symmetry Breaking 

In the SM, the gauge group is SU(3) C x SU(2) L x U(1) Y with SU(2) L x U(1) Y 
spontaneously broken to U(l) em . We break the symmetry via a single complex spin zero 
SU(2) doublet, the Higgs field, which acquires a vacuum expectation value (VEV). One 
combination of gauge fields remains massless, and is identified with the photon, and the 
other combinations of gauge fields develop masses through the Higgs mechanism. The 
interaction Lagrangian is then rewritten in terms of the mass eigenstate gauge fields, and 
Yukawa interactions between the Higgs bosons and the SM fermions are used to include 
fermion mass terms. 

In the MSSM, we use the same mechanism to develop gauge boson and SM fermion 
mass terms. Assuming the matter scalars do not develop non-zero VEVs, we minimise the 
scalar potential with respect to the neutral Higgs fields and define = v u and (hP^ = vj, 
where we take v u and v& to be real. The requirement that there is a stable minimum with at 
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least one non- vanishing Higgs field VEV leads us to the following minimisation conditions: 



(m^+m2 d + 2/i 2 )sin2/3 



b = ± '- and (2.9a) 



M " tan 2 /3-l 2 ' iA9bj 



in which tan/3 = v u jv&. Mz is the mass of the neutral combination of gauge fields 
left-over after construction of the photon, and is easily found to be given by M| = 

(g' 2 + g 2 ) K + vj) /2. 

We now have all the information we need (aside from certain subtleties) to con- 
struct the full Lagrangian for the MSSM complete with interactions and masses for all SUSY 
particles. Although this is certainly possible, it would clearly result in a lengthly formula. 
We will return to certain terms in this Lagrangian in later chapters, but before closing this 
chapter we take a brief diversion to consider how SUSY may make its first appearance on 
the experimental stage. 



2.4 Diversion: Experimental Signatures of SUSY 

It would be foolish to close this chapter, having set forth a complex theory full of 
many additional parameters, without at least some comment on the potential for discovery. 
Here we consider two methods of attack. 

2.4.1 Dark Matter Detection 

As previously mentioned, the assumption of i?-parity conservation in the MSSM 
leads us to the conclusion that there is an i?-odd particle which is stable, the lightest 
supersymmetric particle, or LSP. The LSP is either one of the matter sfermions or mixed 
higgsinos/gauginos at a mass smaller than ~ 1 TeV, and its presence is another promising 
aspect to SUSY. SUSY particles produced in the Big Bang, which would eventually decay 
to SM particles and the LSP, may leave us with a relic abundance of LSPs that can account 
for some, if not all, of the dark matter. We are led to ask whether there is any way to narrow 
down the suspects on the list of possible LSPs by making the requirement that they are 
good dark matter candidates, i.e., in agreement with current astrophysical measurements. 

Since all particles produced in the Big Bang were in thermal equilibrium, we can 
calculate their density at the time they ceased to interact with each other. This "decoupling" 
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took place at different times for different particles, and when it happened the total number of 
these particles was fixed. We can reduce their decoupling density in line with the expansion 
of the Universe, providing us with an estimate of their relic density — the level observed 
today. 

If the LSP is a charged particle it will be present in exotic atoms. These atoms 
would look like a very heavy isotope, and should be detected in experiments searching 
for anomalous isotopes. Searches for such isotopes [31] do not find abundances in line 
with the expectations for a TeV scale LSP [32]. We can carry out the same argument for 
a coloured LSP, so that the only viable candidates remaining are the sneutrino and the 
lightest neutralino. The neutralino is currently favoured since sneutrinos are bounded by 
two complementary experiments. Doublet sneutrinos heavier than 25 GeV would already 
have been observed in direct DM detection experiments, and a sneutrino lighter than this 
would have been observed in detailed observations of the properties of the Z boson at LEP. 

There is still much excitement about the possibilities for a neutralino dark matter 
solution. In one way, we can think of the evidence for dark matter as being a strong 
signal pointing towards a SUSY model with a neutralino LSP, and it is fortunate that the 
neutralino is indeed the LSP in a wide range of models. The WMAP experiment has made 
a very good estimate of the dark matter relic density which strongly constrains the MSSM 
parameter space. Perhaps this will be the area of experimental enquiry that will be the first 
to provide compelling evidence for physics beyond the SM. 

2.4.2 Colliders 

Colliders have so far been unable to offer evidence for the presence of TeV scale 
SUSY. They continue instead to place bounds on sparticle masses and reduce the available 
parameter space for the various models of SUSY breaking. Of course, the extended reach 
of the LHC and any possible linear collider will either push these bounds further or offer 
our first glimpse of supersymmetric particles. 

The cosmological arguments from the previous section, in which we used the as- 
sumption of .R-parity conservation in the MSSM, strongly point toward the presence of a 
neutral LSP. In addition, imparity conservation also implies that sparticles can only be 
produced in pairs, and only then if there is sufficient energy in the collision. The current 
bounds on SUSY particle masses would suggest, therefore, that the energies available at 
current or future colliders will only be sufficient to produce two sparticles per event. These 



20 



two particles will undergo cascade processes and eventually decay into two LSPs, which, 
only interacting with normal matter through the exchange of other heavy (SUSY) particles, 
will escape detection — the LSP will behave like a heavy neutrino. 

Since (even in hadron colliders) we know that the total transverse energy for the 
initial state particles is zero, a clear indication of a SUSY event would be an imbalance in 
the transverse energy of the final state. Although a missing transverse energy signal can 
also be due to SM processes such as neutrino production, in this case the amount of missing 
energy is usually significantly smaller. Detection of missing transverse energy events is 
considered an important ingredient in SUSY searches at colliders. 

It should be noted that if the LHC is unsuccessful in finding SUSY, some of the 
elegant arguments in favour of supersymmetry at the weak scale will no longer apply. This 
would not rule out SUSY as a candidate for new physics at some high scale, but its rich 
phenomenology would be out of reach for the foreseeable future. For an overview of the 
current state of collider searches and the discovery potential at future colliders, the reader 
is encouraged to see Chapter 15 of Ref. [6]. 

If the LHC does indeed produce SUSY particles, the next step would be to begin 
to measure the particle spectrum and coupling strengths, and eventually try to uncover the 
origin of the SSB parameters. This would be a complex, but rewarding challenge. 
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Chapter 3 

Flavour-violating Interactions 



If the soft SUSY-breaking matrices of the MSSM (written in the basis in which 
the corresponding SM fermions have definite mass) contain large off-diagonal entries, they 
lead to flavour-changing effects that are not present in the SM. Up to now we have avoided 
any mention of the structure of these terms, but their inclusion in the MSSM is potentially 
a major problem for SUSY model builders. An attractive feature of mSUGRA is that the 
GUT scale inputs in (2.7) contain no flavour off-diagonal couplings, effectively reducing any 
additional flavour-violating effects at the weak scale to acceptable bounds. We shall see that 
we can adapt the mSUGRA framework to carefully control the sources of flavour- violation 
in the weak scale effective theory. 

This chapter briefly reviews flavour-changing effects in the quark sector before 
describing in more detail the situation in the MSSM and construction of the squark mass 
eigenstates, which dictate the size of flavour-violation in the squark sector. 

3.1 Quark Flavour Mixing 

It is well known that in the Standard Model, the current eigenstate quarks (the 
states that take part in weak interactions), with fixed quantum numbers, are not simulta- 
neously mass eigenstates. The K + , for example, is a combination of a u quark and an s 
anti-quark, and is observed to decay 63% of the time to p + v^. Since the interaction La- 
grangian of the SM only contains transitions between quarks of differing flavour from the 
same generation (i.e., u d, c <-> s and t b), we know from the observation of both 
ir + (ud) and K + (us) decays that the quark mass eigenstates must be a superposition of the 
current eigenstates. 
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When we write down our interaction Lagrangian, there is no principle to prevent 
the Yukawa matrices, which couple left- and right-handed quarks with Higgs bosons, from 
being off-diagonal in the current basis. Therefore, our Lagrangian may contain transitions 
between the various current eigenstate quarks of the same type, which is exactly what we 
need to reproduce the above behaviour. In order to find the mass eigenstates we write down 
all the mass terms in our theory and diagonalise them using a bi-unitary transformation 
which describes the connection between the current basis and mass basis as follows. 

3.1.1 Quark Yukawa Matrix Diagonalisation 

From the master equation, (2.3), and the superpotential, (2.5), the mass term for 
the up quarks is 

C 3 -u 3 {t u )lh° u P LUl , (3.1) 

where the Higgs field acquires a VEV, v u . The expression is altered slightly in the SM 
since the field is not the Higgs field of the Standard Model, and the subsequent relation 
between the MSSM Yukawa couplings (f u ,d,e) an d SM Yukawa couplings (X u ,d,e) wm be 
discussed in more detail in Ch. 4. 

Eq. (3.1) is written in the current (or interaction) basis, where the Yukawa matrix 
may have large off-diagonal entries. We now diagonalise the Yukawa coupling by writing 

ff as = Vl(u)f u V* R (u) , (3.2) 

where the rotation matrices Vl,r(u) are unitary, and requiring that f„ iag has vanishing 
off-diagonal entries. Denoting the mass eigenstates with a superscript M, it is clear from 
(3.1) that the relation between current and mass basis quarks is 

U L i = [V 'L(u)]ijULj 5 U Ri = [VR(u)]ijU Rj . (3.3) 

We can write similar rotations for the down- type quarks using rotation matrices V L,R,(d) 
to obtain an equivalent relation to (3.2) for f^. 

All we have done is to change the basis. To see whether the matrices V l,r{u,<I) 
are physically meaningful, we must consider each term in the Lagrangian. For the SM, 
the Higgs interactions are diagonalised simultaneously with the quark masses. In addition, 
the neutral current interaction, which is by definition flavour diagonal in the current basis, 
remains flavour diagonal in the mass basis due to the unitarity of the rotation matrices. It 
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is only the charged current interaction which causes us trouble. We write 

C 3 --^urfW+SjiPLdi + h.c. (3.4) 
v2 

which only allows transitions in the current basis between flavours of the same generation. 
We transform to the mass basis using (3.3), and the equivalent transformation for the 
down-type quarks, to obtain 

C 3 -^f^W+lV L (u)}] k [VUd)UPLdf + h.c. (3.5) 

Since the right-handed quarks do not take part in the charged current interaction, the 
right-handed rotations do not enter. In addition, only one combination of the left-handed 
matrices, namely V| y (u)Vi(c/), is physical. This combination is usually referred to as the 
Kobayashi-Maskawa (KM) matrix [33], i.e., 

K = v[(u)V L (d). (3.6) 

The SM Lagrangian in the mass basis only depends upon the V-matrices via the 
KM matrix. We can therefore see that all the flavour structure in the quark sector of the 
Standard Model is dictated by this combination alone. Since it is formed of two unitary 
matrices, it is itself unitary and would therefore be completely described by three angles 
and three phases. Using our freedom to make field rotations, we can reduce this to three 
angles and a single phase, and the KM matrix is usually parametrised as [34] 

/ C12C13 S12C13 s 13 e~ l5ls \ 



K = 



-S12C23 - Cl2S23Sl3e J<513 Ci 2 C 2 3 ~ S 12 S 2 3Sl3e lSl3 S 2 3C 13 



(3.7) 



V S12S23 ~ Ci2C23Sl3e i<513 -Ci 2 S 2 3 ~ «12C23Sl3e l<513 C23C13 / 

For the numerical results that we present, we take s\ 2 = 0.2243, S13 = 0.0037 and S23 = 
0.0413, and S 13 = 60°. 

We note for later reference that since flavour-violation in the SM enters only via 
the KM matrix, physical observables can only depend on Vl(u) and Vx,(d) through K as 
specified by (3.6) and must be independent of Vr(u) and V#(d). 

3.1.2 Choice of Basis 

Although it is tempting to choose to work from now on in a basis in which all the 
Yukawa matrices are diagonal^, this approach is not a practical option. Since the quark 



'We shall see that if the Yukawa coupling matrices are taken to be diagonal at some specific scale, for 
example m t , they will not in general be diagonal at other scales. 
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mass basis is defined with different transformations for the up-type and down-type quarks 
this choice would not preserve the SU(2)l symmetry. Instead, we choose a current basis in 
which only one of the Yukawa matrices, f u or i d , is diagonal. 

In order that the underlying supersymmetry is preserved, we must rotate the 
squark multiplets by the same transformations as the quarks, although this does not mean 
that the squarks will also be in their mass eigenstate basis. As a consequence, the SSB 
masses and a-parameters transform from the general current basis to the basis in which one 
of the Yukawa matrices are diagonal as 



where we use a (q) to denote the rotation of the SU (2) l doublet squarks, and set q = u(d) 
in the case that the up (down) Yukawa couplings are diagonal. 

With a view to our study of the flavour- violating decay of the t\ in Ch. 8, we choose 
here a "standard" current basis in which it is the up-type Yukawa coupling matrix that is 
diagonal. In so doing, we have fixed ~Vl,r(u) = Vfl(d) = 11, and Vi(d) =Kas required by 
(3.6). The relation between the various parameters in the "standard" current basis and the 
basis in which the down-type Yukawa matrices are diagonal is relatively straightforward to 
derive. We drop the superscript 'M', instead using (u) and (d) to signify whether f u or f d 
is diagonal, and find that 



Remember that the ~Vl,r(u, d) diagonalise the quark Yukawa matrices and are fixed by the 
relation between the current basis and the quark mass basis. In this regard, the final relation 
above follows from the fact that the right-handed squarks are singlets under SU(2). 



(a Uid ) T = V R (u,d)( a ™ d ) v{(u,d) 
m^=V L (g)K) M V^) 
m u,D = Vr(u, d) (ml !D ) M \ j R (u, d) , 



(3-8) 
(3-9) 
(3.10) 



(f u , d y (d) = (f u , d y (u)k 

(au,df(d) = (a M ) T (n)K 
m 2 Q (d) = Kt m ^(n)K 

m uM d ) = ™-u,d( u ) ■ 



(3-11) 



(3.12) 
(3.13) 
(3.14) 
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3.2 Squark Flavour Mixing 



In the MSSM the squarks have couplings to the Higgs sector equal to those of the 
quarks. They therefore likewise have flavour off-diagonal Yukawa terms in the Lagrangian. 
In addition, since we know SUSY is broken, the squarks also have flavour-violating couplings 
coming from SSB squark mass terms and a-parameters. If there is no new source of flavour- 
violation in the MSSM, the KM matrix still governs the size of the mixing, and the theory 
is still independent of the specific choices for and Vr. However, the squarks will not 
necessarily be simultaneously diagonalisable with the quarks due to the SSB terms in the 
Lagrangian. 

Large flavour dependence in the sfermion mass matrices could lead to large contri- 
butions to flavour-changing neutral current interactions, for example b — > S7 and fj, — > ej. 
The most stringent constraints on the flavour structure of supersymmetry come from studies 
of K° — K° and D° — D° [35] mixing. If we naively estimate the off-diagonal entries of the 
SSB squark mass matrices to be of the same order as the diagonal entries, we would achieve 
SUSY contributions to the mixing that violate current bounds from experiment by many 
orders of magnitude unless sparticles are very heavy — excluding the model as a solution 
to the hierarchy problem. This is known as the SUSY flavour problem. 

In addition to contributing to SM processes, the off-diagonal elements of the SSB 
masses and a-parameters will result in squark eigenstates which can couple to all three 
quark flavours. In particular, the squarks may experience flavour-violating decays including 
our sample decay in which the squark eigenstate that is predominantly stop can decay to a 
charm quark (or an up quark) and a neutralino. 

In order to calculate the size of flavour- violating decays of squarks, we must work 
with mass eigenstates of both quarks and squarks. We construct a matrix of all mass terms 
for squarks in the Lagrangian in the basis where the up-type quarks are diagonal, and then 
diagonalize this matrix to find the squark eigenstates. Note that in the next section, we drop 
the basis labels on the couplings and fields on the understanding that we write everything 
in the "standard" current basis. This choice allows us to easily identify the coupling of the 
c quark mass eigenstate to the up- type squarks. 
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3.2.1 The Squark Mass Matrix 

The current basis Lagrangian for the MSSM contains several mass terms for the 
up-type squarks that can be collected into a (6 x 6) matrix, {■M-'i) ^ ia ^j b y such that 

£3-(4^L)(Ml) M0 - 6) ^2) ' (3 ' 15) 

where the indices i,j label left-handed squarks and a, b label right-handed squarks. These 
indices can be expanded out to show the current basis squarks to be uu = (ul, cl, £l) and 
URa = (uR,CR,tR). Considering the individual contributions to this matrix, we can write 
{^■Vj {ia){jb) m terms °f x 3) sub-matrices 

(*a),- w -»=C {MIl) ? iMla) A, (3.i 6 ) 

and find that 

(Ml L ) tJ =(m£).. + vl{rX) l3 + (g - f ) (vl - vj) 6 l0 (3.17a) 
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(M RR ) ab = (™u)ab + <{*uK) ah V d) S «* ( 3 - 17b ) 

{•M 2 LR ) ib = -v u (& u )* ib + v d (fi*f u )* ib . (3.17c) 

The Yukawa coupling terms will be (almost) diagonal in the quark mass basis'^ 
and therefore not contribute significantly to the flavour- violating terms. If uiq, mfj and 
a M are all diagonal in the same basis as the Yukawa couplings there is no mixing between 
generations, but we still have mixing between sfermions with the same charge in the same 
generation. In this scenario, we can write the mixed sfermions in terms of the unmixed 
sfermions via a (2 x 2) mixing matrix 




cos9f sin9f 
— sinOf cos9f 




(3.18) 



which depends on a single rotation angle, Of. 



^We will discuss the extent to which the Yukawa matrices deviate from diagonal in later chapters. 
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As mentioned in Ch. 2, there is a further simplification in many models where the 
terms A u ^ t are related to the matrix a u by the following relation: 

/ KU o o \ 

a u = A c f c • (3.19) 
V A t ft J 
Using this substitution, the mixing angle for up-type squarks is given by 



t&nOf = 



m 2 j +m 2 f + M| cos 2/3 (§ - | sin 2 <%) 



to- 

71 



rrif (—Af + [i cot /?) 
and the masses of the eigenstates fa and fa are 

™l = ^ f™l + m 7 ) + T M i cos 2/3 + to 2 , 
/i,2 2 \ jl jrJ 4 J 



=F 



sin 2 9 W 



(3.20) 



1/2 (3.21) 



+ mj(-Af + ficotPY 



2 V "7l ' A V, 1 3 

It is straightforward to derive similar equations for the masses and mixings of the down-type 
squarks and sleptons. 



3.2.2 Flavour-violating Decays of Up-type Squarks 

It should be clear that in general, if hIq, and &u are not all diagonal in the 
same basis as the Yukawa coupling matrices, all three flavours will mix with each other. In 
this case we label the mass eigenstates with the name of the current basis squark to which 
they are most closely aligned, and with a number indicating their relative mass. For the 
up-squarks, t\ is defined as the lighter of the two squarks with the most amount of stop, c\ 
is the lighter of the two squarks with the most amount of scharm and so on. 

To find the mass eigenstates, we diagonalise the (6 x 6) matrix given in (3.16) by 
rotating the squarks by a unitary matrix, U. This matrix can be split into left and right 
sections 

U R 

where both IAl and 14 r are non-unitary (3 x 6) matrices. Denoting the squark mass 
eigenstates with the superscript 'M' and the label a, we can write the generalisation of 
(3.18) as 

/ n, r ■ \ I 7 J r . \ 

(3.23) 



U = 



(3.22) 
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so that u„ = tf Lai ULi + U' Raa ujia- When we transform the Lagrangian into the squark 
mass basis, we will use the relation: 

(fiti.fili.)=fi? t ("i«'"*«) • (3 ' 24) 

and define the squark mass basis to be in the order (ti,t2,ci,C2,ui,U2) so that the t\ is 
the state uf 1 (which is not necessarily the lightest of the up-type squark mass eigenstates). 
Eq. (3.23) is essentially the inversion of (2.8), which symbolically lists the components of 
one of the mass eigenstates, the t±. We can confirm that if the up-type squarks are not 
simultaneously diagonalisable with the up-type quarks, i.e., the squark mass matrix has 
non-zero off-diagonal entries in the "standard" current basis, the t\ will contain a certain 
amount of c (in addition to t and u) and therefore be able to couple to, for example, a c 
quark, corresponding to non-zero values of either one or both of 02 and 05 in (2.8). 

In the mSUGRA framework, the SSB mass matrices are proportional to the unit 
matrix at some high energy scale (Mqut) and the a-parameters are proportional to the 
corresponding Yukawa coupling matrices, such that these Yukawa couplings are the only 
source of flavour-violation in the theory. However, the form of the SSB parameters will be 
altered by radiative corrections, as embodied by their renormalisation group evolution, and 
the resulting weak scale terms will not generally be simultaneously diagonalisable with the 
corresponding Yukawa matrices. 

In order to estimate the effect of radiative corrections on the squark mass terms, 
and calculate the squark mass eigenstates, we must consider the dependence on scale of the 
various couplings in our Lagrangian. Our discussion therefore turns to the renormalisation 
group equations, which will tell us how the relevant weak scale couplings are related to our 
GUT scale inputs. 
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Chapter 4 

Renormalisation Group Equations 



With the MSSM, we have put in place a model with mass scales that are stable 
under radiative corrections [10], which, however, contains a large number of additional 
parameters. We can hope that this proliferation of parameters simply reflects our ignorance 
of the mechanism of SUSY breaking, and that these arise from a simpler underlying theory 
valid at some high scale. In the absence of knowledge of such a theory, our approach, 
therefore, is to introduce more manageable frameworks, valid at the high scale, which 
predict the MSSM masses and couplings. 

It is tempting to imagine that we could use the GUT scale parameters of a model 
such as mSUGRA to directly calculate the interactions and decays that take place at the 
weak scale. However, we are unable to carry out such a calculation due to the dependence 
on scale of the fields, couplings and masses of the theory. This scale dependence means 
that any perturbative calculation in terms of some coupling, g, renormalised at the high 
scale, would contain large logarithms, In (Mqut/Mz), which would invalidate perturbation 
theory through the appearance of powers of |- In (Mqut /Mz) ~ 1 for a fixed value of g 2 . 
We must therefore obtain the weak scale parameters before computing the various transition 
amplitudes, since in this manner we will sum all the large logs that would have invalidated 
the perturbative expansion. 

Scale dependence in the parameters arises as a result of our "renormalisation" 
procedure. Although the so-called 'bare' parameters are perturbatively divergent quantities, 
we shift the divergence into counterterms. Once we have reparametrised the theory, the 
dependence of the renormalised quantities on the scale, Q, is described by the Callan- 
Symanzik (3- function, defined as 
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and we calculate the /3-functions for the various parameters by considering all diagrams that 
contribute to the renormalisation. We then say that the coupling constants and masses in 
the theory 'run' with scale according to the renormalisation group equations, which provides 
us with a convenient method to obtain the weak scale parameters from our high scale inputs. 

Given the RGEs for all dimensionless (gauge and Yukawa) and dimensionful (soft 
mass and trilinear) parameters we will be able to study the flavour-violation in models 
with arbitrary values of high scale SSB parameters. The RGEs for the SM [13, 14] and 
the MSSM [15, 16] are already known at the two-loop level. If we insist on achieving two- 
loop accuracy, however, we must also be prepared to model the transition from a fully 
supersymmetric theory, with RGEs as in the MSSM, to a SM theory with SM RGEs. The 
effect of this transition in the one-loop RGEs can be as large as two-loop effects and perhaps 
even larger if there is a sizable mass difference between the various SUSY particles. The 
task of estimating the effect of including this transition has been attempted in existing 
literature [17, 18]. We detail our extension to this work by introducing the procedure for 
deriving the RGEs in a general field theory and then moving on to consider the problems 
encountered when applying these RGEs to the MSSM-SM transition. 

4.1 Gauge coupling running and unification 

The derivation of the RGEs for the gauge couplings in a general theory is laid out 
in many places [4, 13]. For a theory with a number of scalars and fermions in different 
representations, the one-loop general result is [36] 



W &li-]oop = V 



Y c ^~t £ s(Rf)-\ £ S(R S ) 

fermions scalars 



(4.2) 



where C{G) is the quadratic Casimir for the adjoint representation of the associated Lie 
algebra, and S{Rf) and S(Rs), respectively, are the Dynkin indices for the representations 
Rf, Rs under which the fermions and (complex) scalars transform. For the Lie algebra of 
SU(N), C{G) = N, while S(R) = 1/2 for the fundamental iV-dimensional representation, 
and S(R) = N for the adjoint representation. For the U(l)y gauge coupling g', C{G) = 
while S(R) = (Y/2) 2 . 

Of course, we obtain different RGEs for the gauge couplings with SM field content 
compared to MSSM field content and at some scale our theory must change from pure SM 
running to pure MSSM running. With i running from 1, ... 3 to indicate the different gauge 



32 



60 



50 



40 



30 



52 20 



10 - 




10000 



_L 



le+08 le+12 
Scale, Q (GeV) 



le+16 



Figure 4.1: The running of the inverse gauge coupling strength, a~ , for all three gauge groups, %. 
The solid lines contain a switch to MSSM running at 1 TeV, the dashed lines show SM running over 
all scales. Two lines are plotted for aj 1 to indicate the limits of experimental uncertainty. 
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(4.3) 



(4.4) 



(4.5) 



In (4.3) we have used a scaled U(l) gauge coupling^, g\, which is related to g' by g\ = 
s/5/3g'. 

In Fig. 4.1 we give an example of the one-loop running of the gauge couplings in 
two simple scenarios. Beginning with the experimental values of the gauge couplings at 
Mz, we evolve their values up to the GUT scale using the RGEs. We plot the running 
of a^ 1 where a» = gf/Air, and have drawn two lines for g% to indicate the relatively large 
experimental uncertainty in this value. 



^This change is motivated by grand unification considerations which fix the normalization of the (7(1)} 
generator. 
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The dashed lines show the evolution in the SM all the way to the high scale. In 
contrast, the solid lines use a transition at an approximate SUSY scale, Msusy = 1 TeV 
where SM running gives way to the MSSM. Remarkably, the couplings appear to unify (or at 
least come extremely close to unification) at a scale somewhere in the range Q = 10 16 GeV 
[37]. This outcome is clearly dependent on where we choose to implement the change to 
SUSY running, but the general result strongly suggests some sort of unification of the 
couplings and hence new physics at this high scale. 

Note that even with the uncertainty in a s (Mz) taken into account, there is no hope 
of unification within the SM. The appearance of (near-)unification only when we introduce 
MSSM effects at Msusy is regarded by some as indirect evidence of SUSY at the weak 
scale. 

We mention briefly here that the one-loop RGEs (in the full MSSM) for the gaugino 
couplings, Mi, M2 and M3, are proportional to the /^-functions for their respective gauge 
couplings. Since the constant of proportionality is the same in each case, this has the 
consequence that, in models where both the gauge couplings and the gaugino masses unify 
at the GUT scale, we obtain the RGE invariant relation 

Mi ~ M 2 ~ M 3 ' [ ' ' 

which is valid at one-loop order. In Ch. 7 we shall consider deviations from this relation in 
the two-loop MSSM with thresholds. 



4.2 RGEs for a General Field Theory 

Our starting point for including threshold effects in the RGEs for the MSSM 
will be the general results in the seminal papers by Machacek and Vaughn [13] for the 
dimensionless couplings, and the paper by Luo, Wang, and Xiao [38] for the dimensionful 
couplings. These RGEs are written in a general gauge field theory, which is not necessarily 
supersymmetric, and uses two-component notation for fermions and real scalar fields. Since 
most phenomenologists are more familiar with four-component spinor fields, we convert 
their work to apply to four-component fermions (Majorana or Dirac) and complex scalars. 

We denote our four-component spinors by ^ Di for Dirac fields and \& Mi for Ma- 
jorana fields. The indices, k, ... are fermion field indices which carry information about 
field type (quark, lepton, gaugino, higgsino) and also carry information about flavour and 
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colour. Similarly, complex scalar fields are denoted by $ a where a, b, ... are scalar indices 
running over squark, slepton, and Higgs fields. In the context of the MSSM with con- 
served .R-parity we will write our general Lagrangian density without including baryon or 
lepton number violating Yukawa terms, or terms which only appear in the MSSM after the 
spontaneous breakdown of electroweak symmetry. With this in mind, the four-component 
Lagrangian takes the form 



£ (4) = ^rfD^j + (D^ a y (D»* a ) - -.F^aF^ 



1 



(m x ) jk ^Mj^Mk + i (m' x ) k ^Mjlb^Mk 



+ 



Bab^a^b + h.C. 



™abK®b 



+ (Va) jk ^DjPL^Mk^a + (Wa) jk *MjPL*Dk*l 

+ \ ( X a) jk ^MjPL^Mk^a + \ (X-l) jk ^MjP^Mk^a + ^-C 



+ 



2! 



& a 1I abc <S> b <S> c + h.c. 



1 1 
2! 2! 



Aabcd$t®t®c$d 



3! 



h.c. 



(4.7) 



In (4.7), there are a number of matrices with a variety of symmetry properties under in- 
terchange of indices. The matrices U*, U 2 ,, V a and W a have no particular symmetry or 
Hermiticity properties, but X* and X 2 are symmetric under j <-> k because of the symmetry 
properties of the Majorana spinor bilinears [6] . Finally among the collection of dimension- 
less parameters, the A a h c d and A.' abcd are symmetric under a <-> b and/or c <-> d for A and 
under interchanges of b, c, d for A'. 

The dimensionful parameters in (4.7) also respect certain symmetries. The scalar 
mass squared matrix, m 2 , is Hermitian, while the scalar bilinear matrix, K, is symmetric. 
The CP-conserving and CP-violating fermion bilinears, mx and respectively, are both 
real and symmetric, and finally the trilinear coupling H a fe c is symmetric under b <-> a We 
continue by writing the four-component general RGEs for the dimensionless parameters 
from (4.7). 
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4.2.1 Dimensionless Parameters 



After converting the general RGEs in [13] to four-component notation, the (3- 
functions for the Yukawa coupling matrices in (4.7) are found to be [39], 



( 47r ) 2 ^Lioo P =2 



(usujt + ugu? + v 6 v!)ui 



+uifuj t u 6 1 + uj t ug + wjw { 



+ 2 



uju^ug + ugu^tuj + v b x 2 Jw b 

+ VlTr I (Uj+Ui + UfU 2 ) + I (xjtxi + X^X 

+ UfTr { (u?Ui + U?Uj) + \ (x*X\ + X?Xj) } 
- 3 5 2 [U^(F) + C* (F)Vl] 

+ \ v a (ytvl + x 2t w 6 ) + 2ul wt w b 

+ U 6 1 Tr{wtw 6 + vtV (l } ; 



1 



w ii-ioop =2 k u w + U » U * T + v * v s J u ° 



T2 TT 2t 



+U 

+ 2 



2 (u b lf uj 



+ l#U? + w. 



+ ^Tr{(u^ + uM) + l(v^ 

+ U 6 1 Tr{(u^ + uM)+^ 
- 3 5 2 [U^C^F) + Cf(F)V 2 a ] 

+ ^ (u^wj + V 6 X^) W a + 2v 6 vtu; 

+ U 2 Tr{w|;W a + Vtv fe } . 



a ^bj 



x,'x^ + x 



(4.8) 



(4.9) 
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Since the Yukawa matrices represent the coupling of a scalar to two Dirac fermions, we will 
use these equations to derive the familiar f u ^ e Yukawa coupling matrices. In addition 

1 



(4tt) /3xi I 



+ 2 
+ X 6 



w 6 w|; + vjv; + xjxj* + x 2 x 2t J xi 
+xj (w 6 * + vtv 5 + xj'xj + x 2t x 2 ) 
w 6 u 2t v fe + v^u 2 *w^ + x£x 2 tx 2 + X 2 x 2 t X i 

^Tr { (Ujtui + U 2 tU 2 ) + I (Xjtxi + X 2 t X 2 ) ^ 



X 2t Xi + XifT X, 



XgTr{(u^ + uM) + ^ 
3 5 2 [X^(F) + Cf(F) X i] 

2 [( W b U 6 2t + X 6 V 1) V * + V « (U^*W^ + V b *X*) 
W 6 Wt X i + XJW:W[] + X^Tr {Wt W 5 + VjV a } 



+ 2 
+ 2 



(4vr) 2 /3x2 L loop =\ [(W 6 WJ + V£Vj + XjxJ* + X 2 X f) X 2 



-X 2 (w 6 *w^ + vjv 6 + X^Xl + X 2t X 2 



+ 2 



w b u^t v fe + vlul*wl + x^tx 2 + x 2 x^ 

+ X 2 Tr { (U 2 tu 2 + UM) + \ (X^X 2 + XK) } 
+ X^Tr { (l^U 2 + U^U 2 ) + \ (x^X 2 + X^X 2 ) } 
-3 5 2 [X 2 C^(F) + Cf(F)X 2 ] 



+ 



1 



2 
+ 2 



(4.10) 



(4.11) 



(vlul* + x 2 w b *) + w a (u^Vj, + w+x 2 

V^V:X 2 + X 2 VtV b ] +X 2 Tr{wtw a + Vtv b } , 

which are the couplings of scalars to two Majorana fermions and will therefore describe 
the RGEs of Higgs-higgsino-gaugino interactions. Note that the property X„ 2 = (Xa' 2 ) T , 
exhibited by the Lagrangian (4.7), is observed by /3 x i,2. Also, in (4.8)-(4.11) we have 
separated out terms which vanish in the context of the imparity conserving MSSM and 
written them on the final two lines of each equation. These terms vanish because W a and 
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V a are zero for a values for which Uj' 2 and xi' 2 are non-zero. Finally, 
(4vr) 2 Pv^-^v =\ [( U W + U 2 U 2t + V b Vj) V a 

+v a (w* b w{ + vlv b + x?xl + X 2t X 2 )" 

+ 2 [uj WtX 2 + U 2 Wtx^ + V 6 W:W^ 

+ V 6 Tr { Wtw 6 + V+V a } - 3<? 2 [V a C#(F) + Cf (F)V a ] 

+ ^ [(u 2 Wj + VftXjt) (Xi + X 2 ) + (Ui + U 2 ) (u^Vft + W+X 2 ) 
+ 2 [u 2 (U^t + U 2 t) V b + V 5 (X*t + X 2 t) Xj 
+ V 6 Tr{ 



itf fui + u 2 )+ uit + u 2 .t U 2 



+ : 



xjt (xi+x 2 ) + (xit + x 2 t)x 



(4.12) 



(4^) 2 ^Ji-ioop =^ [(W 6 WJ + +XiX 6 lt + X 2 X 2t ) W a 

+w a (uj+uj + u 2t u 2 + wlw b ) 



+ 2 



v[v;w 6 + xM^ + xMuJ 



+ W 6 Tr { WjW a + Vt V b ] - 3g 2 [W a C%(F) + Cf (F)Wj 

+ ^ [(w b u 2t + xjvj) (ui + u 2 ) + (xi + x 2 ) (vju^ + X 6 % 6 ) 

+ 2 [w 6 (itf + U 2 t) + Xj (x% + X 2 t) W 6 " 
+ W 6 Tr { [U 2t (Ui + U 2 ) + (l# + U 2 t) U 6 1_ 

* x 2 t(xi + x 2 ) + (xit + x 2 t)x( 



(4.13) 

are the RGEs for squark-quark-gaugino and squark-quark-higgsino interactions. In this 
instance, the last three lines of (4.12) and (4.13) contain terms which will be zero when 

12 12 

applied to the -R-parity conserving MSSM because each of U ' and X a ' vanish for the 
values of a for which V a and W a do not. 

Eqs. (4.8)-(4.13) contain the quadratic Casimirs C^iF) and Cf(F). Along with 
C<jf(F), to be used in the next section, they are given by C2(F) = t A t A , with t A the group 
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generators for the reducible representation that includes all fermion fields. The superscripts 
L, R, M indicate the contributions from left-handed, right-handed and Majorana fermions 
respectively. Flavour is carried not only by some of the indices of the various matrices, but 
can also be carried by a, which may also carry a colour index. As a result, it is not necessarily 
true that a trace in (4.8)-(4.13) (which is a trace over fermion types) will translate into a 
trace over flavours in the MSSM. 

The process of converting the general RGEs to four-component notation involves 
significantly increasing the number of terms. Writing the Yukawa part of the general two- 
component Lagrangian as 

£ (2) 3 i^yD^ p + X -D^ a D^ a - \f^ a F% - QY^VjCV'^a + h.c.) , (4.14) 

the two-component RGE for the Yukawa couplings [13, 38] is* 

(4vr) 2 (3^ loop =\ [Y%(F)Y a + Y a Y 2 (F)] + 2Y 6 Y a tY 6 

r / m ( 4 - 15 ) 

+ Y b Tr | \ (Y bt Y a + Y at Y 6 J j - 3g 2 {C 2 (F), Y a } . 

where, Y 2 (F) = Y 6t Y 6 . Clearly, the matrices Y a , Y 2 (F) and C 2 (F) all have the same 
dimensionality, determined by the total number of two-component fermion fields in the sys- 
tem. Although the four-component RGEs of (4.8)-(4.13) look longer and more cumbersome, 
much work needed to obtain the MSSM RGEs has been done and our general equations 
therefore represent an intermediate step in this process. The reader is strongly encouraged 
to consult Ref. [39] , not least because it contains a detailed discussion of the method used 
to convert between notations. 

Before concluding this section, we note that we have not listed the RGEs for the 
dimensionless A a f, c rf and A' abcd , all of which (in the MSSM) correspond to SUSY counterparts 
of standard SM terms. We anticipate that threshold corrections to these terms will be 
unimportant for our analysis since they will not enter two-body (flavour-violating) squark 
decays except at the loop-level. In what follows we have set all these quartic couplings equal 
to the "square" of either the usual Yukawa couplings or the corresponding gauge couplings, 
their value in the SUSY limit. 

f Eq. (4.15) is slightly modified from that in Ref. [13]. We have written the Y^(F) instead of Yj(F) in 
the first term and symmetrized the trace with respect to a and b. The second modification also appears in 
Ref. [38] , while the first one preserves the symmetry of the Yukawa coupling matrix. 
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4.2.2 Dimensionful Parameters 

As with the dimensionless parameters, we use the general equations in Ref. [38] to 
derive the full threshold RGEs for the dimensionful parameters. In the general Lagrangian 
of (4.7) we write the fermion mass terms as m^, with the subscript 'X' indicating that they 
are mass terms for Majorana fermions. We do not write terms like m^, my or m[y since 
gauge invariance precludes the corresponding fermion bilinears in the -R-parity conserving 
MSSM. 

Applying a similar conversion from two- to four-component fermion spinors and 
from real to complex scalars, we find [40] (writing t = In Q) 

(47r)2 =\ [( WfeW » + vIv * b + X W + X W) ( mx ~ 

+ (m x + *m x ) (W b wl + v£V£ + XjxJ* + X 2 X & 2t ) 

+ (m x - im ' x ) (w;wZ + v*v b + xjtxj + x 2t x 2 ) 

+ (W 6 *W^ + VjV 6 + Xj'xJ + X 2t X 2 ) (m x + im! x ) 
+ [Xj (m x + iui' x ) X 2 + X 2 (m x + im! x ) Xj 

+X 2t (m x - zm' x ) Xj* + Xj* (m x - im' x ) X 2t ] (4.16) 
+ 1 [x^Tr {X^t ( mx _ im ' x ) + ( mx + im ^) X 2 } 
+ X 2 Tr {X 2t (mx - *m x ) + (m x + im' x ) Xj} 
+ X^Tr { (mx + im' x ) X, 1 + X 2t (mx - im' x ) } 
+X 2t Tr { (mx + im' x ) X 2 + Xj* (mx - im' x ) }] 
- 6 5 2 Cf (F)mx , 
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and 



W b wt + V^V 6 * + XjxJ* + X 2 X 2t ) (m x - zm' x ) 

- (m x + zm' x ) (W 6 WJ + V*V£ + XjxJ* + X 2 X & 2t ) 
+ (m x - *m x ) (W 6 *W^ + Vjv 6 + Xj'xJ + X 2t X 2 ) 

- (W 6 *W^ + VJV 6 + Xj*Xj + X 2t X 2 ) (mx + im! x ) 
+ i [Xj (mx + im' x ) X 2 . + X 2 . (m x + zm x ) X 6 X 



-X? (mx - im' x ) Xjt - Xjt (mx - im x ) X 
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+ 



(4.17) 



X*Tr {xj* (mx - im! x ) + (m x + im' x ) X 2 } 
+ X 2 Tr {x 2t (mx - im^) + (m x + im^) X*} 
- X^Tr { (mx + im' x ) X, 1 + X? (m x - im' x ) } 
-X 2t Tr { (mx + im' x ) X 2 + Xj* (mx - im' x ) } 
- 6 5 2 Cf (F)m x , 

where, as mentioned after the dimensionless RGEs, C\ l (F) is the quadratic Casimir for the 
Majorana fermions. Note that, as with the RGEs for xi' 2 , (4.16) and (4.17) are symmetric 
in line with the property of m^ in (4.7). Also, in the MSSM with i?-parity conservation, 
the trace terms in (4.16) and (4.17) vanish. This is because, as seen in Eq. (4.7), X only 
connects Higgsinos to gauginos, while m^ never connects Higgsinos to gauginos. Therefore, 
a trace of the multiple of these two matrices is always zero. 
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Writing the scalar quadratic Casimir as (^(S 1 ) = t A t A , the RGE for the trilinear 
couplings is 

(4?r)2 dtu, = [2 + A ^ R ^ + KhefU:ef + {b ^ c)] 

+ A e/bc H ae/ + 2A' efec/ (H* a/ + H /ae ) 

+ 2Tr { (m x - im! x ) [(v{ {u^Wj + V 6 X 2 t} + (x*t + X 2 t) X*X 2 t 

+ w^w^x^t + x 2t {w a wt + (xi + X 2 ) x c 2 t} 
+w 6 *uf v; + x 2t { v^v: + x£ (x^t + x 2 t) }) 

+ (6-c)]} 

+ 2Tr { (Wj (m x + im',) { V^U 2 * + X'W*} 

+ (Xi + X 2 ) (m x + im' x ) X^X 2 t + V b (m x + im' x ) W a U 2 t 
+ X, 1 (m x + zm' x ) {w a Wt + (X* + X 2 ) X 2 t} 
+ V 6 (m x + im x )Xivt 

+x\ (m x + im y {v^v: + Xj (Xjt + X 2 t) }) 

+ (&~c)} 

+ Tr {U 2 | (Ui + U 2 ) + (Ujt + U 2 t) Uj, + Vt V a , + 

+ 1 {x 2 ] (Xi + X 2 ) + (Xjt + X 2 t) Xi,}} H a , 5c 
+ Tr {ujjuj + U 2t U 2 + vl,V b + W+W 6 , 

+^{x^+X 2 tx 2 }}H af/c 

+ Tr ju*Uj + U 2t U*, + I {x*Xj + X*Xj,} } (H 6 , ac + H^,) 
+ Tr {u^Uj + U 2 tu 2 , + V+,V C + WjW c , 

+^{x^Xi+X 2 tx 2 }|H ac , fe 

+ Tr |u 2t Uj + U 2 tuJ, + 1 {x 2t Xj + X^X',}} (H c , a6 + H* bac ,) 
-3g 2 [C 2 {a) + C 2 {b) + C 2 (c)]U abc , 

(4.18) 

where, again, it can be seen that the b <-> c symmetry is preserved by the RGE as it must 
be. Note that in the .R-parity conserving MSSM we can choose (without loss of generality) 
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for the first index of H a b c to apply to a sfermion. This means that Ul' 2 and X„ 2 are always 
zero, and the corresponding terms drop out from the equation. 

Finally, since in the -R-parity conserving MSSM we never have non-zero entries in 
the Lagrangian corresponding to both m 2 ab and B a b for the same a, b, we write the RGE for 
a combination of these two terms. 



dt 



{2 (A afbe + A' fabe ) m 2 ef - KbefKf ~ (A e/ a6 + Kaef) &ef} 
+ {H ae/ H^ e/ + 2 (H* a/ + Hfae) H efe/ } 



2 Tr 



{ 



-2f 



v^v: + w a wt + (xi + x2) x 2 

+XJ (X^t + X 2t)] (m x - imy (m x + i m ' x )} 

+ Tr { (Xi + Xl) (m x + im' x ) X* (m x + im' x )} 
+Tr { (X^t + X 2 j) (m x - im' x ) X 2t (m x - im^) }] 
- 3<? 2 [C 2 (a) + C 2 (b)] (m 2 ah - B ab ) 

+ [ Tr {u? (ui + U 2 ) + (i# + u 2 t) ui, + vt V a , + wl W a 
+ix 2 t (xl + XD + 1 - (xit + x^t) xi,} ^ 

" Tr { U a' ( U « + U a) + + U 2 , 

+ vt ; v a + wtw a , 

+±x;j (xl + x2) + \ (x^t + x?) X 2 , } B a , h 

+ Tr{u 6 1 , t Ui+U 2t u2 +VJ,V 6 

+Wtw 6 , + ix^Xj + ^X 2t X 2 , | ( m ^ - Bav) 
+ Tr jujuj + U*Uj, + ±X*Xj + ^X*Xj,j 



(4.19) 

1 2 

In a similar manner to the RGE for H a bc above, when a and b are sfermion indices, U a ' b 
and X*' 2 are zero. Likewise, when a and b label Higgs fields, V aj b and W aj b are zero. The 
derivation of the MSSM RGEs is, therefore, considerably less cumbersome than it appears 
at first sight. 
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It is worth repeating before ending this section that we have written the RGEs 
for both the dimensionless and dimensionful couplings with regard only to the imparity 
conserving MSSM. As such, we have removed terms which are zero in this context. Complete 
RGEs derived from the Lagrangian in (4.7) can be found in Refs. [39] and [40] which also 
contain a more detailed description of the procedure for obtaining our general RGEs from 
the two-component ones. 



4.3 Choice of Renormalisation Scheme 

The removal of divergences as part of the renormalisation process requires us to 
choose a regularisation scheme that is consistent with all the symmetries in our theory. 
Although the modified minimal subtraction (MS) scheme is most commonly used in the 
Standard Model, the modified dimensional reduction (DR) scheme is more appropriate for 
supersymmetric theories. 

As we will describe in Ch. 6, we use the MS gauge couplings at the weak scale as 
our boundary conditions, and these will therefore need to be converted to DR, using the 
prescription in, for example, Ref. [41] : 

1 1 



a™ s a DR 



(4.20) 

l 

+ ^ (4.21) 



a MS a DR 

-L= =-L= + -J- . (4.22) 

a MS a DR 47T V > 

In addition to converting our weak scale inputs, we must also consider the effects of the 
relations (4.20)-(4.22) on the RGEs. We find that the additional contributions that appear 
when we convert the one-loop RGEs are of two-loop order and have the effect of changing 
the various coefficients in the two-loop RGEs. Therefore, since we wish to keep all two-loop 
sized effects, we must ensure that all RGEs we use are converted to the DR scheme. 

The invariance of the one-loop RGEs themselves under MS-DR conversion means 
that we are free to use our adaptation of the general RGEs from Refs. [13, 38] (which use 
the MS scheme) to derive new one-loop RGEs with threshold effects included. However, 
our two-loop RGEs are taken from existing literature, and the SM RGEs must therefore 
be converted. In addition to the gauge couplings, the Yukawa couplings also differ between 
the two scehemes. We use the symbol X Uj d,e (as opposed to f Uj d, e ) to differentiate the SM 
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Yukawa matrix from the MSSM one''' and convert the two-loop RGEs [17] using the following 
relations [41] 

A « " A « I 1 12016*2 8167T2 + 3167T2/ ^ 

\MS_-vDiif-i 13 9\ 3 g% 4 _gi \ 0/1 n 
~~ t 1 " 120 16^ "816^ + 3 16^/ ^ 

a ms _ A rm f 1 , A_g!_ _ 3_gl_l a 251 

Ae " Ae \ 1 + 40 16vr2 8 167T2/ ' (4 ' 25) 

where we omit scheme labels from the gj terms since we are only keeping two-loop sized 
differences. The RGEs for two-loop running of the gauge and Yukawa couplings are listed 
in Appendix C, and since the two-loop terms of the Yukawa couplings depend on the SM 
Higgs quartic coupling, A, its RGE is also given. We do not convert the A RGE to DR since 
its appearance at the two-loop level will mean corrections would only be of three-loop order. 
A simple check of the g% coefficient in the DR RGEs for Yukawa couplings is outlined in 
Appendix D, where we compare our result with that given in Ref. [42]. 

Having developed our one-loop RGEs with threshold corrections, we augment these 
with two-loop terms without threshold corrections from the RGEs for the SM [17] and the 
MSSM [15]. We change the two-loop RGEs from MSSM to SM at a single threshold, m H , 
so that our running consists of one-loop RGEs with full thresholds and two-loop DR RGEs 
with a single threshold. Threshold effects in the two- loop RGEs are expected to be of 
three-loop order and therefore small enough to be neglected in our analysis. 



4.4 The MSSM to SM Transition 

We now turn to consider the transition from SM to MSSM evolution. In Fig. 4.1 we 
plotted the SM running below 1 TeV and MSSM running above. This is a single 'threshold' 
where our effective theory transitions from being the SM, with SM field content, to being 
softly broken SUSY. The coincidence of unification coupled with the fact that SUSY is 
stable to radiative corrections strongly suggests that softly broken SUSY may be the correct 
effective theory all the way to ~ 10 16 GeV. 

Considering that we wish to integrate the RGEs to full two-loop accuracy, we are 

led to consider the errors introduced by our approximation of a single threshold. We first 

^The two Yukawa matrices are of course related, something we will discuss in Sec. 4.4.1 when we consider 
the Higgs boson mass eigenstates. 
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write the change in some quantity, q, due to the one-loop running over the whole range 
M z < Q < M GUT as 



\ M z J 



where (3 q represents the RGE for q, and then compare order of magnitude estimates of the 
corrections to this change from two- loop terms versus threshold corrections. 

(2) 

Two- loop: The two-loop contribution to the weak scale value of q, S q , is smaller than 
the one-loop contribution to the running, 8q~\ by approximately a factor of 1/167T 2 , 
so that the size of the two-loop contribution can be symbolically written as 

q 16vr 2 q 

Thresholds: Smeared out thresholds will only affect the running at most over the range 
from the heaviest SUSY particle (mnsp) to the lightest SUSY particle (tolsp)- In the 
intermediate range, the /3-function for the one-loop running will be neither equal to 
the MSSM value nor the SM value but will be the same order of magnitude. We can 
therefore estimate the correction to the one-loop weak scale value of q, 5q h \ to be 

~ /Vn ( ^) ~ 5?hn ( ^)/ln ( . 
y V m LSP/ y V m LSP/ \ M Z J 

We can see that if In (mnsp/m-LSp) ~ 1, which is a very likely scenario, these two contribu- 
tions will be comparable. We therefore find that we must include more detailed threshold 
effects to one- loop order, since these can easily be of the size of two-loop corrections, or 
even more sizable if the mass spectrum of SUSY is spread over a large range. 

One method to improve on the single threshold approximation is to include Higgs 
boson and SUSY thresholds as step functions in the RGEs, an approach which has been 
considered before [17, 18] and is also used here. We follow Ref. [17] and make the assumption 
that each particle with mass Mj is included in the effective theory if Q > M{ and excluded 
if Q < Mi. Using this approach, the /3-functions are as in the MSSM at a scale above 
the masses of all SUSY particles. As we move down in scale, the particles are decoupled 
individually until we have only SM particles in the theory and the /3-functions reduce to 
those of the SM. The reduction to the RGEs of the MSSM and SM in the appropriate limits 
will be a check on our work. 

Our RGEs must be able to describe the slow transition from the MSSM to the 
SM via a sequence of effective theories with changing numbers of particles. To this end, we 
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separate out the contribution to the RGEs from each particle and identify the terms with 
a 9-p for each particle V [17, 18], such that 

9 V = 1 if Q > M v 
if Q < M v . 

In this manner, we obtain the MSSM limit of our equations by setting all 9-p = 1 and the 
SM limit when all 9-p corresponding to decoupled SUSY particles are set equal to zero. 

If the separation between mass scales is not large, the errors introduced by decou- 
pling the various particles at the same point, say Msusy, are small. Since the threshold 
effects enter only logarithmically in the RGE solutions, we have seen that any error would 
be of O l^jln^ 1 ! and therefore less important as Mj — > Msusy- As a result, we 
see that it is permissible to find only approximate mass eigenstates, as long as we ensure 
that the size of the splitting between masses is small enough to introduce only a negligible 
amount of error compared to the two-loop corrections. 

Attempting to follow this route uncovers two issues: 

1. How do we define the 'mass eigenstates' that are decoupled? 

2. What happens once SUSY is explicitly broken by the decoupling of the heaviest SUSY 
particle? 

The following two sections consider each of these questions in turn. 

4.4.1 Mass Eigenstates 

Although it is possible to develop a scheme in which particles are all decoupled 
from the effective theory at a given point, this procedure involves its own complications. If, 
for example, we had run the gauge couplings in Fig. 4.1 all the way to Mz according to the 
MSSM RGEs, we would have to apply a correction to the weak scale couplings proportional 
to the logarithm of Msusy/Mz. In the more realistic scenario where the SUSY particles 
have a variety of masses in the range of 1 TeV, it would be necessary to apply a number 
of different corrections to the weak scale values. Also, the larger the difference between 
MgusY and Mz the less accurate this procedure becomes. 

In many cases it is therefore simplest to decouple the various SUSY particles at 
the scale of their mass. If we are to do this we must rewrite the Lagrangian in terms of 
fields which are in their (approximate) mass basis, so we continue by considering the mass 
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eigenstates of each sector of the MSSM in turn. The one exception to our general approach 
is the top quark, whose mass is ~ 172 GeV. Since we wish to maintain SU{2) invariance, 
and decoupling the top would break this, we decouple this quark at the scale Mz and 
apply a fixed correction to the gauge couplings and Yukawa matrices to compensate [43]. 
This allows us to keep our RGEs SU(2) invariant over the whole range we consider, i.e., 
M z <Q < M GU T- 

Higgs Sector 

In the Higgs sector the rotation to the mass eigenstate basis is particularly impor- 
tant. Since the MSSM contains two Higgs doublets, h u and ha, we have a number of mass 
eigenstates, labelled h, H, A, and H ± , which are combinations of the h u and fields. It 
is not clear, therefore, where we would place the transition from a theory including the 
field, to one where this field is absent. In contrast, it is clear where to decouple, for example, 
the A field, since it has a definite mass. Moreover, we can write our fields in such a way 
that the lightest field, h, does the job of the SM Higgs boson and remains in the theory all 
the way down to the weak scale. 

In view of the discussion just before the beginning of this section on the definition 
of the approximate mass eigenstates, we can simplify our treatment of the Higgs sector 
considerably. We find that threshold effects are only important if the mass scale of the non- 
SM Higgs bosons is much larger than m^. In this scenario, diagonalization of the various 
mass matrices in the Higgs sector leads us to conclude that — ran — m H ± >> m h with 
h approximately the SM Higgs particle. We rearrange the doublets h u and hd into two new 
doublets which transform as 2s under SU (2) and have positive weak hypercharge: 

' I k " ! d-20) 

-s\ kd |, (4.27) 
where the electrically neutral, complex fields h and TL are given by, 
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and s = sin (5 and c = cos [3. Here, G° and G + are the would-be Goldstone bosons which 
are incorporated into the longitudinal components of the heavy gauge bosons as a result of 
the Higgs mechanism. 

The heavy doublet (containing A, H and H ± ) will be decoupled at a common 
scale, ma- At this stage some terms in the Lagrangian drop out, and to continue running 
we must switch from considering the Yukawa matrices, f u ,d,ei which couple to both h u and 
hd, to running only that component of the Yukawa couplings which couples to the remaining 
field, h. A sample term in the MSSM Lagrangian which couples the up-type quarks to the 
Higgs doublets via the Yukawa matrices is 

C 3 -Uj{i u )%hlP LUi = -%(f u )£(sh + cH)P LUl . (4.29) 

We can see that at scales below run, the coefficient of the operator that remains in the 
theory is sf„, which we will identify as the Standard Model Yukawa coupling, denoted 
previously by A n . Similarly, the other SM Yukawa matrices are A^ = cf^ and A e = cf e . We 
will write our RGEs in such a way so that as we pass the threshold for the heavy Higgs 
particles we switch from evolving the MSSM Yukawa couplings (f u ,d,e) to evolving the SM 
Yukawa couplings (X u ,d,e)- 

Note that the rotation defined in (4.26) and (4.27) depends on a fixed angle, tan/? 
— the same angle as in the standard mSUGRA parameters. We therefore define the relation 
between \ u ,d,e and f Uj d,e to be valid only at the scale Q = ran, along with the mSUGRA 
input for tan/3. At this scale, tan/? = v u /vd- 

Our rotation to the approximate Higgs mass basis is unimportant when Q > 
itih since all we have done is make a field redefinition. The point is that by making the 
rotation we are able to identify the Higgs field that remains in the theory below ran- As a 
consequence, we are also able to determine the specific combinations of original Lagrangian 
parameters which remain in the effective theory. We will return to this point when we 
discuss the derivation of the MSSM RGEs in more detail. 

Neutralinos and Charginos 

Since decoupling effects are unimportant if the scale of the decoupling particle 
is of the order of Mz, threshold corrections due to charginos and neutralinos will only be 
important if either or the SSB gaugino masses are significantly larger than the weak scale. 
In this scenario, higgsinos and gauginos will have a small mixing and the mass eigenstates 
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are well approximated to be the bino, three winos and two higgsino doublets given by 

~ hl , 2 = . (4.30) 

Both higgsinos have the same mass, decoupling at the scale Q = \fi\, and therefore the 
rotation of the fields is largely irrelevant for this study. Having said this, it is important to 
note that the rotation that we perform on the higgsinos is not the same as that for their spin 
partners. This appears to be different from the approach used elsewhere in the literature 
[17, 18]. 

Finally, note that there is the potential for some of the fermion fields to have 
negative mass eigenvalues. In this case, we must redefine the fields [6] — ► («7s) 6 ' fe Vl/fc, 
where 6^ = 0(1) for positive (negative) eigenvalues, so that the mass term in the Lagrangian 
has the standard form. We have checked that the RGEs are independent of these additional 
75 factors. 

Squark (and Slepton) Sector 

We have already described in some detail the potential for mixing among squarks 
(and sleptons) of equal charge. Ignoring the small effect on the eigenstate from flavour- 
violating contributions, the left- and right-handed squarks will mix with each other so that 
the mass eigenstates of the up- type squarks, for example, are formed from ul — ur, &l — c~r 
or II — Ir. This mixing is intrinsically an SU{2) breaking effect, since the SU{2) doublet 
superfield Q contains both ul and di which will form parts of mass eigenstates with different 
mass. If we are to truly describe the effects from such mixing, we would have to introduce 
new SU(2) x U(l) violating operators that would appear when the squarks acquire mass, 
and derive corresponding RGEs for these terms. 

We argue, however, that to the extent that our calculation only depends logarith- 
mically on the location of the thresholds, we can ignore left-right mixing in the squark sector 
and decouple them at the scale of the eigenvalues of the respective SSB mass matrices. This 
approximation only breaks down when the off-diagonal entries in the squark mass matrices 
result in large cancellations, causing the physical mass of the lighter squark eigenvalue to 
become significantly smaller than either of the diagonal entries. We therefore do not claim 
to be accurately describing such a situation, which would require full consideration of the 
SU(2) breaking terms referred to above. 
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As with the other particles, we must identify the approximate mass scale at which 
to decouple each squark. Since we are neglecting SU(2) x U(l) breaking effects when 
pinpointing these thresholds, they will be located solely by the eigenvalues of the various 
SSB mass matrices. Approaching the highest squark threshold from above, we will reach 
a point where the largest eigenvalue becomes less than the scale, Qq. At this point, we 
decouple the corresponding squark, keeping all others in our effective theory below Qq. To 
implement the decoupling we rotate to the basis in which the corresponding SSB matrix is 
diagonal at Q = Qq using the appropriate one of 





u R = R u u|f , (4.31b) 

cIr = ~RddR , for squarks, or (4.31c) 
( e M \ 

Rl ^j, (4.31d) 

R e e^ for sleptons. (4.31e) 

The unitary rotation matrices, R,, are chosen to diagonalize the Hermitian SSB squark (or 
slepton) mass matrices, for example: 

(R^R Q )„ = «)* ag = {ml)p8 i3 . (4.32) 

With this in mind, the RGEs will be developed such that the contributions from each 
squark, q\ (with q = Ql, ur or cIr) with mass basis index k, will be identified with a 9q k . 
This identification requires that we write our RGEs in the current basis where the various 
SSB mass matrices are diagonal. In any other current basis, the 6g k become matrices in 
flavour space, 0^. Taking 6q^ for example, we can rotate back to the original current basis 
using 

(e Q )ij = (R ©J iag R^ 3 = 0Q fc (Ro) ifc a w (Ri)y (4.33) 

It should be clear that when Oq^ = 1 for all k, Gq = 11 so that in this regime we can use our 
RGEs in the current basis of our choosing without first rotating to the squark mass basis. 

Below the scale of decoupling of the heaviest squark, we continue to evolve the 
RGEs in our new current basis in which the heaviest squark decoupled. The eigenvector and 
mass of the decoupled squark are frozen and stored for use in connection with the physical 
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squark states, but the Oq k for the set of three squarks in question will ensure that the 
decoupled squark no longer contributes to the running. Although we continue calculating 
the running based upon three active squarks, the decoupled squark does not contribute, and 
the next squark will decouple at the heavier of the two eigenvalues of the (2 x 2) sub-matrix 
in the space orthogonal to the decoupled eigenvector. 

Finally, once two of the three squarks in the set have decoupled, the eigenvector of 
the remaining squark is fixed to be orthogonal to the other two. Thus only the eigenvalue 
evolves, and the squark is decoupled from the effective theory at the scale Q = Qo for 
this eigenvalue. Naturally, once all the squarks have decoupled, 6q k = for all q\ and 
equivalently all matrices & q = so that any rotation from the original current basis is no 
longer necessary. 

4.4.2 'Broken' SUSY 

The second question on our list of issues was that of the effect of explicitly breaking 
SUSY with the decoupling of the heaviest SUSY particle. Aside from the issues with 
changing operators in the Lagrangian as mentioned above (see the paragraph on page 49 at 
the end of the Higgs sector discussion), there is an additional effect which must be accounted 
for in a full two-loop analysis of the RGEs. This is the fact that when the heaviest particle 
is decoupled from our effective theory there is no longer any requirement that the couplings 
of SUSY particles remain the same as those of their SM counterparts. In fact, over the 
scale that they remain in the theory before they themselves decouple, their couplings have 
the potential to become significantly different, and through the RGEs appreciably affect 
the running of their SM counterparts. 

With this effect in mind, we must carry out our derivation of the RGEs allowing 
for differences in the RGEs for the SUSY couplings compared with those for their SM 
counterparts. In the SUSY limit, the scalar-fermion-fermion couplings of squarks and quarks 
to gauginos and higgsinos are equal to the gauge couplings and usual Yukawa matrices 
respectively. However, once the heaviest SUSY particle decouples, we find that the effective 
theory 'knows' about the SUSY breaking and as we shall see these SUSY 'Yukawa' couplings 
retain some contributions which are removed from their SM counterparts. 

For the gaugino couplings we must allow for flavour off-diagonal terms to develop, 
and we therefore write the couplings, which were diagonal and equal to the corresponding 
gauge coupling in the SUSY limit, as the matrices gf , where $ denotes the squark that 
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takes part in the interaction and i denotes the gauge group. Similarly, the higgsino couplings 
become different to their SM counterparts, f Uy d,e, and we label them f* de , again denoting 
the relevant scalar with a <£. 

The situation is best illustrated with an example. The coupling of a squark to a 
quark via an SU (2) gaugino is restricted in the MSSM to be equal to the quark-quark-W 
coupling. They are both proportional to the SU(2) gauge coupling, #2, with some fixed 
coefficient. In addition, the RGE (to one-loop) for §2 depends only on §2 and nothing else. 
However, there is nothing conceptually to prevent the squark-quark-gaugino coupling from 
being different to 02 once the effective theory contains broken SUSY. In this regime we 
therefore have a new coupling which we will label (following the discussion above) with a 
tilde, and also a superscript $ to signify the scalar field in the interaction. For the SU(2) 
coupling of left-handed squarks and quarks, we will call this coupling g^. 

In the notation for the additional coupling, we use a bold font since there is also 
the possibility of this coupling becoming a matrix in flavour space. The MSSM limit of 
this coupling is #2, which couples quarks and squarks of the same flavour. However, once 
we are in a broken SUSY regime, the RGEs allow for the possibility of flavour off-diagonal 
terms to appear. To see why, consider the RGE (A. 17) which will be derived in Ch. 5. 
The Feynman diagrams in Fig. 4.2 show all entries in the RGE that depend on the (flavour 
off-diagonal) Yukawa couplings. We find that these four contributions to the one loop RGE 
for g*^ cancel in the MSSM limit since they appear as 

RGE 3 (4.2a) + (4.26) - 4(4.2c) + 2(4.2d) = . 

However, if we are in a regime where our effective theory contains all the SUSY particles 
except higgsinos, we see that only (4.2a) remains in the RGE, which now contains flavour 
off-diagonal Yukawa coupling terms. We see similar features in the RGEs for the other 
sfermion-fermion-gaugino couplings . 

Off-diagonal terms in the squark-quark-gluino couplings induce flavour-violating 
decays of squarks, q\ — > qjZ^, even if the only source of flavour- violation in our theory is 
from the KM matrix via Yukawa couplings. Although the partial widths for these decays 
will be small, they may be important for light squarks for which standard tree-level decays, 
i.e., cji — > qjZi and % — > q'jWi, are kinematically forbidden. Therefore, the decay to a bino- 
like lightest neutralino, which is the LSP in many phenomenological models, is of particular 
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(a) 



(b) 




q 



q C3~ q 

q 




(c) 



(d) 



Figure 4.2: Feynman diagrams which contribute to the RGE for the squark-quark-gaugino coupling, 

5-Q 



interest. Although in this case it may seem that g' would be the most relevant coupling 
this will depend on the size of off-diagonal couplings from other sources. 

To illustrate the potential for sizable off-diagonal terms in these tilde-couplings, we 
preempt our discussion in Ch. 5-7 and examine the solution of the RGEs to two loop order 
using a simplified scenario for sparticle decoupling. In this scenario, we artificially place 
all thresholds at one of two points: the heavy Higgs bosons and gluinos at 2 TeV, and the 
sfermions and electroweak gauginos at 600 GeV. This helps to express the general features 
introduced by keeping the broken SUSY couplings separate from their SM counterparts. 
The theory we present is supersymmetric above 2 TeV, contains SM particles, squarks, 
sleptons, charginos and neutralinos above 600 GeV, and reduces to the SM below this scale. 

With these simplifications, the magnitudes of the couplings g /( 2 and g' UR are plot- 
ted in Fig. 4.3. The arrow at 2 TeV indicates the first threshold where the heavy Higgs 
particles and gluinos decouple from the theory. At this point, the boundary condition on 

-lib ■ 
g IS 

gi*(Q = 2 TeV) = g'Sij . 

The curves end at 600 GeV since this is the point where the squarks and gauginos decouple 
from the theory and the couplings no longer exist. As the vertical scale is logarithmic, and 
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Figure 4.3: The evolution of the magnitudes of the off-diagonal elements of g"2 and g,' UR for our 
simplified threshold scenario where the thresholds (denoted by arrows) are clustered at 2 TeV and 
600 GeV. Also shown for comparison is the running of the (2, 3) element of the up Yukawa coupling 
(see Sec. 7.2). Displayed in the "standard" current basis as discussed in Ch. 3, with the legend in 
the same order (top to bottom) as the curves. 

varies over a large range, the (2,3) element of the up-quark Yukawa coupling matrix f n 
is also plotted. t The (2,3) element is the largest of the off-diagonal entries of and we 
can see that there is potential for the off-diagonal entries of g to considerably exceed this 
coupling. It is therefore essential that the tilde-couplings are included in any discussion of 
flavour physics, especially when the KM matrix is the only source of flavour- violation via 
the Yukawa coupling matrices. 

Also seen in Fig. 4.3 is that the off-diagonal elements of the gaugino coupling 
matrices are symmetric under interchange of the two indices. This symmetry is not exact, 
but deviations are not visible on the plot. In fact, the elements in question are equal to 
within a few parts per mille. This approximate Hermiticity is observed in the gaugino 
coupling RGEs when the diagonal elements of the gaugino coupling matrices are large with 
respect to the off-diagonal elements, and the difference between the Higgs and higgsino 
coupling matrices, f u ,d, e and f* rfe , can be ignored. 

^Note that the features of the running of f u will be discussed in more detail in Ch. 7. 
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We move on to detail how the concepts in this chapter are applied to the MSSM 
and the way in which the subsequent RGEs are used to calculate the running for the 
couplings as in Fig. 4.3. After considering these points, we will return (in Ch. 7) to further 
examine specific examples of the running of the various parameters of the MSSM in both 
our simplified spectrum and more realistic scenarios. 
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Chapter 5 

Application to the MSSM 



The previous chapter lists the general RGEs for both the dimensionless and di- 
mensionful parameters in our general Lagrangian of Eq. (4.7). This general Lagrangian was 
constructed with a view to the i?-parity conserving MSSM which will be our effective theory 
between the highest SUSY threshold and Mqut- The next step is to write a Lagrangian, 
using (2.3) as a guide, for the specific particle content and gauge groups of the MSSM. 

Due to the splitting, below the highest SUSY threshold, between couplings that 
are otherwise equal in the MSSM, we must place different couplings for these terms in the 
Lagrangian on the understanding that they will be equal to their usual SM counterparts 
above the highest threshold. This requirement will serve as a boundary condition for some 
of the dimensionless SUSY couplings, as described later. In addition, keeping in mind that 
we wish to derive RGEs for the MSSM with full thresholds for SUSY particles, we must be 
careful to retain information about the source of each contribution by introducing a 9-p for 
each SUSY particle as described in Sec. 4.4. 

5.1 Interactions 

Before deriving the RGEs, we first write out all interaction terms in the imparity 
conserving MSSM applicable to our analysis, using the general supersymmetric Lagrangian 
in (2.3). We begin by obtaining the scalar-fermion-fermion interactions, which make up 
the U^' 2 , V$, W$ and X^' 2 matrices. The superpotential, (2.5), is used to evaluate the 
dimensionless couplings that govern the interactions of quarks and squarks with Higgs 
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(5.1) 



bosons and higgsinos, which are given by 

C 3 - [u^ u )lhlP LUi - UjiQlh+PLdi + d^lh-PLUi + d j (f d )T./ l 0p L d i 
+e j (f e )j i h-p L v i + e,(f e )J^P Lei + h.c] 

- u^lduP^rt + * h -/ Rj (f d d «)lP LUl + d^)] t u Lt P L * h - 

+e^)lD Ll P L ^ h - + * h0 e^<g*)lP L ei + ej(f e %e Li P L ^ h o + h.c. 

Notice that we write the interactions of quarks with the Higgs bosons as the usual superpo- 
tential Yukawa matrices. In contrast, the sfermion interactions have a tilde and are labelled 
with the sfermion that enters the interaction. Above all SUSY thresholds f* de = f u ,d,e, but 
below this scale, all the Yukawa' matrices may have different values. 

In addition, (2.3) describes scalar-fermion-fermion interactions involving gauginos, 
which are given by, 

(4,4.) GgpJ^+^gy G L P L 



V2 




+ Ujig^^ii-DumPLXo + djig'^jiiDdmPLXo 

+e,(g ,e «) jj (2)g ffi P L A + h.c.} (5-2) 

-v^{(gf%(-o* 5 4^¥ p £*-(^ 



h ~ d \ + li c 

"8 



V2 

where we have defined the following matrices: 

G = (g Q ^A 3 + |(g'^A (g«)ii (Ai - <A 2 ) 

Q (g Q ),* (Ai + 1A2) -(g%A 3 + i(g'%,Ao 

Gi = (g L )^A 3 - (g' L )^A (g L h (Ai - tA 2 ) 

(g L ), 4 (Ai + a 2 ) -(g L )^A 3 - (g ,L )^A 

^A 3 + ^"A ^(Ai-iA 2 ) , 

=| _ _ I ( 5 - 5 ) 

g h - (Ai + iA 2 ) — sr^Aa + g' hu Xo 
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G h d = 



(5.6) 



g hd (-\i + f\ 2 ) g hd ~\3-g' hd Xo 
and hypercharge and SU{2) gauginos are denoted Ao and Ai 5 2,3 respectively. Below the 
highest SUSY threshold, not only can the g* develop different values to their corresponding 
gauge boson couplings, but they are also able to obtain a non-trivial flavour structure, as 
pointed out in Sec. 4.4.2. In a similar manner to the higgsino couplings, above all SUSY 
thresholds the gaugino couplings must equal their gauge boson counterparts, so that the 
off-diagonal elements are zero, i.e., g'* = g' x 11, etc. 

The interaction terms for the dimensionful parameters originate mainly in the 
SSB part of the Lagrangian, (2.6). In addition to these terms there are mass terms for the 
higgsinos which arise from the superpotential via (2.3) 



C B - 
+ 



1 f 1 



(/x + n*) tt^o + ^o^o + + V h -* h± 



- (// - fj) [*fc0 75 * ft o + ^075^0 + ^+75% + ^ h -l^ h i 



(5.7) 



As mentioned previously, both higgsino mass states are taken to have a mass, |//|, and are 
therefore written with a single threshold, 6j L . 

Finally, there are there are additional terms derived from the following parts of 

(2.3): 



C 3 



2 4^ 



df 



dS r 



(5.8) 



s=s 



These result in scalar bilinear, trilinear and quartic terms. In particular, there are bilinear 
Higgs interactions with coupling strength \ arising from the superpotential, which is given 
a tilde to differentiate it from the n in the higgsino mass terms. As an illustration, if we 
take the second term in (5.8) and choose to differentiate with respect to the up-type Higgs 
superfield, i.e., S = H u , we see that 



(5.9) 



£ 3 - \fl\ 2 h°Jh° d - U ] Rk U ] Ll {i u )l n (f u )* m URmULn 

-(ui k (^y k u R ih° d +h.c) . 

In accord with our discussion at the end of Sec. 4.2.1, we have not placed additional labels 
on the Yukawa matrices that appear in the quartic part of (5.9) as we are not keeping 
track of the difference between this term and the "square" of the corresponding Yukawa 
coupling. Conversely, (/}*f^") fe/ is a trilinear coupling, and as such we will derive an RGE 
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for this term independently of li or f u (the usual up quark Yukawa coupling). This term 
enters the up-squark mass matrix as a left-right coupling term, and is therefore relevant for 
determining the squark mass eigenstates. Below the heaviest SUSY threshold, the coupling 
will evolve differently to the multiple (//* x f u ) due to differences between the respective 
RGEs, and in this case, where the term had its origin in the derivative of the superpotential 
with respect to H u we label the f u with h u . 



5.2 Gauge Coupling RGEs 



The one-loop RGEs for the gauge couplings including threshold effects are well 
known. We rewrite them here without quark thresholds since we will decouple the top 
at Q = Mz, and always work in the range Q > Mz- This means that we avoid the 
complications associated with SU(2) breaking, which was also our motivation for ignoring 
left-right mixing when deciding upon the location of the squark thresholds. 

Since there is no flavour structure in the one-loop gauge RGEs, the squark thresh- 
olds can be simplified using 



i=i 



Applying (4.2) to the MSSM particle content we obtain the familiar RGEs 



(4vr) 2 /3 Sl | 1 „ loop =5i 3 



A + *o N Q + ls N ** + h N ** 



~{9 h + e H ) + \e- h 



+ r N L + l N ** 



(5.10) 



22 1 1 1 2 4 

-Y + 4+ 2 N Q + 6 NL + 6 {9h + ° H) + 3 e ~ h + fw 



(4vr) 2 p. 



9a I 1 — loop 



-9! 



3 Q 6 R 6 dR 



+ 29- Q 
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(5.11) 



(5.12) 



Here, as in Sec. 4.1, g\ is the scaled hypercharge gauge coupling that unifies with the SU(2) 
and SU(3) couplings when the MSSM is embedded in a GUT: g' 2 = ^g 2 . It should be clear 
that the one-loop RGEs for the gauge couplings are a set of three differential equations, 
which form a closed set. As we introduce the RGEs for each group of couplings, we will 
increase the size of our closed set of differential equations until we have the whole system 
of RGEs for the MSSM. 
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5.3 Yukawa and Yukawa-type Couplings 



There are many RGEs for scalar-fermion-fermion couplings in the MSSM. Glancing 
at our Lagrangian terms for the softly broken MSSM, (5.1) and (5.2), we see that we must 
find RGEs for: the usual Yukawa couplings to Higgs bosons, f n , f<2, and f e ; the couplings 
of higgsinos to the various fermion-sfermion pairs, f-u •> 

if, f e L , fj", f e e «; hypercharge 
gaugino couplings, g /Q , g' L , g' UR , g' dR , g' eR , g' hu , g' hd ; the SU{2) gaugino couplings, g Q , g L , 
g hu , g hd ; and finally, the gluino couplings, gf , g" fl , g dR . The RGEs for all these couplings 
must be taken together, along with those for the usual gauge couplings in Sec. 5.2, to form 
a closed set. 

In order to apply the general RGEs to the particle content of the MSSM we must 
construct the various scalar-fermion-fermion coupling matrices from (4.7). Not only do we 
need to keep track of particle type, we must also remember the presence of three particle 
flavours for the Dirac fermions and sfermions, and it should be stressed that in this context, 
counting of different flavours is included in our count of particle types. We will write our 
final RGEs as matrices in flavour space, in accord with the notation we have used in our 
MSSM Lagrangian. 

We use a basis for the fermion content such that the Dirac fermions are 
{ui,di,Vi,ei} and the Majorana fermions are {h\, h^, hf, Ao, Ai, A2, A3, cja}, where i = 
1, . . .3 is the flavour index for the Dirac fermions, and h\ 2 an d hf 2 are the neutral and 
charged components of the Majorana spinor introduced in (4.30). In the MSSM, \j\ 
and U| will be (4 x 4) blocks of (3 x 3) matrices when $ is one of the Higgs bosons in (4.26) 
or (4.27). Similarly, since flavour is carried by the sfermion scalar index, $ = /j, Vj. will 
be a (4 x 9) matrix where the number of rows can be further expanded to show each of the 
three matter fermion flavours. Likewise, Wj. a (9 x 4) matrix where now the number of 
columns can be similarly expanded to exhibit these flavours. Thus, when fully written out, 
Vj. is a (12 x 9) matrix, while Wj. is a (9 x 12) matrix. Finally, X$ and X| will both be 
(9 x 9) matrices. In working out the size of these matrices, we have suppressed the colour 
index, A, which otherwise expands the size of the matrices due to additional quark entries. 

12 12 

Within the i?-parity conserving MSSM, the matrices U^' and X $ ' are non-zero 
only for = {b,7i,G + , H + } and the V$ and W<j> matrices are non-zero only for <I> = 
{uL,dL,e-L,VL,UR,dii,eR}. The matrices can be readily worked out from (5.1) and (5.2), 
and the derivation is tedious but straightforward. More detail can be found in Ref. [39] . 
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5.3.1 Comparison with Earlier Literature: Part I 

Before moving on to calculate the RGEs for the usual Yukawa couplings, we pause 
here to compare our method to that used in the previous literature. As these earlier studies 
did not keep the distinction between our tilde-couplings and the standard SM couplings, we 
perform an initial derivation for the up-type Yukawa coupling in the same manner. Having 
rotated the Higgs fields in the Lagrangian, the coupling of the up quarks to the light Higgs 
field, h, is s{ u and it is this RGE that we derive next. 

Note that while we find it instructive to write the RGE including separate thresh- 
olds for A, H, H^, and the four higgsinos, h\ 2 and hf 2 i we require that SU{2) remains 
unbroken. This corresponds to placing all the Higgs boson (and separately, higgsino) thresh- 
olds at the same scale, since it is only when we enforce this restriction that the couplings 
related by SU(2) have the same RGE, as we would expect. 



62 
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dt 



2 2 

y (o h + # G o + # G +) + j {Oh + e A + e H+ ) 



Ignoring the separation of tilde terms in the Lagrangian, we use (4.8) to obtain 
2 d(sf u )^ _3s 
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+53 2 {8^-4 



(5.13) 

In all cases of a repeated index, k, there is an implied sum, including the case of three 
repeated indices where one is a 9j k , as in, e.g., 9Q k (f u fl) ik (f u ) k y Since we have a collection 
of effective theories valid at various scales, with the SM and the MSSM as the limits at the 
weak scale and the high scale respectively, the RGE in (5.13) must reduce to the SM and 
the MSSM under the appropriate conditions. We will consider each case in turn: 

MSSM The MSSM limit is achieved when the effective theory contains all SUSY particles 
in addition to their SM counterparts. This corresponds to setting 0-p = 1 for all V . 
Since the rotation angle, tan {3, is fixed' , we can safely remove it from the differential 



^n the case that the rotation is scale-dependent, we see from (4.26) and (4.27) that the rotated fields 
themselves have an explicit scale dependence. Although our starting RGEs from Refs. [13, 38] do not take 
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on the left-hand side. This means that, since the right-hand side is proportional to 
sin P (in the limit we are considering) , the sin (5 cancels and we recover the equation 
for f„. 



SM The SM limit is achieved when the effective theory only contains SM particles (includ- 
ing, of course, the SM-like Higgs particle, h). To recover this limit, we therefore set 
all 6-p = except Oh, Oqo and 0q+. When we do this, we see that all dependence on 
the MSSM Yukawa coupling matrices, f u ,d,e, gives way to the SM Yukawa couplings, 
\ U} d,e, and we recover the SM RGE for X u [14, 44]. 

Before proceeding to our comparison with Refs. [17, 18], we note that if the higgsino 
threshold is at a lower scale than the heavy Higgs boson threshold, so that 0^ = 1 at the 
same time as Oh = 0, the RGE cannot be written only in terms of SM Yukawa couplings. 
The higgsino couplings, which we will soon write as f* de , remain in the theory and enter 
the RGE without a dependence on the Higgs rotation angle. 

Turning to our comparison with Ref. [17], we first note that our MSSM and SM 
limits agree, as expected. To compare the full RGE structure, we take a common threshold 
for the higgsinos and find agreement will all terms apart from the SU{2) and U{\) gauge 
couplings. We find, 

1 „ 4„ 
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with the down-type Yukawa coupling RGEs exhibiting a similar disagreement. 

In contrast to our method and that of Ref. [17], Ref. [18] appears to write the 
RGEs in a basis in which the Higgs fields are unrotated. It is therefore difficult to construct 



scale dependence of the fields into account, it should be clear that the rotation cannot change the RGEs 
above all thresholds since it is just a field redefinition. We can therefore safely conclude that the modification 
of the RGEs in the case of scale dependent fields would be precisely that necessary to cancel the contribution 
from removing the sin (3 from the differential on the left-hand side. 
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the appropriate relationship between their Higgs thresholds and the fields h, H, G + and 
H + . As a consequence it is hard to see how to reduce the RGEs in Ref. [18] to the SM, or to 
carry out a comparison with our result. The RGEs do, however, reduce to the MSSM when 
the appropriate limit is taken. We will return to compare our results for the dimensionful 
parameters to Ref. [18] in Sec. 5.7. 



5.3.2 Full Yukawa Coupling RGEs with Thresholds 

We must emphasise that (5.13) is not the correct RGE since it does not retain 
information about the splitting between the tilde-couplings and their usual SM counterparts. 
Using (4.8) to carry out the full derivation once more, retaining information about the 
tilde-couplings and setting single thresholds for both the heavy Higgs bosons, 9h, and the 
hie gsinos, 9^, the RGE for the coupling of the up-type quarks to the Higgs bosons becomes 

N 2 d(sf u )ij 
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(5.14) 
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In a similar manner to (5.13), three repeated indices indicate a sum, as in, for example, 

the term s(f u )jfc0^Qj(^)L(^)y - ^ n ^ s * s tne same f° r the RGEs in Appendix A where we 
repeat the i u RGE along with the RGEs for f^, f e , all their tilde counterparts, f, and the 
gaugino 'Yukawa' couplings, g. In addition, we repeat that the decoupling of the matter 
sfermions takes place in the (approximate) sfermion mass basis. The RGE is written with 
this in mind, i.e., using Of , so that, when the sfermions have partially decoupled, it should 

Ji 

be applied in the current basis where the specific sfermion soft mass matrix is diagonal.''' 

The reduction of the RGE in (5.14) to both the MSSM and SM limits is clear, with 
A n = si u the corresponding coupling below Q = iuh- If some SUSY particles remain in the 
theory after the heavy Higgs bosons have decoupled, we find that the factors of sin/? and 
cos/3 are exactly those needed to combine with the usual Yukawa matrices-'-, f u ,d,e, to form 
\i,d,e- It should be noted, however, that this is not true of the tilde-couplings, f, which 
would remain unchanged in this case. Again, the term s(f u )ik9j l 6 ] q (?u)\i(?u)ij illustrates 
this effect, since below m# we can write it as ( A u ) ^0^0^ (f„ )\i(?u 

Before moving on to consider the higgsino and gaugino couplings, we repeat that 
the reader should be vigilant when applying the general RGEs in (4.8)-(4.13). The matrix 
multiplication that takes place in these equations is over the full matrices which include all 
possible fermion types. In contrast, the matrix multiplication in (5.14) and in Appendices A 
and B takes place over the subspace of flavours. This means that not all trace terms in (4.8)- 
(4.13) lead to traces in the RGEs of the MSSM, and, conversely, not all traces in the RGEs 
are a result of traces in the general equations. For elaboration of this difference, see the 
comments below (36) of Ref. [39] . 

5.3.3 RGEs for Higgsino and Gaugino 'Yukawa' Couplings 

The RGEs for the couplings of squarks to quarks via gauginos and higgsinos are 
contained in the general RGEs of (4.12) and (4.13). In eq. (5.14) we allowed for the 
possibility that this coupling (and the other gaugino couplings to matter fermions) can 
develop a non-trivial flavour structure, as discussed in Sec. 4.4.2. This, of course, only 
happens at scales below the highest SUSY threshold, since above this scale we apply our 

*As stated in Sec. 4.4.1, when the three flavours of sfermions in a group are either all present or all absent 
from the theory, it is unnecessary to rotate to the applicable mass basis. This is because if all are present 
the rotation is just a field redefinition, and therefore cannot affect the RGEs, and if all are decoupled the 
whole term is absent from the RGE. 

*The reader will also note that there are additional factors of sin/3 and cos/3 which combine with Higgs- 
gaugino-higgsino couplings, as we will discuss in Sec. 5.3.3. 
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boundary condition: g* = gi x 11, where in this instance i indicates the various gauge 
groups. It is essential to keep track of such non-trivial flavour structure at the weak scale 
since these couplings are exactly the ones that enter the decays of squarks and sleptons. 

As an example, we write the RGE for the coupling of the right-handed up-type 
squark to up quarks and the £7(1) gaugino, g' u R; 



[s 2 e h + c 2 9 H ] (fZr u ) ik (g' u ») kj 
+ 



ki 



+ -^- lfc (g^)L(g /Ufl )^+^4(g /<iil )L(g^ 



+2^(g' efl ')L(g ,efl )^ 



+ ~ d Ms IUR h \ [*X + c 2 e H ] ~g' h - + [ c 2 e h + s 2 e H ] 



it 



-Jh n 



(5.15) 



3e s e- h [s 2 e h + c 2 e H ] ~g ,h HQl(^ R )l 3 



^h^Q k ^u)Jk { 



+ 0u k (s' UR )ik 

~4 



^g^)L(g' Ui %- + le~ 9 (g R )lM R h 



- (s IUR ) 



~,9 



/2 



4<?I 



As with the RGE for f u , care must be taken with the trace terms, which may not be 
immediately obvious in (5.15). An example of such a term is 



^B(g' UR h 



in which the trace is not explicitly written since we need to keep information about the 
position of the squark thresholds. Note that when we are checking reduction to the MSSM, 
the trace terms contribute an additional factor of three to account for the sum over flavours. 

In all the RGEs we take the approach that we do not insert explicit thresh- 
olds for external particles. This means that for the g' UR RGE we do not include a 0„ 
or 9yy when the term involves external forms of these particles. For example, the term 
-°h e Q k (^)J k (g' Q )ki(fu R )iji which comes from 2V 6 W*Wjf in (4.12), contains no fij , even 
though we have the coupling (f„ fl )^. When solving the RGEs, we do, however, freeze each 
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coupling at the scale where the heaviest particle decouples. Thus, for g /UR , if the wino is 
heavier than the up-type squark in question, the value of the coupling stops running when 
the wino decouples, and vice-versa. 

The RGEs for the couplings of higgsinos and gauginos to Higgs bosons, i.e., the 
interaction terms (5.5) and (5.6), are derived from the general expressions in (4.10) and 
(4.11). Since they are couplings with Higgs bosons, the rotation to the Higgs boson mass 
basis introduces a sin/3 or a cos (3 in the same way as the usual Yukawa couplings. As a 
result, we find that below mn, when the heavy Higgs fields decouple from our theory, we 
must use the combination of the coupling with the appropriate s or c, e.g., sg hu . Glancing 
back at (5.14) and (5.15) we see that the Higgs-gaugino-higgsino couplings always appear 
with the appropriate angle. This change to the subset of couplings that are associated with 
the Higgs field remaining in the effective theory is not only important for the Yukawa and 
ghu,h d couplings, but will appear again when we turn to discuss the trilinear couplings and 
Higgs boson scalar mass terms. 



5.4 RGEs for Gaugino Mass Terms 

To derive the RGEs for gaugino mass parameters we begin by constructing the 
matrices from the m['\ 3 and fj, terms appearing in our MSSM Lagrangian. The higgsino 
mass term, fi, should not be confused with the Higgs parameter, fx, which we will consider 
shortly. 

For the rotated Lagrangian, with the higgsinos defined as in (4.30), the matrices 
are diagonal and the derivation is remarkably simple. The RGE for M 2 , for example, 

is 



dt 



w 



3^(g«)«(g°)L + e Lk (s L )M(g L y ik + e- h ~ (s 2 e h + cH H ) 



+0l 



-Ji-d 



(c 2 e h + s 2 e H ) 



(5.16) 



2sc(-e h + H ) h ~g hd ^t u + (~9 hd Y >(^T - V9 w M 2 g 2 2 , 



r,h u \ 



where, once again, we write 6*'s for all internal particles. The 6-^ in the first and last terms 
may seem superfluous, but they are included to indicate to the reader that these terms are 
the result of loops which contain an internal W Majorana fermion. Note that the RGEs for 
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the gaugino and higgsino mass terms, when added to those for the gauge couplings and the 
full set of Yukawa couplings, once again form a closed set. 

The reader may be struck by the second term on the right-hand side, proportional 
to /x, whose appearance seems odd at first sight. It originates in the terms in (4.16) where 
is sandwiched between two X^' 2 matrices. Since X^' 2 connects gauginos to higgsinos, 
and the external fields are themselves gauginos, the is necessarily a higgsino mass 
term. Notice that this term is a threshold effect: it only enters the running of M2 below 
ma and only if the higgsinos decouple below this scale. The light and heavy Higgs boson 
doublets in (4.26) and (4.27), respectively, make equal but opposite contributions to the 
evolution of M2 which indeed cancel above all thresholds. The appearance of /x is a general 
feature of the electroweak gaugino mass RGEs — listed in full in Appendix B, (B.2)-(B.5). 
For much the same reasons, the RGE for the complex parameter fi, (B.l), develops a 
dependence on the electroweak gaugino mass parameters, M^p and if the masses are 

ordered appropriately. Finally, we point out that although the RGE in (5.16) develops 
a dependence on the complex parameter fx, it nevertheless preserves the reality of M2 as 
specified in the Lagrangian of Eq. (2.6). 

5.5 Trilinear Coupling RGEs 

In the i?-parity conserving MSSM, with interactions listed in (5.1)-(5.8), our col- 
lection of complex scalar fields consists of |h,H,G + , H + ,ULi,du,eu,^Li,URi^, where we 
have written the Higgs fields in their mass basis and the index i runs over all three flavours 
(the colour index is suppressed). We will use all these fields when constructing the the 
matrices H a6c , m 2 ab , &ab, A- a b C( i and A.' abcd . The resulting H a j, c is a (25 x 25 x 25) array, 
whose entries are mostly zero. 

Following on from our discussion regarding the Higgs rotation in Sec. 4.4.1, we 
must take a few moments to consider how the rotation determines the combinations of pa- 
rameters that remain in the theory below mu , when the heavy Higgs fields have decoupled. 
In addition to requiring us to redefine the quark Yukawa couplings (f u d,e ~^ ^u,d,e) an d 
the Higgs-gaugino-higgsino couplings, we also discover that we must restrict the trilinear 
couplings and Higgs boson mass terms. For example, when the two trilinear terms which 
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couple ul and ur to the Higgs bosons are rotated into the (h, W) basis, we see that 

T 



£ 3 u Rk (aI)klunh° u - u Rk (iffr 



kl 



uuh d 



u Rk {a T u )kiuu {sh + cH) - u Rk (^u u ) kl uli (ch - sH) , 



(5.17) 



which form the entries of the trilinear matrix^ H 



and H 



Upon decoupling 



the field 7i, we are left with only H^^^, and this part of the Lagrangian becomes 



C B u 



t 

Rk 



ki 



(5.18) 



Although in the complete MSSM (indeed, at all scales above win) we can sensibly talk 
about the evolution of the constituent pieces (a^)fcz and (ji*f^ u },, separately, for Q < ran 
it is no longer possible to write RGEs for both a u and /^*f„". We must instead talk only 
about the single combination ila Rk ,h,u L i that remains in the effective theory below Q = ma- 

Before using (4.18) to construct the RGE, we must also compose the quartic cou- 
pling matrices, A a b c d and A.' abcd , which are (25 x 25 x 25 x 25) arrays whose entries are 
given by (5.9). Our task is greatly simplified by noting that the A' a6cd -type couplings are 
absent above all SUSY thresholds in the i?-parity conserving MSSM since in this case the 
number of daggered and undaggered fields is always the same (remember that these origi- 
nate from the absolute square of a quadratic operator) . We find that A'-type couplings only 
arise because of the way we define the linear combinations of the Higgs rotations, (4.26) 
and (4.27), in order to take into account the decoupling of the heavy Higgs bosons. 

It is now straightforward to use (4.18) to derive the various trilinear RGEs, and 
for the coupling of up-type squarks to the light Higgs boson we obtain: 
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^It should be noted that we have chosen the index, a, of H a bc to be a sfermion as per the comments 
below (4.18). 
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(5.19) 



The RGE for the corresponding coupling, c(a u )jj + s (fi*i^ u )^ , to the heavier Higgs dou- 
blet can be obtained in the same manner. By taking linear combinations of these RGEs, we 
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can obtain the separate RGEst for a u and /2*f„". Since the coupling jl*f^ u always occurs 
as a product, it is not possible to obtain the RGEs for the individual factors. Of course, 
above all thresholds, we must have jl = \i and f£ u = f u . We have checked that with these 
replacements, our RGEs reduce to the MSSM RGEs [15] if we put all Oi = 1 and take care 
to sum over all internal flavours. 

Before closing this section we point out that in the -R-parity conserving MSSM 
the one-loop RGEs for the trilinear couplings do not contain dependence on the SSB scalar 
masses, a fact which is still true in the two-loop RGEs. This means that the complete set 
of coupled RGEs now includes the trilinear couplings, but the scalar mass RGEs are not 
needed except when setting the location of the various scalar thresholds. 

5.6 Scalar Mass RGEs 

Once we have derived the RGEs for the scalar mass parameters we will have 
completed the full set of RGEs for the -R-parity conserving MSSM (with the approximation 
that the A a b c d and A.' abcd couplings are equal to their usual SM counterparts). The scalar 
mass RGEs depend on all the interaction terms in our general Lagrangian, (4.7), including 
the matrices m a i, and lS a b- Most of these entries can be written down directly from the SSB 
Lagrangian, (2.6), but there are also |/i| 2 bilinears coming from the Lagrangian in Eq. (5.8) 
including the term written out in Eq. (5.9). Note that when the Lagrangian is written in 
terms of the rotated Higgs fields, there are no B a b terms, so that the RGE in (4.19) contains 
fewer terms. 

5.6.1 Mass Terms for Higgs Bosons 

As mentioned in connection with the trilinear RGE derivation, the mass terms for 

Higgs bosons in the Lagrangian become restricted when we are in the regime Q < run and 

the heavy Higgs fields have therefore decoupled. The complete set of mass terms in the 

^The s and c can be taken out of the derivatives on the left-hand side as noted in Sec. 5.3.1 when 
discussing reduction to the MSSM. 
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rotated Higgs basis are 



C 3 



s 2 (m 2 Hu + m 2 ) + c 2 (m 2 Hd + |/i| 2 ) - sc (b + b*) 



c 2 (m 2 + |A| 2 ) + s 2 (m 2 + |/z| 2 ) + sc(b + b*) 



sc 



sc 



c 2 b + s 2 b* 



(5.20) 



[ m H u + IAI 2 ) - sc (jn 2 H d + l/^l 2 ) - 
[m 2 Hu + |£| 2 ) - sc (m 2 Hd + |/i| 2 ) + s 2 b - c 2 b*] H*h . 
out of which only the first will remain in the effective theory below run- For Q > run, 
we can obtain RGEs for all four couplings and thereby derive separate RGEs for the real 
parameters {m 2 Hu + \ jl\ 2 ^j and (jn% d + |/2| 2 ^, and the complex parameter, b, by taking the 
appropriate linear combinations. 

We see from (5.20) that the term \jl\ 2 only appears in the Lagrangian in combi- 
nation with either m 2 Hu or m 2 Hd , so that we cannot derive an RGE for them separately. 
Naturally, above all thresholds supersymmetry requires that ji = //, so that we can use 
the RGE for fi, (B.l), to derive the RGEs for m 2 H ^ and m H , but these will no longer be 
valid below the scale of the heaviest SUSY particle. We can, however, use this to check the 
reduction of our RGEs to the MSSM limit when all 9i are set equal to unity. 

Applying (4.19) to the light Higgs boson soft mass term, we obtain the RGE: 
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(5.21) 

2 H ^ + |/i| 2 ^ and = (jnjj d + |A| 2 )- A discussion of the 

origin of each of the terms in this RGE can be found in Ref. [40] . It is perhaps informative 
to point out here, however, that all terms in (5.21) up to the trilinear terms derive from the 
single quartic, A a jbe m e/ m (4.19). All other quartic terms are zero either because A.' abcd 
vanishes when a = b or due to 3 a b having all vanishing entries in the Higgs mass basis. 
Also of note is that the terms in this group that are proportional to g' 2 contain the so-called 
S-term from the RGEs in Ref. [6, 15] and elsewhere. The term is noteworthy due to the 
fact that if it is zero at one scale (which is true at Mqut in many models) it will remain 
zero at least in the one-loop renormalisation. 
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5.6.2 Squark and Slepton SSB Mass Terms 

The RGEs for the sfermion soft mass terms are somewhat simpler on account of 
the fact that they have no Higgs fields in their operator in the Lagrangian. Otherwise, 
the derivation follows in a similar fashion to the mass terms for Higgs bosons above. We 
point out that all sfermion soft mass RGEs, (B.16)-(B.20) and (B.25)-(B.29), contain a 
dependence on the usual Yukawa matrices without a corresponding s or c to use in the 
conversion to X Ut d,e below m#. These terms have their origin in the quartic couplings, 
which we have set equal to the corresponding square of f u ,d,e, and would therefore have 
their running described by the quartic RGEs, not considered here. In solving the RGEs, 
since we do not have a valid RGE for f u d e , we use the value of the Yukawa couplings frozen 
at ran when Q < mn- 

5.7 Comparison with Earlier Literature: Part II 

As a short conclusion to this chapter we note that Ref. [17] only includes RGEs 
with thresholds for the Yukawa couplings, f u de- It is possible, however, to carry out a final 
comparison with Ref. [18], which writes out the RGEs for the Higgs boson, higgsino and 
gaugino mass terms, and the (3, 3) component of all the trilinear couplings and SSB scalar 
masses. 

We have checked that, except for terms involving couplings of Higgs boson fields, 
the RGEs that we obtain agree with those in Ref. [18]. We found it difficult to compare 
contributions involving Higgs boson fields since, as noted in Sec. 5.3.1, the RGEs of Ref. [18] 
appear to have been written without any rotation of these fields, making it impossible to 
abstract the relationship between the thresholds for the Higgs bosons h, A, H and in 
Ref. [18] and the thresholds in the RGEs listed in Appendix A and Appendix B. 
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Chapter 6 



Deriving the Weak Scale Couplings 
from the RGEs: RGEFLAV 

Now that we have the one-loop RGEs for the MSSM with full threshold effects from 
Chapter 5, we are in a position to combine these with the two- loop RGEs for the MSSM 
[15] and the SM [14] converted to the DR-scheme as in Appendix C. As described at the 
end of Sec. 4.3, our two-loop RGEs do not contain (numerically higher order) threshold 
effects, and, for the gauge and Yukawa couplings, we use the MSSM form above Q = mu 
and the SM form below Q = run- We solve this system of equations given our weak scale 
Standard Model inputs and GUT scale ansatz using FORTRAN computer code to iteratively 
reach a numerical solution. 

The code has been designed to be incorporated into the public release software, 
ISA JET [45], which is a Monte Carlo program that simulates pp, pp, and e + e~ interactions 
at high energies. RGEFLAV is the main program, which sets up the problem and calls the 
subroutines necessary to iteratively solve the RGEs. It will be employed as a subroutine of 
ISASUGRA, a section of ISAJET that can currently be used to solve the flavour-diagonal RGEs 
with a variety of different boundary conditions and then, through a call to an additional 
subroutine, ISASUSY, calculate the decay modes. 

This chapter aims to provide an overall description of the code and the various 
issues encountered in the procedure. After outlining the broad approach, we will consider 
each segment of the code individually, from the input file, through the choice and application 
of the various boundary conditions, to the final output. 
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General Outline 



Since RGEFLAV has been designed to be used in combination with ISASUGRA, this 
programme must be initialised first with equivalent mSUGRA inputs, and RGEFLAV is called 
at the end of the ISASUGRA code. A general outline of the order in which RGEFLAV carries 
out the various steps is: 

1. Read the input file INRGE.DAT that contains all the choices available to the user. 

2. Introduce gauge couplings, light quark Yukawa couplings, third generation Yukawa 
couplings from ISAJET, and the SM Higgs field VEV all at M z . 

3. Run the gauge couplings and diagonal Yukawa matrices to mt. At this scale, insert 
ISAJET's value for ft(m t ) and rotate the Yukawa coupling matrices into the current 
basis using user defined rotations. 

4. Run up to the high scale using appropriate thresholds, whose initial values are ob- 
tained from ISAJET, to transform SM running into MSSM running. 

5. Insert high scale boundary conditions. It is after this step that we begin the iterative 
loop. 

6. Run back down to the weak scale, decoupling particles as their thresholds are passed. 

7. At ran, change over to the new set of RGEs, with restricted Lagrangian terms, as 
a result of the decoupling of the heavy Higgs boson terms. Save the values of the 
couplings at mu for use when running back up. 

8. At MsusY) defined as yJfni L m t R > a PPly the electroweak symmetry breaking conditions 
and radiative corrections to the third generation Yukawa couplings [46] . 

9. Once the running reaches m t rotate back to the basis in which the Yukawa matrices 
are diagonal, reset ft(nit) to the value obtained from ISAJET and continue running to 
M z ■ From now on ft remains in the theory all the way to M z and will therefore change 
the running of the various couplings from that obtained during the first upwards run. 

10. Reset the values of the gauge and all Yukawa couplings other than ft. 

11. Run back up to the high scale, applying the rotation (and boundary condition on f t ) 
at mt and the thresholds as before. 
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12. Reset the high scale values of the a-parameters and SSB mass matrices. Iterate until 
a solution with pre-selected accuracy is found. 

13. On the final run, stop at mn, the heavy Higgs threshold, and output the couplings at 
this scale. 

14. At present, the code then calls STDEC, which is discussed in Sec. 8.4. 

The code must run as part of the ISAJET distribution because there are some common 
blocks which must be filled for RGEFLAV to execute properly. These common blocks are: 

• SUGPAS and SUGMG: Contain the initial values of the various thresholds. 

• WKYUK: Contains the weak scale third generation Yukawa couplings. 

• BSG: Contains the radiative corrections to the Yukawa couplings. 

• RGEMS: Contains the value of MgusY from ISAJET and the value of fi at this scale. 

The subtlety in initialising ISAJET is that equivalent mSUGRA inputs are needed 
since the interface has been developed to work with the mSUGRA-based executable ISA- 
SUGRA. The user may use non-universal GUT scale boundary conditions according to the 
choices given by ISASUGRA and available in the input file (discussed next), but the user 
must also enter approximate mSUGRA conditions before RGEFLAV can begin the iterative 
process. 

Once ISASUGRA has called RGEFLAV the signs of the numerical values of Aq and /U 
are flipped. This is because RGEFLAV has been developed for inputs that use the notation 
found in Ref. [6], and the sign changes are required to fix notational differences. Some signs 
are again flipped in the stop decay routine (STDEC) and this will be discussed later. 

6.1 The Input File 

Once RGEFLAV is called at the end of ISAJET, it reads the input from INRGE.DAT, 
a sample of which is shown in Appendix E. The choices, numbered in the same order as 
the input file, are: 

1. Whether to use the complete two-loop equations. If the following line is '0', only one- 
loop RGEs will be calculated, but threshold effects will still be included. Choosing to 
only calculate one- loop running speeds up the programme by around 15%. 
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2. This item allows the user to choose between complex inputs and real inputs. If the 
next line is '1', all inputs must be in complex form as in Appendix E and the KM 
matrix will include a phase. If '0' is entered, any optional inputs which follow must 
be single numbers representing one real entry each. Otherwise the programme will 
print an error message and continue no further. Choosing to use a real KM matrix 
and real boundary conditions dramatically speeds up the running of the RGEs. 

3. When the answer to Question 2. is '1' the user may wish to enter a phase for [i as 
opposed to a simple sign. If the user answers yes, the optional setting for THETA on 
the next line is used, otherwise this value is ignored, and a sign for \x consistent with 
the ISAJET input is used. 

4. If the response to Question 4. is '1' the programme will not attempt to change the 
thresholds from those obtained by ISAJET. Otherwise, thresholds for: the left- and 
right-handed squarks and sleptons (i.e., the eigenvalues of the SSB sfermion mass 
matrices), the higgsino (fi), the bino (Mi) and the winos (M2), will be set during the 
running. 

5. The mSUGRA GUT scale inputs used by ISAJET are automatically passed to RGE- 
FLAV, but they are only all used if the user sets this entry to '1'. If so, reading of the 
input file moves to Question 17. Otherwise, reading continues to the non-universal 
inputs, 6. — 16. 

Inputs 6. — 16. are the alternative GUT scale inputs. We will deal with these inputs 
in more detail in Sec. 6.4. After taking care of the GUT scale inputs, the file continues with 
the choices relating to the weak scale rotations Vl,r(u, d). 

17. The first choice in the final section is the basis in which to output the results. Since 
SU(2) is broken in the quark mass basis, we only use this basis for the weak scale 
inputs. Instead, the output is in a current basis in which either the up- or down- 
type quarks are diagonal at m t , thereby retaining SU(2) invariance. Note that when 
interfacing with STDEC, described in Sec. 8.4, the answer to Question 17. must be '1' 
since we wish to associate the cl,r squarks with the c quark mass eigenstate. 

18. This question asks whether the user should enter their own choice of general rotation 
matrices, as described in Sec. 6.2. These matrices will be used to rotate the mass 
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basis quark Yukawa matrices into the current basis and therefore this corresponds to 
a choice of current basis. Note, however, that the output is always in a current basis 
where either the up- or down-type quarks are diagonal at m t . If the user chooses '1', 
the inputs for Question 21. will be read next. If not, the file proceeds to the next 
question, which allows for more basic choices for the rotation matrices. 

19. Here the user chooses the form of Vj,(u). Either they can choose to have V l{u) = K, 
or V/,(u) = 11. The programme will ensure that \ l{<1) is fixed so that the appropriate 
combination of these two matrices is the correct KM matrix, according to (3.6). 

20. If the KM matrix is the only source of flavour- violation, the form of the Vr(u, d) 
matrices is unimportant. In this entry the user is given a restricted choice for the 
matrices Vr(u,cI) which serve as default entries if the response to Question 18. is '0'. 
Either both V^m) and Vr(cT) are the unit matrix or Vr(u) = and V^(cT) = K. 

21. The final entries are the choice of parameters which will be used to define the rotation 
matrices if the response to Question 18. is '1'. If the user has chosen to have real 
running in Question 2. the phase in these inputs is ignored. For more detail on the 
way these inputs are converted into unitary matrices, see Sec. 6.2. 

We now move on to consider some of the above inputs in more detail. 

6.2 Entering a General Unitary Matrix 

As mentioned above, Question 21. of INRGE.DAT allows the user to define their 
own Vl 5 #(u, d) matrices. Since the KM matrix must still be correct, Vi(ti), V#(u), Vr(cT) 
are taken in as inputs and V^(d) is then fixed by requiring that (3.6) is satisfied. 

In order to guarantee the unitary nature of the rotation''', the input file reads three 
angles (a, (3 and 7) and a phase (5) using a similar parametrisation to that used for the 

^The matrices ~Vl,r(u) should be numerically unitary to high accuracy. Otherwise, inverting the up 
Yukawa couplings from this new current basis to our "standard" basis will leave residual off-diagonal elements 
rather than zero even at Q — m t . If the size of these elements is comparable to the values of the smallest off- 
diagonal elements at values of Q substantially away from mt, it is clear that our solutions will be dominated 
by the error from the non-unitarity of the ~Vl,r(u) matrices. A similar consideration applies to Vr(cI). 



81 



KM matrix as given in (3.7). Each unitary matrix, U, is therefore given by 



U 



/ 1 \ / cp spe~ iS \ ( 

c 7 s 7 10 

\ -s 1 Of J \ spe iS cp J \ 



s« 0\ 
c a 
1 / 



(6.1) 



\ -spe 1 " U cp 

where s a = sin a, c a = cos a, etc. Note that this is not the most general unitary matrix 
possible, which would include an additional two phases. However, it is considered that for 
the time being this choice of inputs will suffice for most practical applications, and the 
addition of two phases would require only a minor change to the code. 



6.3 Weak Scale Boundary Conditions 

Once ISAJET calls the main program, and the input file has been read, RGE- 
FLAV fixes the weak scale boundary conditions and runs them up to the GUT scale using 
the subroutine UPMZMHIGH. For the gauge sector, we take as our input the current PDG 
values [34] 

a'^Mz) = 127.925 ±0.016 ; a s (M z ,MS) = 0.1176 ±0.002 ; 
sin 2 e w (M z , MS) = 0.23119 ± 0.00014 . 



These are the couplings extracted using the effective theory with the electroweak gauge 
bosons and the top quark integrated out at Q = Mz- In order to use the SM evolution 
for Q > Mz, we must match these couplings to those of the full SM, which to two-loop 
accuracy implies that the SM gauge couplings in the MS scheme are given by [43], 



where 
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Notice that in order to preserve the SU(2) symmetry of the effective theory down to Q = 
Mz, we have, as mentioned in Sec. 4.4.1, integrated out the top quark at Q = Mz rather 
than at its mass as we do for all other particles. This is the origin of the ln(m t /fi) terms 
in the matching conditions for the gauge couplings above. Since we decouple all SUSY 
particles as well as the additional Higgs bosons at the scale of their mass, we do not get 
corresponding jumps in the gauge couplings as these decouple. Our method — which is 
also used in ISAJET — has an important advantage that it "sums the logs of the ratio 
of any large sparticle or Higgs boson mass to Mz" , in contrast to the frequently used 
procedure that applies MSSM evolution all the way down to Mz, and then corrects for this 
via a "single-step evolution" (between the heavy scale and Mz) to take into account the 
difference between the running in the MSSM and in the SM. 

Next, we convert the values of these gauge couplings in the MS scheme to their 
corresponding values in the DR scheme using the relations (4.20)-(4.22) and use the results 
as boundary conditions at Q = Mz when solving the RGEs. 

For the Yukawa couplings, we begin with the quark masses at Q = Mz (the masses 
of the light quarks and leptons at Mz can be found in Ref. [47]), and convert to SM Yukawa 
coupling matrices using vsm = 248.6/\/2 GeV as in Ref. [46]. The masses of the first two 
generations of quarks have substantial error, which leads to a corresponding error in the 
Yukawa coupling. The third generation quark masses are more precisely known — we take 
the top pole mass nit = 172 GeV — and in practice the values of the third generation 
Yukawa couplings are taken from ISAJET, with bottom and tau Yukawa couplings at the 
scale Q = Mz and the top Yukawa coupling at Q = nit. In extracting these Yukawa 
couplings we include SUSY radiative corrections [46] at Msusy obtained by ISAJET during 
its execution, with inter-generation quark mixing neglected. 

These diagonal Yukawa couplings, which are in the "quark mass basis" are run to 
nit with no flavour structure since it is a good approximation that the running is mainly 
due to the strong coupling. At nit ah three Yukawa matrices are rotated to the user's choice 
of current basis using the SM version of (3.2) and the corresponding relation for A<f. 

Running then continues to the GUT scale with a basic RGE subroutine, RGE215, 
which only contains the RGEs necessary for running the gauge and Yukawa couplings, and 
implements only rudimentary thresholds. When we reach the scale Q = ran m the course 
of running up, we switch from SM Yukawa matrices (X u ,d,e) to MSSM Yukawa matrices 
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(fu,d,e) an d from the SM VEV (vsm) to the VEVs in the two Higgs doublet model: 



VSM = 




The subroutine HIGHIN takes care of the boundary conditions at the high scale as 
discussed next. 

6.4 Boundary Conditions at the High Scale 

The running is deemed to have reached the GUT scale when a\{Q)— a!2(Q) becomes 
negative, unless the user responded to Question 6. with a fixed Mhigh, in which case the 
running terminates at the value that was chosen in INRGE . DAT. 

Since our purpose is to simulate flavour physics of sparticles in as general a way as 
possible, subject to experimental constraints that seem to suggest that that flavour physics 
is largely restricted by the structure of the Yukawa coupling matrices, we thought it would 
be useful to first seek a general parametrisation for SSB parameters that does not introduce 
a new source of flavour-violation, but allows for non-universality of model parameters. A 
different source of flavour-violation can easily be incorporated by allowing for additional, 
arbitrary contributions to the SSB mass and trilinear parameter matrices. We use the 
(s)quark sector to illustrate our arguments, but almost identical considerations will apply 
to (s)leptons (except that in this case we would also have to include lepton number and 
lepton flavour- violating matrices in the singlet (s)neutrino sector). 

Within the framework of the i?-parity conserving MSSM, the SSB matrices m^-^ q 
and a Ut d potentially include new sources of flavour-violation, not included in the superpo- 
tential Yukawa couplings. In order not to introduce a new source of flavour- violation, these 
SSB matrices must be diagonal in the same superfield basis (where the SM fermions and 
their scalar superpartners are rotated by the same matrices) that the superpotential Yukawa 
interactions (renormalised at the same high scale as the SSB parameters) are diagonal. t Of 
course, it is impossible to simultaneously diagonalize f u and f^, but what we mean is that 
the SSB mass matrices that describe the mixing of both left- and right-handed up type 
squarks, and their trilinear couplings must be diagonal in the basis that the matrix f u is 
diagonal, and likewise for the down sector. However, since SU (2) symmetry dictates that 

^This is not equivalent to the requirement that the Yukawa coupling matrix commute with the corre- 
sponding a-parameter matrix because these non-Hermitian matrices are diagonalised by bi-unitary and not 
by unitary transformations. 
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the m~ and rrij SSB matrices must be identical, this matrix must be proportional to the 
unit matrix (11) in order to remain diagonal, both when the up- or the down- type Yukawa 
coupling matrix is diagonal. In contrast, the matrices mfj and a u (ui 2 D and a^) can be func- 
tions of the up (down) type Yukawa coupling matrices (and their Hermitian adjoints) chosen 
in such a way that these matrices are simultaneously diagonal when we transform these to 
the basis where the corresponding superpotential Yukawa coupling matrix is diagonal. 

To find the most general parametrisation of the mfj D and a. u ^ matrices of the type 
that we are looking for, we first note from (3.8)-(3.10) that under field "rotations", these, 



where n is any integer. Thus any linear combination of these matrices (with n = 0, 1, 2, ...) 
is guaranteed to be diagonal in the basis that f u ,d is diagonal (at the high scale at which we 
enter the SSB parameters). The only question, then, is just how many terms we need to 
allow in the linear combination to guarantee the most general form for the SSB matrices, 
so that flavour- violation enters only through the superpotential Yukawa coupling matrices. 
This is easiest to see in the diagonal basis for the Yukawa couplings. The SSB matrices are 
also diagonal in this basis, and so are completely specified by n g diagonal elements, where 
n g is the number of generations. Transforming to a general basis does not alter the number 
of parameters that we need: we thus know that we must have n g terms in each of the linear 
combinations for mfj D and for &u,d that we discussed above. For the MSSM with n g = 3 
generations, we thus parametrize the SSB sfermion mass and a-parameter matrices at the 
high scale as, 



where cjj,d,e = or 1, and f u ,d,e are the superpotential Yukawa coupling matrices in an arbi- 
trary current basis at the same scale at which the SSB parameters of the model are specified. 
The matrices Tq : l,u,d,e and Z Uj d,e have been introduced only to allow for additional sources 
of flavour-violation not contained in the Yukawa couplings. Setting Tq,l,u,d,e = ^u,d,e = 
gives us the most general parametrisation of the three- generation .R-parity conserving MSSM 
where the Yukawa coupling matrices are the sole source of flavour-violation. 

Note that this collection of inputs (excluding the matrices Tq,l,u,d,e and Z Uj d,e) 
is a special case of the minimal flavour violation principle [48] in that we assume the GUT 



respectively, transform in the same way as the matrices 




m U,D,E 
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) 2 ] + Tu,d,e (6.4b) 
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scale SSB parameters do not introduce new sources of flavour- violation. Since our conditions 
are not invariant under renormalisation group evolution, the weak scale form of the SSB 
parameters will not be the same as in (6.4), but will, however, be described by the more 
general minimal flavour- violation ansatz. 

Questions 9. to 16. in INRGE.DAT allow the user to choose arbitrary values of all 
the input coefficients above, subject to the constraint that Tq,l,u,d,e are Hermitian. The 
familiar universal mSUGRA boundary conditions are reproduced by setting cu,d,e = 1; 

m {Q,L}0 = m lu,D,E}0 = m h A {u,d,e}0 = A 0] Ru,D,E = ^U,D,E = W uAe = X Ujdfi = 0; 

Tq,l = Tu,d,e = Z„,d,e = in (6.4). 

The remaining GUT scale inputs — namely, those for the gaugino and Higgs boson 
scalar masses — are simple numbers, which can be either given by the mSUGRA parameters 
passed from ISAJET, or entered by the user in Questions. 7. and 8. of INRGE.DAT. 

6.5 Electroweak Symmetry Breaking 

After fixing the high scale parameters, the programme proceeds to the subrou- 
tine DOWNMHIGHMZ, which runs the entire collection of RGEs contained in the subroutine 
RGE646. At each step, the code checks to see whether the scale is below Msusy and, at the 
point Msusy is passed, applies the electroweak breaking conditions using the subroutine 
D0WNMSC0ND. 

As discussed in Sec. 2.3.4, it is common practice to use the observed value of M§ 
to determine the value of via (2.9b), but our inclusion of threshold effects comes with 
additional complications in this case. Since Msusy is at the scale y/ m t L m i R i t ne electroweak 
symmetry breaking conditions will always be applied at a scale smaller than the mass of the 
heaviest SUSY particle. As a result, fi, the higgsino mass parameter, is no longer equal to jl, 
the parameter that enters the Higgs potential. Moreover, the Higgs potential depends only 
on M 2 ^ = {rn?H u + |/^| 2 ^ and M^ d = {jn\j d + |/i| 2 ^, so that it is not possible to separate 
\jl\ 2 from the SSB parameters rnjj^ and fn"jj d that we specify at the high scale. Notice, 
however, that we can define the relations 

{M 2 Hu + M 2 Hd ) = m 2 Hu +m 2 Hd +2 |/i| 2 and 
(M 2 Hd -M 2 Hu )=m 2 Hd -m 2 Hu . 



86 



so that the tree-level minimisation conditions of the Higgs potential can be written as, 

b = sc{M 2 Hu + M 2 Hd ) (6.5) 

«. + M h) = " " M ^ " \ {9 ' 2 + ^ ^ + ^ ■ (6 - 6) 

The second of these fixes the sum (M 2 ^ + Mjy ) in terms of the difference. Since we know 
this difference at the high scale, we can evolve this down to Msusy (along with other SSB 
parameters) during the iterative process that we use to solve the RGEs. At Q = Msusy we 
use (6.6) to solve for (M 2 ^ + M# ), which can be evolved back to the high scale. We then 
use the sum to fix ft = ji at the high scale, reset the difference, (M 2 ^ + Mj^), to its input 
value, and iterate as per the discussion in Sec. 6.6. The value of the higgsino parameter ji 
can then be obtained at all scales using (B.l). 

Finally, with a small adjustment to the usual relation in Eq. (2.9a), the 6-parameter 
can be eliminated in favour of tan/3 using (6.5). Although the 6-parameter is complex in 
general, our decision to make v u and Vd real and positive^" requires that 6 also be real and 
positive at the scale Msusy- Our parameter does, however, retain the ability to develop 
complex parts as a result of the running, and will not necessarily remain real at all scales. 

Up to this point, we have ignored another potential complication that arises if 
Msusy < tuh- In this case, the heavy particles of the Higgs sector have decoupled by the 
time we apply the electroweak symmetry breaking conditions, and we only have the light 
doublet in the effective theory that we use to calculate the RGEs. In this case, the heavy 



Higgs doublet mass term 
terms 



c 2 ( m 2 + |/I| 2 ) + s 2 (m 2 + \tf )+sc(b + b* 



and the mixing 



2 Hu + |/2| 2 ) - sc (m 2 Hd + \fi\ 2 ) + s 2 b - c 2 b* 



(and its complex conjugate) from 



Eq. (5.20), together with tan/3, are frozen at their values at Q = run, while the light dou- 



blet mass parameter, 
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H u + W) +c 2 {mj Hd + \tf)-sc{b + b* 



, along with i>sm, 



continue to evolve to Msusy- The three frozen coefficients together with the evolved mass 
term for the light doublet must therefore be used to solve for {jn 2 }j d + |/^| 2 ^) {j^Hu + \^ 
and the complex 6-parameter. We can then find a solution in the same manner as for 
Msusy > m H - 

Before closing this section, we should add that although we have discussed EWSB 

conditions only at tree-level, in practice we minimize the one-loop effective potential in- 

can always make a gauge transformation such that just the lower component of H u has a VEV, and 
that this VEV is real and positive. Then the minimization of the scalar potential in the Higgs sector requires 
that the VEV of Ha is aligned; i.e., is also in its lower component. This alignment is a result of dynamics. 
Finally, we can redefine the phase of the doublet superfield Hd so that va is real and positive. This is not 
compulsory, but is the customary practice that allows us to define tan/3 to be real and positive. 
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eluding effects of third generation Yukawa couplings, but ignoring all flavour-mixing effects. 
These corrections, which effectively shift the Higgs boson SSB mass squared parameters 
by T, u and respectively, are evaluated by replacing ft t b, T m the standard relations by 
the (3,3) element of the corresponding Yukawa matrices in the basis where they are diago- 
nal at m t , and with the dimensionful parameters also replaced by the (3,3) element of the 
corresponding matrix (or the appropriate frozen value). 



6.6 Iterative Stage 

Now that the boundary conditions are defined at each of the three relevant scales, 
the iteration can begin. The subroutines DOWNMHIGHMZ and UPMZMHIGH2 implement the 
running in each direction. The iteration takes place a fixed number of times, chosen so that 
the RGEs reach a stable solution. 

Unless the user has answered '1' to using fixed thresholds from ISA JET in Question 
4. of INRGE.DAT, the programme will alter all the thresholds on each downwards run, except 
those for the gluinos and heavy Higgs fields, which are fixed at the locations obtained by 
ISAJET. 

The running is carried out as follows: 

Downwards running 

The subroutine first runs from the GUT scale to the highest threshold with the 
number of steps given by the variable NSTEP. During the ith iteration, this variable is 
has the value 100 x (1.6)*, until we reach iteration number 5 at which point it becomes 
fixed at 100 x (1.6) 5 . If this is the first iteration, the thresholds are taken to be the ISA- 
JET thresholds. Running continues between thresholds (inserting the boundary conditions 
at mu and -Mstjsy when necessary) with the number of steps given by 

Number of Steps = - l 1 "^*/^)! x ( 2 5 x NSTEP) , (6.7) 
log (M HIGH /m t ) v ' v ' 

where Q\ and Q2 are the scales of the two thresholds between which the running is being 
carried out. The factor 25 was chosen to ensure enough sampling between the thresholds 
without unnecessarily slowing down processing time. 

At each step, the SSB sfermion mass matrices are diagonalised. If the user has not 
fixed the thresholds to be the same as those obtained from ISAJET, the derived eigenvalues 
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are checked to see if any matter sfermions have decoupled. When one of these sparticles 
does decouple, the subroutine CHDEC carries out the following procedure: 

1. Set the corresponding 6>? to zero. 

Ji 

2. Store the eigenvectors of the mass matrix so that the rotation between the diagonal 
basis at the decoupling scale and the original current basis is saved. 

3. Store all the entries of the mass matrix itself so that they can be used as a boundary 
condition when running up. 

4. Call the subroutine REMSF so that, in the basis where the mass matrix is diagonal, 
the entry corresponding to the decoupled particle is set to zero. This ensures that 
the eigenvector for this particle is removed from the original current basis matrix and 
cannot influence further downward running. 

The subroutine continues to run down until it reaches nit, where it rotates from the current 
basis back to the basis in which the Yukawa matrices are diagonal. ^ Running resumes using 
the SM DR RGEs in SMRGEDR, without decoupling the top quark, to the scale Mz- 

The integration of the RGEs is carried out in the same way as ISA JET, i.e., by the 
CERNLIB routine RKSTP. The RGE subroutine, RGE646, contains RGEs for all couplings with 
and without tildes and with full thresholds for the one-loop running. The quartic couplings 
are entered separately, but they are set to be equal to their SM counterparts since the 
RGEs for the quartics are unavailable at this time. In addition, RGE646 contains the two- 
loop terms from the RGEs, which depend only on the MSSM values of the couplings. In 
order to obtain an estimate of the two-loop contributions, the pure MSSM RGEs are solved 
even below all thresholds and these MSSM couplings are used for the two-loop level running 
of the SUSY couplings. This is acceptable since we are only trying to achieve two-loop level 
accuracy, and threshold effects in the two-loop terms are numerically much smaller. 

Once we have decoupled at least one of the matter sfermions, the right-hand side 
of the RGEs are calculated in the basis where the squark/slepton matrices are diagonal at 
the decoupling scale. This is still a current basis, since the quarks and leptons are rotated 
by the same amount as the squarks and sleptons. RGE646 rotates all couplings into this 

^Rather than diagonalise the Yukawa matrices at this scale, RGEFLAV uses the rotation matrices, ~V l,r(u, d) 
that were used on the first upwards running. Practically, this means that the rotation matrices are one of 
the boundary conditions on our iterative running. 
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basis when calculating the right-hand side of the RGEs and, at the end of the subroutine, 
rotates the result back to the original current basis. 

Note that if the location of the thresholds is altered every iteration, the Yukawa 
couplings are unable to reach a convergent solution. This is because moving the thresholds 
can disrupt the fine cancellation that is required to obtain vanishing values at m t for the 
off-diagonal elements of the Yukawa matrices in their mass basis. We therefore only allow 
RGEFLAV to change the locations of the thresholds for the first five iterations. This ensures 
that the Yukawa couplings converge as closely as allowed by the numerical accuracy of the 
machine. 

Upwards running 

Before commencing the run back up to the GUT scale, the boundary conditions 
at Mz — the gauge couplings and Yukawa coupling matrices in the quark mass basis - 
are reset. Running then continues to nit, where the top quark Yukawa coupling is reset and 
the Yukawa matrices are rotated into the current basis, and then continues again until the 
first threshold above m< is reached. 

The upwards running makes no changes to the thresholds. At each step, the 
subroutine checks whether a threshold has been passed. If so, the 6 for the threshold in 
question is set to 1 , and if this is a matter sfermion threshold the soft mass matrix is set to 
the value which was saved at this point during the downward run. 

The RGE subroutine RGE646 is used just as with downwards running. When we 
have some but not all matter sfermions present in the theory, we rotate to the basis where 
the soft mass matrix for the sfermions in question is diagonal at the scale of decoupling. 
We know what this rotation is since it too was saved in the previous run down. 

Once all the thresholds have been passed, the RGEs are equivalent to the standard 
MSSM RGEs and running continues in a straightforward manner until the high scale is 
reached. Residual inaccuracies in the running mean that the tilde-couplings are not precisely 
equal to their non-tilde counterparts once we have passed the highest SUSY threshold. Since 
these differences can feed back into the other couplings via the RGEs, we set the tilde- 
couplings (except fi) equal to the usual SUSY couplings once the highest SUSY threshold 
has been passed. Also, if the answer to Question 6. is '0', since the scale at which the gauge 
couplings unify may be altered by the location of the thresholds, we allow the running to 
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continue past the unification scale from the previous iteration, and increase the number of 
steps for this iteration so that the step-size remains constant. 

We have checked that if we do not reset the weak scale boundary conditions, and 
instead use the final values from the previous call to DOWNMHIGHMZ, the upwards running is 
precisely the same as the downwards running. 

6.7 Technical Aside: Numerical Accuracy of the Rotation 

Before concluding this chapter with a discussion of the output, we digress momen- 
tarily to draw attention to a technical issue which is important if the squark thresholds are 
not entirely degenerate. As discussed earlier, for values of Q below the highest sfermion 
threshold we need to rotate at each step to the basis where the SSB squark mass squared 
matrices are diagonal, which, in turn, requires us to obtain the unitary matrix R, that 
relates our "standard" basis to this "squark mass basis". After evaluating the right-hand 
side of the RGEs, we always rotate back to our "standard" basis using R.I. While this 
is straightforward in principle, the practical problem is that when two squark eigenvalues 
become very close — this is always the case in mSUGRA — R.Rl has non-zero off-diagonal 
elements due to numerical noise at the 10~ 10 level. This ruins the delicate cancellations 
necessary to obtain the tiny magnitudes of some of the off-diagonal elements of the squark 
mass matrices that we will see in the next chapter. The subroutine ORTH is called in several 
places to fix this problem. A further discussion of this effect and the procedure we use to 
mitigate against the numerical instability can be found in Appendix F. 

6.8 RGEFLAV Output 

The iterative section exits at the high scale after resetting the GUT scale boundary 
conditions. In order to provide useful output, the code makes one final downward run to 
run- This scale was chosen due to its significance as a point where a number of the operators 
in the Lagrangian change, however, the output scale could have been chosen to be anywhere 
between the two extremes of the running. It is expected that mu will be fairly close to the 
scale at which the user will be using the couplings for their calculations. 

Before passing on the output to an external routine, all couplings are rotated to 
the current basis chosen by the user in Question 17. of INRGE.DAT, where either the up- 
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or the down- type Yukawa couplings are diagonal at m t . For any specific calculation the 
user can then simply evolve these couplings to higher or lower scales as desired, without the 
need to iterate. We will use the output, in the basis where the up-type Yukawa coupling is 
diagonal, in the subroutine STDEC described in Sec. 8.4, to calculate the t\ decay rate. 
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Chapter 7 

Sample Solutions to the RGEs 



Having constructed the RGEs for all dimensionless and dimensionful parameters 
of the MSSM (except the scalar quartic couplings), and developed a program to solve these 
RGEs subject to SM and GUT scale boundary conditions, we are now in a position to 
consider their solutions. 

We will focus on mSUGRA as a simple model for GUT scale inputs, and discuss the 
consequences of splitting between SUSY couplings and their SM counterparts. In addition, 
since our eventual goal is to study flavour-changing decays of squarks, we will consider the 
flavour structure of the couplings, and how certain aspects of our general GUT scale inputs 
in (6.4) can affect the off-diagonal elements of the various couplings renormalised at the 
scale of the sparticle mass, i.e., at scales Q typically between ~ 100 GeV and a few TeV. 

To simplify the discussion further, we begin by reintroducing our basic model of 
decoupling from Sec. 4.4.2, in which there are only two thresholds. We assume in this simple 
model that the sfermions and the electroweak gauginos are all at a mass scale ~ 600 GeV 
while the heavy Higgs bosons and gluinos have a mass ~ 2 TeV. In other words, our theory is 
supersymmetric for Q > 2 TeV, contains only sfermions, charginos and neutralinos together 
with SM particles for 600 GeV < Q < 2 TeV, and is the SM for Q < 600 GeV. We use this 
model to examine the gauge and gaugino couplings before turning to quark Yukawa and 
Yukawa-like higgsino couplings where we will eventually consider more realistic scenarios. 

7.1 Gauge and Gaugino Couplings 

The RGEs of the gauge couplings, as in (5.10)-(5.12), are flavour diagonal under 
renormalisation group evolution, even at higher loops. Although the Yukawa couplings enter 
at two-loop level, they only appear as 'flavour-blind' traces. As a result, notwithstanding 
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Figure 7.1: The evolution of the magnitudes of the off-diagonal elements of for our simplified 
threshold scenario where the thresholds (denoted by arrows) are clustered at 2 TeV and 600 GeV. 
Also shown for comparison is the running of the (2, 3) element of the up Yukawa coupling. The 
legend is in the same order (from top to bottom) as the curves. 

that the gauge couplings are affected by the location of the thresholds via the factors N-p 
and 8-p, these effects cannot contribute to additional flavour- violation. 

On the other hand, we saw in Sec. 4.4.2 that the form of the RGEs for the (Yukawa- 
like) squark-quark-gaugino couplings means that they develop a non-trivial matrix structure 
once the heaviest SUSY particle has been decoupled. Although these couplings — and 
similarly the higgsino couplings — cease to have any meaning below the lower threshold, 
there is still the opportunity in the intermediate regime for these couplings to develop 
significantly different values from their SM counterparts. Contributions to the g* running 
from terms involving Higgs boson and higgsino couplings, which cancel above all thresholds, 
result in non-zero off-diagonal elements when Q < ran. 

As a supplement to Fig. 4.3, where we showed the off-diagonal elements of the 
U(l) gaugino couplings g"^ and g' UR , we display in Fig. 7.1 the off-diagonal elements of g^ 
(of course, the coupling g UR does not exist). To construct these curves, we have solved the 
RGE in (A. 17) along with the rest of the gauge and Yukawa couplings, which form a closed 
set. As in Sec. 4.4.2, the boundary condition on g*^ is 



gJ(Q = 2 TeV) = </<%, 
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and the curves terminate at 600 GeV when the squarks and gauginos decouple from the 
theory. Once again, we show the (2, 3) off-diagonal element of f u for comparison. 
We note the following features: 

• The variation of the vertical scale is over a much smaller range than Fig. 4.3. We 
see that the approximate scale of the elements of g^ is the same as for g /( ^ and 
therefore, as with the U(l) coupling, the elements are of similar size to the largest 
of the off-diagonal Yukawa couplings. In models where the Yukawa couplings are the 
only source of flavour-violation, we would conclude that a study of flavour should 
necessarily include the matrix structure of the gaugino couplings. 

• The overall size of the left-squark gaugino couplings are several orders of magnitude 
larger than those of g /n «. We see from (A. 10), (A. 12) and (A. 17), that the left-handed 
gaugino coupling RGEs depend on both the up-type and down- type Yukawa couplings, 
whereas the right-handed RGE depends only on the up- type Yukawa couplings. Since 
the scale is set by the Yukawa matrices, which (unlike the gaugino coupling matrices) 
have non-zero flavour structure at Q = m#, we find that in our "standard" cur- 
rent basis (which was defined in Sec. 3.1.2) the down-type Yukawa terms have larger 
contributions to the running. 

• Finally, as pointed out in Sec. 4.4.2, the magnitudes of the off-diagonal elements of 
the gaugino coupling matrices are (almost) symmetric under interchange of the two 
indices. If we ignore the difference between the Higgs boson and higgsino coupling 
matrices, to be discussed next, we find that the right-hand side of the RGEs preserves 
the approximate hermiticity of the gaugino coupling matrices, which are proportional 
to the unit matrix at Q = run- 

Before moving on to the Higgs boson and higgsino couplings, we briefly mention 
the diagonal entries of the gaugino couplings. Of course, these are equal to the usual gauge 
couplings above Q = inn, but, as with the off-diagonal elements, can evolve differently below 
the scale of the highest threshold. Indeed, we find that their evolution below Q = mu in 
our simplified scenario is in the opposite direction to the gauge couplings [39] because, for 
Q < nig, the evolution of the diagonal gaugino couplings now depends on the much larger 
gluon coupling, even at the one loop level. When the gluino decouples from the theory, we 
can see from, for example, Eq. (A. 10), that a cancellation which was taking place between 
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Figure 7.2: Evolution of the magnitudes of the elements of the up-quark Yukawa coupling matrix, f„, 
for the MSSM with thresholds, shown by arrows, clustered at 600 GeV and 2 TeV as discussed in the 
text. Above m# (= 2 TeV) we plot |(f u )jj| whereas below m#, where the effective theory includes 
just one scalar Higgs doublet, we plot |(A u )y| / sin/3 which is equal to |(fu)y| at Q = inn = 2 TeV. 



gs and g s terms no longer occurs, causing a large change in the /3- function as 8g — ► 0. We 
note that, as can be seen in Fig. 4. of Ref. [39], although the splitting between the thresholds 
in this scenario is not large, the gaugino and gauge couplings develop a difference of ~ 4%. 
The existence of a difference in these couplings has been discussed in Refs. [49, 50] where the 
possibility for this splitting to give an insight into the sparticle mass spectrum is explored. 



7.2 Quark Yukawa and Higgsino Couplings 



In this section we consider both the usual quark Yukawa couplings f u ,d,e an d the 
higgsino Yukawa-like couplings, fjfj e - As an introduction to the more realistic scenarios to 
follow, we begin by exploring further our simplified threshold model. 

The running of the elements of f u is shown in Fig. 7.2. In order that all entries 



of the Yukawa matrix be displayed for the whole range of the running, we plot 



(fu) 



for 



Q > run and 



(A u )j ■ / sin/3 for Q < rriH- These two lines join continuously at Q = ran 



and, since we are plotting this figure in the basis where the up quark Yukawa couplings are 
diagonal at mt, the off-diagonal entries are seen to go to zero at this scale. 
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The most striking feature of Fig. 7.2, which we have seen before in Figs. 4.3 and 7.1 
but not specifically addressed, is the presence, in the off-diagonal entries, of a sharp dip in 
the range Q ~ 850 GeV, which similarly appears in the higgsino couplings. This feature is 
particularly important for our purposes since (in a more realistic scenario) it may appear 
close to the sale at which we decouple the squarks and directly enter the amplitudes for 
flavour-violating squark decays that will be calculated at this scale. 

A detailed discussion of this feature can be found in Ref. [39]. We simply note 
here that the off-diagonal up-type Yukawa couplings in our "standard" current basis of 
Sec. 3.1.2 obtain their overall scale from the 'seed-terms' in the RGEs that depend on 
id or ff\ In contrast to the up- type Yukawa terms, for which all off-diagonal entries are 
small, these terms contain the KM matrix, whose presence allows for non-zero off-diagonal 
contributions. Given that the contribution to the RGE from these id terms is dominant, 
we find that the /^-function undergoes a sign change in morphing from pure SM to pure 
MSSM running. This sign change ensures that the coupling in question alters direction and 
heads back towards zero, which was its starting value. Due to the fact that — as is the 
case in Fig. 7.2 — the real and imaginary parts of each specific off-diagonal coupling both 
head towards zero simultaneously, the magnitude of the Yukawa matrix entry will vanish 
at some specified point, producing a pronounced dip since we are plotting the logarithm 
of this magnitude. We have shown that in a more simplified threshold scenario, where all 
SUSY thresholds are located at a single point, the dips for all the off-diagonal entries occur 
at a single scale [39] , but the blurred out location of the zeros is a result of our using a more 
realistic spectrum. 

Having considered the running of the usual Yukawa couplings, we now turn to 
consider the off-diagonal elements of the f* matrices. We know that these will deviate 
from the elements of the corresponding Yukawa coupling matrix only below the thresholds 
at Q = ran, but have magnitudes similar to the corresponding Yukawa coupling matrix 
elements from Fig. 7.2. It is the difference between the couplings of Higgs bosons and 
higgsinos that will be the main focus of our attention. In Fig. 7.3, we show the evolution of 
the real and imaginary parts of the (1,3) element of (i) the Yukawa coupling matrix, i u for 
Q > mu and \ u j sin/3 (which connects continuously to i u ) for Q < mn, (H-) the higgsino 
coupling matrices i® , and f" fl whose evolution is given in (A.4) and (A. 7), respectively, 
of Appendix A. There is no particular reason for our choice of the (1,3) element (which 
happens to have comparable real and imaginary pieces) for the illustration in the figure. 
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Figure 7.3: The real (R) and imaginary (I) part of the (1, 3) element of the up-quark Yukawa coupling 
matrix, f„, along with the corresponding elements of the matrices and for the MSSM with 
the spectrum as in Fig. 4.3. For the solid black and red lines, we plot the elements of the matrix 
X u / sin j3 below Q = m#, as in Fig. 7.2. The main figure zooms in on the low end of the range of 
Q where the Higgs boson and higgsino couplings differ from one another, while the inset shows the 
evolution all the way to Mqut- 

Here we have focussed on the lower end of Q where the and f u couplings are different, 
while the inset shows the evolution all the way to Mqut- Several points may be worthy of 
notice. 

• For Q > ran where the effective theory is supersymmetric, we see that the real and 
imaginary parts of the Higgs boson and higgsino couplings separately come together 
as expected. 

• For Q < rriH, the higgsino couplings are split from the corresponding Higgs boson 
couplings as well as from one another by a factor of several. For instance at the scale 
of squark masses, the real parts of the (1,3) element of both f„ and f" fl are quite 
different from the real parts of (f n )i3, and likewise for the imaginary parts. It seems 
to us that the use of the evolved Higgs boson coupling in place of the corresponding 
higgsino coupling could be a poor approximation. 
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• Notice that while the real and imaginary parts of (fu)\3 and (f„ H )i3 come to zero at 
the same point, the position of the zero differs for the two couplings. If the magnitude 
of the higgsino couplings were plotted on a logarithmic scale, as for the usual Yukawa 
couplings in Fig. 7.2, we would see dips in their values at ~ 1400 GeV and ~ 1000 GeV 
respectively. 

The variation of the diagonal entries of the usual Yukawa matrices and their hig- 
gsino counterparts can be found in Fig. 5. of Ref. [39]. We find a large difference between 
the calculation of the complete RGEs in our simple threshold scenario compared to the 
two-loop calculation without thresholds. In addition, we see that the higgsino couplings 
can evolve to significantly different values from the usual Yukawa couplings over the range 
600 GeV < Q < 2 TeV. As with the diagonal entries of the gauge and gaugino couplings, 
this difference is due to the appearance in the RGE of a dependence on the strong couplings 
when certain cancellations between SUSY and SM terms are no longer complete. 

7.2.1 mSUGRA 

Having considered the qualitative aspects of introducing thresholds into the dimen- 
sionless RGEs for the MSSM, we move on to continue some more specific models, where 
GUT scale parameters dictate the overall spectrum. Since the RGEs for the gauge and 
Yukawa couplings depend on the dimensionful parameters through the thresholds, we must 
use more realistic inputs constructed from the relations in Sec. 6.4. 

To this end we choose an illustrative mSUGRA point with tuq = 200 GeV, m 1 / 2 = 
—400 GeV, Aq = —200 GeV, tan (3 = 10 and ^ > for which the particle spectrum is as in 
Table. 7.1 and begin by again showing the magnitude of the up Yukawa coupling matrix^, in 
Fig. 7.4. In a similar manner to Yukawa matrix in the simplified spectrum (Fig. 7.2), there 
are (approximately) aligned dips at Q ~ 650 GeV, common to all the off-diagonal elements, 
that occur because of the change in the sign of the coefficient of the f^fj-type terms that 
drive the growth of these off-diagonal elements from zero at Q = rat. The presence of the 
second zero at the higher value of Q, just in the (2,3) and (1,3) elements, is accidental. It 
occurs because of conspiracies between terms in the corresponding /^-function as the left- 
type squarks are decoupled. Notice that the lowest four curves, though they do not have 

^We have broken the scale at Q ~ 2 TeV to more clearly show features of the running at small Q. Note 
that this gives the impression of a change in slope of all the curves between the high Q and low Q regime, 
but this is just an artifact of the dramatic change in scale. 
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Msusy 703 GeV 

Higginos (/i) 538 GeV 

Gluinos (m~ g ) 941 GeV 

H, H± (m H ) 631 GeV 

Bino(|Mi|) 166 GeV 

Winos (\M 2 \) 315 GeV 

(u L , d L ), (c L , s L ), (h, b L ) 837 GeV, 837 GeV, 763 GeV 

u R , c R , t R 809 GeV, 809 GeV, 645 GeV 

d R , s R , b R 806 GeV, 806 GeV, 801 GeV 

(j>eL,e L ), {vtL,t l ) 331 GeV, 331 GeV, 329 GeV 

Sr, ji R , f R 249 GeV, 249 GeV, 245 GeV 

Table 7.1: The approximate location of the thresholds for the mSUGRA case in Fig. 7.4 and most 
subsequent mSUGRA examples. 



mSUGRA: m ( =200 GeV, m 1/2 =-400 GeV, A () =-200 GeV, tanp=10, \i>0 




Q (GeV) 

Figure 7.4: Evolution of the magnitudes of the complex elements of the up-quark Yukawa coupling 
matrix in the mSUGRA model with m = 200 GeV, m 1/2 = -400 GeV, A = -200 GeV, tan (3 = 10 
and /i > in the basis where this matrix is diagonal at Q = mt- For Q > ran we plot |(f u )jj | whereas 
for Q < rriH, we plot |(A u )ij| / sin/? which is equal to |(f u )y| at Q = mn- In all the figures we take 
m t = 172 GeV. 
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mSUGRA: m =200 GeV, m 1/2 =-400 GeV, A ( =-200 GeV, tanp=10, \i>0 
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Figure 7.5: Evolution of the magnitudes of the complex elements of the down-quark Yukawa coupling 
matrix in the same mSUGRA model as Fig. 7.4. In the left-frame, we show the elements of the matrix 
in the "standard" basis, whereas in the right frame, we show these elements in the basis where the 
Xd is diagonal at Q = m t . As with Fig. 7.4, for Q > ran we plot |(fd)ij| whereas for Q < ran we 
plot \(Xd)ij\ / cos/3. 

this additional dip, show kinks at these same values of Q, corresponding to the decoupling 
of these squarks. For several other mSUGRA cases, we have checked that while squark 
decoupling causes kinks in the curves, the coupling does not drop to even close to zero for 
a second time, in contrast to the behaviour in our illustrative example in the figure. 

The evolution of the down-quark Yukawa coupling matrix for the same mSUGRA 
point is shown in Fig. 7.5. We have shown these couplings both in our "standard" current 
basis from Sec. 3.1.2 (where the up-type Yukawa couplings are diagonal at Q = mt) in the 
left frame, and in the basis where the down-type Yukawa coupling matrix is diagonal at 
Q = uif in the right frame. The matrices in the two bases are connected by the KM matrix, 
as given by (3.11). 

The curves in the left frame are all smooth (except for the small kink in the curves 
for | (fd)is I that occurs because of the SUSY correction to the bottom quark Yukawa coupling 
[46]), and do not show the dip to zero that appeared in the previous figure. This is not 
surprising because underlying the explanation of this dip was the fact that the off-diagonal 
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elements evolved from zero at Q = mt [39]. Notice also that in this frame the off-diagonal 
elements are not necessarily smaller than the diagonal elements even for Q below the region 



The magnitudes of the off-diagonal matrix elements in the frame on the right, 
which do start at zero at Q = mt, show the anticipated aligned dips, except that the 
location of the dip is shifted considerably to the right relative to Fig. 7.4. This shift is 
not difficult to understand. The large top quark Yukawa coupling in the SM governs the 
evolution of the off-diagonal elements of the down- type Yukawa couplings, causing them 
to evolve much more rapidly from zero in Fig. 7.5. Therefore, in order to evolve back to 
zero after the sign flip in the /3-function due to the additional Higgs and SUSY particles, 
a longer evolution distance is needed. Furthermore, beyond the Higgs boson threshold the 
off-diagonal elements of f M in Fig. 7.4 are accelerated to zero on account of the fact that the 
down-type Yukawa couplings f<i, that enter in the evolution of these elements, are enhanced 
by a factor ~ 1/ cos/3, pushing the dip in this figure to a low value of Q. 

7.2.2 Additional GUT Scale Flavour- violating Terms 

We now turn to briefly discuss what happens when we allow non-universality of 
GUT scale squark mass parameters, so as to split the squarks more than in mSUGRA. Recall 
that if the squarks all decouple together, all that happens is a change in the slope of the 
/3- function. If, however, the squark masses are not all the same, our decoupling procedure 
entails an additional rotation to the squark mass basis. If this basis differs significantly 
from our "standard" basis in which the quark Yukawa coupling matrix is diagonal (at 
Q = mt), one may expect considerable deviation in the evolution of the off-diagonal elements 
from the mSUGRA case for Q values in-between the highest and lowest squark thresholds. 
The introduction of non-universal squark mass parameters via non- vanishing values of the 
constants Ru,d,e and Su,d,e in (6.4b) never leads to significant effects because the rotation 
from the "standard" basis to the squark mass basis is small by construction. 

The question then is whether we can have large deviations from Figs. 7.4 and 7.5 
via the Tu,d,e ( or Zu,d,e) matrices. To examine this, we set all GUT scale inputs to be the 
same as the mSUGRA case in Fig. 7.4 except 



of 1 TeV. 



m {U,D}0 - 

T UD = diag {10000, 40000, 90000} GeV 2 . 



(7.1a) 
(7.1b) 
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m () =200 GeV, m 1/2 =-400 GeV, A Q =-200 GeV, tan(3=10, u>0 

m 2 |UDK =0, T UD = diag( 10000,40000,90000) GeV 2 and V R (u,d)=l 




(2.3) " 
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Figure 7.6: The same as Fig. 7.4 except that the GUT scale right-handed squark mass matrices are 
now given by (7.1) in the "standard" current basis. 

We need, of course, to specify the basis in which the squark mass matrix is diagonal. If 
this is the "standard" current basis, with \~l(u) = Vji(u,d) = 11 at the weak scale, we 
see in Fig. 7.6 that although the right squark masses are now significantly split relative to 
the mSUGRA case, the Yukawa matrices are mostly unaltered. That is to say, except for 
detailed changes in the evolution for Q values in-between the squark thresholds (e.g., the 
additional structure in the dip at the higher value of Q in the (2,3) and (1,3) elements), 
the evolution is not different in any important way at large values of Q. By choosing to 
enter Tjj with diagonal (but unequal) entries, we are introducing only a mild amount of 
new flavour structure because the rotation between the "standard" basis and the squark 
mass basis is again small. 

If we take instead that the right-handed squark matrices in (7.1) are diagonal 
in a completely different basis specified by unitary matrices V^(-u), V/j(u) and V^j(d), 
entered in INRGE.DAT using (6.1) with randomly chosen values of a, 0, 7 and 5, we may 
expect the evolution of Yukawa couplings to depart from the corresponding evolution in 
mSUGRA. We emphasize that in this model, in which mfj D are diagonal in the basis where 
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m =200 GeV, m 1/2 =-400 GeV, A =-200 GeV, tan|3=10, n>0 

m 2 =0, T UD = diag( 10000,40000,90000) GeV 2 and V R (u,d)*l 



T» — | 1 1 1 





Q (GeV) 

Figure 7.7: The same as Fig. 7.4 except that the GUT scale right-handed squark mass matrices 
are now given by (7.1) in the basis specified by (7.2) of the text. Note, however, that these matrix 
elements are shown in the "standard" current basis where the matrix is diagonal at Q — m t . 

the corresponding Yukawa coupling matrices are very off-diagonal, includes potentially large 
flavour-violation in the singlet squark SSB mass matrices. 

In Fig. 7.7 we illustrate the evolution of the magnitudes of the elements of the 
up-quark Yukawa coupling matrix with randomly chosen matrices, 

V L (u) : a = 2.053, j3 = 0.254, 7 = 2.030, 5p = 0.4829 (7.2a) 
Vr(u) : a = 1.188, j3 = 2.218, 7 = 0.763, 5p = 0.87 (7.2b) 
V R {d) : a = 1.904, j3 = 2.947, 7 = 1.847, 5p = 1.14 . (7.2c) 

Note that, just as in Fig. 7.4, we have plotted these elements in our "standard" basis where 
the Yukawa coupling matrix is diagonal at Q = mt, and that in this basis the GUT scale 
matrices mfj and will be random Hermitian matrices. We see that the evolution of 
the diagonal elements is not significantly altered from mSUGRA. This is because, although 
there is a large mismatch between the squark mass basis and our "standard" basis, this 
mismatch is operative over the small range of Q between the highest and lowest thresholds 
and so has little impact on the largest elements. The two largest off-diagonal entries {i.e., 
(f«)23 and (f u )i3) similarly do not change significantly from Fig. 7.4, but all other entries 
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are greatly altered. While it may seem that the values of these matrix elements at a large 
scale is quite irrelevant phenomenologically, the altered form of the Yukawa coupling matrix 
at the high scale could be of relevance to model-builders. 

We have seen that the magnitudes of the off-diagonal elements in Fig. 7.7 remain 
small because the splitting in the squark spectrum is limited to 0(100 — 1000) GeV. The 
natural question then is whether we can get these to be larger by choosing extreme splitting 
between the squark eigenvalues. Even aside from potential flavour-changing effects, this 
is not easy. In general, such a GUT scale splitting also has a large value for S = fn\ u — 
fn 2 Hd + Tr m.Q — m 2 L — 2m^ + mjj + m 2 E , which pulls the other squarks also to large 
masses, so the squark mass splitting is reduced by RGE effects. It may be possible to 
obtain split squarks by very carefully adjusting S to be zero, since S is then invariant under 
renormalisation group evolution, but we have not investigated this here. 



7.3 Gaugino Mass Parameters 

We have seen in Sec. 5.4 that the evolution of the electroweak gaugino mass pa- 
rameters acquires a dependence on the ^-parameter (and vice-versa) when threshold effects 
from splitting in the Higgs sector is taken into account. In models where m# » \fi\,M\ t 2, 
the effect of the term explicitly dependent on fi in (B.2)-(B.5) (and the corresponding terms 
dependent on the gaugino mass parameters in the RGE for n) may be significant, so that the 
relation from Eq. (4.6), M2/M1 = a^/cti, which is expected in many models, is modified. 
We should keep in mind that two loop terms will, in general, also alter this relation. Our 
point is that we should expect threshold corrections from the /x term in the RGE, as well as 
from the decoupling of sfermions, to be comparable to the two loop modifications, and so 
need to be included in a quantitative analysis. Within the mSUGRA context we have small 
values of and hence ran ^ | | , in the hyperbolic branch/focus point (HB/FP) region 
[51] which occurs for large values of mo, and is one of the regions favoured by the association 
of a light neutralino with the dark matter relic-density in connection with the comments in 
Sec. 2.4.1. We mention here that the location of this HB/FP region is significantly altered 
by the inclusion of the threshold corrections. 
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mSUGRA: m =3075 GeV, m 1/2 =-600 GeV, a =0 GeV, tan(3=10, u>0 




100 1000 10000 

Q (GeV) 



Figure 7.8: The evolution of the gaugino mass ratio, Af 2 /Af l5 for the mSUGRA point: m = 
3075 GeV, mj/2 = —600 GeV, Aq = 0, tan/3 — 10 and fi > 0. We plot the ratio in three cases: 
one-loop with thresholds; two-loop with thresholds; and two-loop without thresholds. The gauge 
coupling ratio, otz/oti, is also shown for comparison. 

7.3.1 mSUGRA 

The above considerations led us to examine the evolution of the gauge couplings 
and the electroweak gaugino parameters for a relic-density-consistent mSUGRA model point 
in the HB/FP region with m = 3075 GeV, m 1/2 = -600 GeV, A = 0, tan (3 = 10 and 
/i > 0, for which (/i,Mi,M 2 ) ~ (304,-258,-495) GeV at the weak scale." 1 " As can be 
seen in Fig. 7.8, at the two loop level with all threshold effects included, M2/M1 = 1.919 
to be compared to M2/M1 = 1.881 obtained without including threshold effects. Thus, 
although threshold effects actually bring us closer to 02/01 (= 1.964), their inclusion is 
clearly necessary for a quantitative analysis of mass parameters that may be extracted at 
an e + e~ linear collider, where a precision of better than 1% will be possible if charginos are 
kinematically accessible. 

^Without threshold corrections, ISAJET gives a similar spectrum for mo — 3660 GeV. 
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Q (GeV) 

Figure 7.9: The evolution of the gaugino mass ratio, M2/M1, at the one (dashed) and two (solid) 
loop levels, along with the two-loop evolution of c^/ai, for a split SUSY model where scalars masses 
are around 10 7 GeV, with gauginos and higgsinos at the weak scale. Gaugino mass unification is 
assumed and the gaugino mass parameters are negative. The other parameters are as mentioned in 
the text. 

7.3.2 Split Supersymmetry Models 

The threshold corrections to gaugino mass parameters can be much larger in the 
so-called split SUSY model [52] that has received considerable attention in the recent liter- 
ature. In these models, the naturalness of the scalar Higgs sector (which we view as one of 
the primary motivations for weak scale SUSY) is abandoned, while gauge coupling unifica- 
tion and the neutralino dark matter candidate of -R-parity violating models are preserved. 
Gaugino mass parameters and \fi\ are assumed to be at the weak scale, while scalar masses 
are at an intermediate scale. This means that sfermion masses as well as ran are very large 
(with the SM Higgs doublet fine-tuned to be light), so that charginos and neutralinos are 
the only new particles (other than a SM Higgs boson) at the weak scale. 

As an illustration, in Fig. 7.9 we plot the variation of the ratio M2/M1 at the one- 
loop level (dashed), as well as at the two- loop level (solid), with the renormalisation scale 
Q, along with the two-loop value of a^/ai. We show results, first where the value of is 
set exactly in-between \M\\ and \M2\ so that the lightest neutralino acquires a significant 
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higgsino component, to qualitatively mimic mixed higgsino dark matter. t We have checked 
that the M2/M1 values do not change in a significant way for yet larger values of tan/3. It 
is clear that the relation M2/M1 = 02/01 is violated at the several percent level by the 
threshold corrections, without which the M2/M1 lines would have continued to low values 
with the same slope that they have above 10 7 GeV. For this point, the 2sc x /x term that 
explicitly appears in the RGE is very small so that the result is independent of the sign of 
/j,: most of the difference is an effect of the sfermion loop contributions being switched off 
below 10 7 GeV. 

To gain some idea of how large the effect of this /j, term might be, we have also 
shown M2/M1 for fi = — 2 TeV, with the multiple (/XM2) > and tan/3 = 4. Since such a 
large value of fi would be totally incompatible with the measured relic density, and small 
values of tan/? are unnatural in these models without some modification to the EWSB 
sector, the reader should view these curves only as a guide to how much the gaugino mass 
ratio may deviate from its "unification value". The difference between the two cases is 
almost entirely due to the different choice of fi. We see that the gaugino mass unification 
condition will, in this case, be violated by ~ 10%. If instead we choose the opposite sign 
for (11M2), but keep |/i| = 2 TeV, this "large \n\ line" would be lower than the line with 
fi = (Mi + M2)/2 by about the same amount that it is higher than this line in the figure. 
Clearly, increasing the splitting between the scalar and the gaugino/higgsino sector of the 
theory will cause even further violation of the unification condition, and ~ 20% effects 
appear to be plausible if the scalars are instead at the 10 11 GeV scale. 

7.4 SSB a-parameters: mSUGRA and Non-universal Inputs 

Returning to the mSUGRA case considered in Figs. 7.4 and 7.5, we show the 
magnitudes of the individual entries of the trilinear coupling matrices a u and in Fig. 7.10. 
There are three separate plots: (a) | (a u )^ | and (b) | (a^)^- | in our "standard" current basis 
where the up- type Yukawa coupling matrix is diagonal at Q = m t , along with | (a^)^- | in 

^In the absence of a real theory of split SUSY, we should view this figure only as a qualitative illustration of 
potentially large threshold effects. Here, we take the sfermion mass parameters to be 10 7 GeV at Q = Mqut, 
mi/2 = —350 GeV and Ao = 0. Since it is not possible to satisfy the EWSB conditions except when tan/3 is 
hierarchically large - this would cause down-type Yukawa couplings to become non-perturbative - we treat 
fj, and tan/3 as phenomenological parameters, and fix m H to be 10 7 GeV in this figure. The parameters 
m 2 Hu and m 2 Hd (indeed all scalar mass parameters) are never needed since the RGEs for gaugino masses, /u, 
the a-parameters, and the dimensionless couplings form a closed set even at the two-loop level. Sfermion 
masses only enter via the location of thresholds. 
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Figure 7.10: The magnitude of the elements of the trilinear coupling matrix a Mj( f for the mSUGRA 
model in Fig. 7.4. We show (a) | (a K ) tJ |, (6) | (a<j)y \, m the "standard" current basis, and (c) | (a^ ■ | 
in the basis where down-quark Yukawa couplings are diagonal at Q = m t . The curves extend between 
Q = m H and Q = M GUT . 



the basis where the down- type Yukawa coupling matrix is diagonal at Q = mt in frame (c). 
We terminate the curves at Q = mu since below this scale we have a single Higgs scalar 
doublet model, and the trilinear couplings a Ut d evolve only as part of a linear combination 
with A*f^' d as discussed in Sec. 5.5. 

We see that the curves in frame (a) show a simple dip structure indicating that the 
real and imaginary parts of a. u ^ are really monotonic functions of Q that pass through zero 
together, in a manner similar to the elements of the Yukawa coupling matrix in Fig. 7.3. 
The actual location of the zero is somewhat harder to analyse than in the case of the Yukawa 
couplings because even though a u obtains off-diagonal components only because the down- 
quark Yukawa matrix is not diagonal at Q = mt, the matrix a u is off-diagonal even at the 
GUT scale. 

The off-diagonal elements of in frame (b) start off with a much larger magnitude 
in the "standard" basis at Q = Mgut because the corresponding Yukawa coupling matrix 
has large off-diagonal pieces. In this case, the evolution of these off-diagonal elements 
receives significant contributions from all entries in the RGE (unlike the evolution of the 
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Figure 7.11: The evolution of (a) Re(a„) 23 and (b) Re(a u ) 32 for a model with the GUT scale values 
of a-parameters set as in (6.4c) for several values of W and X shown in the same order as the legend 
(with A{ U ( j}o — A))- The GUT scale SSB scalar and gaugino mass parameters are assumed to be 
universal, with the same value as in Fig. 7.10. We terminate the curves at Q = mu- 

off-diagonal elements in frame (a) or of the off-diagonal Yukawa couplings discussed earlier, 
where contributions from the off-diagonal down- type Yukawa matrices govern the evolution), 
and never go through zero; the situation is similar to that in the first frame of Fig. 7.5. 

We see from frame (c) that the magnitudes of (a^)^- in the basis that the down- 
type Yukawa coupling matrix is diagonal at Q = mt again show the characteristic dip 
structure indicating that the off-diagonal elements increase in magnitude from their value 
at Q = Mgtjt to some maximum magnitude at an intermediate scale, but then smoothly 
reverse direction and thereafter evolve monotonically through zero to the low scale. The 
elements in frames (6) and (c) are, of course related by, (3.12). We note that the off-diagonal 
elements in frame (c), because these are driven by the larger "up- type" Yukawa couplings, 
are bigger than those in frame (a) where it is the down-type Yukawa coupling matrix that 
largely determines the entries. 

Finally, in Fig. 7.11 we consider a model with non-universal values of a-parameters, 
but where these are not a new source of flavour-violation. Specifically, we consider a model 
with the same values for mSUGRA parameters as in Fig. 7.10, but with non-zero values for 
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W U: d and X Uj( i in (6.4c) (with A^ ud y = Aq), and illustrate the evolution of (a) Re(a n ) 23 
and (b) Re(a u ) 32 . We have checked that the imaginary parts of these matrix elements are 
about four orders of magnitude smaller. 

The striking feature of frame (a) is that the various curves, which start with very 
different values of Re(a u ) 23 at Q = Mgut, appear to focus to a common value at the low 
scale. We have checked, however, that although they all cross at Q ~ 1.5 TeV, they do not 
all converge at precisely the same value of Q. This apparent convergence, which persists 
for other values of mSUGRA parameters, is sensitively dependent on the special GUT scale 
boundary conditions for a u that we have used. We have checked that if instead we use a 
general matrix Z u in (6.4c), the corresponding evolution is completely different. We do not 
have a good explanation for the seeming convergence in frame (a), and only note that it is 
not generic to all elements of & u as evidenced, for example, by the corresponding evolution 
of Re(a u ) 32 in frame (b) of the figure. 

7.5 Soft Masses 

Our discussion of the evolution of the scalar mass SSB parameters will have a 
similar progression to the previous sections, beginning with an examination of our canonical 
mSUGRA point. We will then proceed to introduce non-universal GUT scale boundary 
conditions, first where there is no additional source of flavour- violation other than the 
superpotential Yukawa couplings, and then when we set Tjj,d / in the same model as 
Fig. 7.7. 

7.5.1 mSUGRA 

We begin by showing in Fig. 7.12 the evolution of the magnitudes of mfj in our 
"standard" current basis for the mSUGRA model with the same parameters as in Fig. 7.4. 
The diagonal matrix elements start from a common value rriQ and increase as we go to the 
weak scale because of gauge (and gaugino) interactions. Although the splitting between 
the (1,1) and (2,2) elements that occurs because of the Yukawa couplings is too small to be 
visible in the figure, the (3,3) element is reduced significantly on account of the large (3,3) 
entry in f u . Notice that the curves become flat once the squarks are all decoupled. 

The magnitudes of the three independent off-diagonal elements of the Hermitian 
matrix start from zero at Q = Mgut, and rapidly rise because f u has off-diagonal 
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Figure 7.12: The scale dependence of the magnitudes of the entries of for the mSUGRA model 
with parameters as in Fig. 7.4 in our "standard" basis. Note the we have broken the vertical scale 
at 0.1 GeV 2 and also used a different scale above this to better show the splitting of the (3,3) entry 
from the other diagonal entries. 

entries at Q = Mqut- Note the break in the vertical scale in the figure. The ordering of 
the magnitudes of the off-diagonal elements of can be simply gauged from the up-type 
Yukawa coupling matrix, since all SSB mass matrices at the GUT scale are proportional 
to the unit matrix. These off-diagonal elements start from zero at Q = Mgut 5 evolve to a 
maximum magnitude, smoothly reverse direction at an intermediate scale and then continue 
to evolve monotonically all the way to the weak scale. The dips in the figure occur where 
the real, and simultaneously the imaginary, part of (m^) changes sign during the course 
of its evolution to Q = mt. 

Fig. 7.13 shows the running of rrig for the same mSUGRA point. In the left 
frame, we show the magnitudes of the elements in our "standard" basis, where the up 
quark Yukawa coupling matrix is diagonal at Q = mt. The frame on the right shows the 
magnitudes of the elements of this same matrix, but in the basis where the down- type quark 
Yukawa coupling matrix is diagonal at mt. We have checked that the large difference in the 
size of the off-diagonal elements in the two frames is indeed accounted for by the fact that 
the corresponding matrices are related by (3.13). Unlike in Fig. 7.12, there is no dip in the 
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Figure 7.13: The scale dependence of the magnitudes of the entries of mfj for the mSUGRA model 
with parameters as in Fig. 7.4 in the basis where the up-type (left frame), or the down-type (right 
frame) , quark Yukawa coupling matrix is diagonal at Q = mt ■ Note the we have broken the vertical 
scales to better the display the matrix elements. 



magnitudes of the off-diagonal elements because they evolve monotonically from zero at the 
GUT scale, until the squarks are all decoupled. 

7.5.2 Additional GUT Scale Flavour Structure 

To understand how the non-universal boundary conditions in (6.4b) impact the 
evolution of the squark mass matrices we examine a model with non-zero values of Rjj,d 
and Su^d, but with universal gaugino masses and a-parameters. In Fig. 7.14 we plot the 
value of (m^)23 for the set of values of Rjj,d and Su,d shown in the figure, with all other 
parameters set as in Fig. 7.12 (including qy^ = 1), so that Rjj,d = Sjj^d = corresponds 
to the mSUGRA model in this figure. 

We see that with non-zero values of Ru,D and Sjj.D: ( m ;y)23 already starts off with 
a substantial value (positive or negative) at Q = MguT; and evolves slowly with Q. This 
situation is qualitatively similar to that of the mSUGRA case shown in Fig. 7.12 once Q 
has evolved away somewhat from Mqut so that the matrix element has had a chance to 
grow from zero. However, because the curves start of with rather large values at the GUT 
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Figure 7.14: The scale dependence of (111^)23 for several sets of values of Rjj,d and Sjj,d that appear 
in (6.4b) for the model with non-universal GUT scale SSB squark mass parameters, but cjjjj = 1. 
The curves are in the same order as the legend, and all the other parameters are set as in Fig. 7.12. 

scale (except for the middle mSUGRA curve) the evolution does not take them through 
zero for any value of Q > writ. As a result, the dip which was the most prominent feature of 
Fig. 7.12 is absent, except in the middle curve which does cross zero for Q ~ 2.5 x 10 5 GeV. 

Finally, in Fig. 7.15 we return to the non-mSUGRA case that we considered in 
Fig. 7.7. We emphasize again that in this case and m? D are diagonal only in the basis 
where the superpotential Yukawa coupling matrices have large off-diagonal elements so that 
we would expect this model to include new and potentially large sources of flavour- violation 
that may well already be excluded by data. As in Fig. 7.12, we show | (m^)^. | in the basis 
where the up-quark Yukawa coupling matrix is diagonal at Q = mt- We see from the figure 
that the matrix has large off-diagonal entries at Q = Mqut- As expected, the gauge 
(and gaugino) interactions cause the diagonal entries to increase, whereas the off-diagonal 
elements which do not "feel" these terms evolve much more slowly with Q. Note also that 
because (m^) n < (m^) 22 at Q = Mqut, Yukawa coupling effects draw them closer as we 
go to low scales. We also remark that the negative GUT scale value of S tends to reduce the 
diagonal elements of mfj as we go to low scales, but pulls up the corresponding elements of 
m 2 D . 
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Figure 7.15: (a) The same as Fig. 7.12 except for the non-universal model considered in Fig. 7.7, 
and (6) the magnitude of the corresponding elements of m 2 D for the same scenario. The elements 
are plotted in our "standard" current basis but are exactly the same in the basis where the down 
quark Yukawa couplings are diagonal. 



Although the renormalisation group evolution increases the gap between the off- 
diagonal and diagonal elements at the low scale, notice that the off-diagonal elements are 
separated by just one order of magnitude from the difference between the diagonal elements, 
so that we may expect large flavour-mixing between the SU{2) singlet up squarks. This 
mixing, if anything, is even larger in the down squark sector as can be seen in frame (6). A 
careful evaluation of inter-generation squark mixing is clearly necessary for any discussion 
of flavour-violation in squark decays, the subject of the next chapter. 
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Chapter 8 

Decay of the Lighter Stop 



As we alluded to in the introduction, the sfermions of the MSSM have interactions 
which permit the t\ to decay according to t\ — > tZi and t\ — > 6Wj. However, the large 
top mass and Yukawa coupling leads to several factors which may preclude the lighter stop 
from decaying through either of these two options. As we have seen, in mSUGRA the large 
top Yukawa coupling depresses the (3, 3) entries of the SSB soft mass matrices rrig and 
mjy during RGE evolution from the high scale, rendering the lighter stop potentially lighter 
than the chargino. In addition, there is the possibility for a large mass splitting between the 
stop states coming from the a-parameters, further reducing the t\ mass. In a scenario where 
the t\ is sufficiently light, the two-body decays listed above will therefore be kinematically 
inaccessible. 

In order to apply our RGE analysis to flavour- violating decays of squarks, we must 
find the relation between the current basis squark states and the mass eigenstates. Using 
the weak scale SSB and 'Yukawa' coupling matrices, we will write out the mass terms in 
the Lagrangian and diagonalise this matrix to obtain the physical squark states. Rewriting 
the various squark interactions in terms of these physical states, we will then be able to find 
the partial width for flavour-changing stop decay [40] . 

8.1 Squark Mass Matrix 

As we saw in Eq. (3.16) of Sec. 3.2.1, the up-type squark mass terms can be 
collected into a (6 x 6) matrix so that the Lagrangian contains 

v ' \ {M 2 LR y aj (M 2 RR ) ab ) \ URh ) 
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where i,j label left-handed squarks and a, b label right-handed squarks. The physical states 
are found by diagonalising the (6 x 6) matrix, and the squark mass basis is defined to be 
(ii, t2, ci, C2, ui,U2) so that the t\ is signified by 

The individual entries of (8.1) were written in Sec. 3.2.1 and we repeat them here 
for convenience 

(ML) tJ =«) y . + vl{rX) i3 + (g - f ) (vl - vj) 8 tJ (8.2a) 

(M 2 RR ) ab =(™ 2 u) ab + v 2 u (%K) ab -^K- «2) ( 8 - 2b ) 

(ML), = - «« K)* 6 + « d (>f^)* b • (8.2c) 

where we now keep track of the possibility for separation between fi and ft, and f u and . 
These relations should be used with caution. Recall from the discussion in Sec. 5.5 that 
we have a restricted group of operators when Q < mn- If we are evaluating {M-\ R ) ib at a 



- S (a u )* 6 + c(A*f>); 



for which the RGE 



scale below mu we must use the combination 
is (B.21) in Appendix B. 

8.2 Calculation of the Width 

The two-body flavour-changing decay of the stop occurs due to non-zero couplings 
of squarks to quarks and neutralinos. These couplings appear in the Lagrangian from both 
gaugino interactions and superpotential terms and are 

C 3 ~ 75 {^ Q)jiXS + l^' Q) ^ X °) PhUi + ^(s' UR )ba U Ra P L X + h.C.j 

" {^o^(fr)^Ln l + n b (f u Q )^ Ll P L ^o +h.c.} . 

(8.3) 



Following the prescription in Ref. [6] , we replace the gauginos and higgsinos with the lightest 
neutralino using 
squared masses 



neutralino using the neutralino eigenvector, and inserting a (275)^ to allow for negative 



A =J>?W)^ (8.4a) 

i 

A 3 =5>?(*75)^ (8.4b) 

i 
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in order to obtain 



C (uLjcZij =u\-Z x (iAjPL - (if'^y^PR) c + h.c. (8.5a) 
C (ujftcZi) =^ Rh Z\ (iB b P R - (-if' )l 2 v?P L ) c + h.c. , (8.5b) 



where 



01-1 



(8.6a) 



B^ — ig'^i^-'v^ . (8.6b) 

Using (3.23) to write the squarks in terms of mass eigenstates and picking out the ti, we 
obtain the Lagrangian of interest: 



C (ticZi) 9 [aF L + (3P R \ c + h.c. (8.7) 

The coefficients of Pl and Pr are 

a s(Wl){,-A,- - {-i)^v^\u R )\ h {^) T b2 (8.8a) 
/? ^(W^B, - (^) el ^ 1) (W L ) t lJ (f M Q )* 2 • (8.8b) 

where we assign the first index of IA^ L R to be a = 1 since this is the index which represents 
the t\ eigenstate. With this Lagrangian we find that the partial width of the t\ to c plus 
Z x is: 

T(h cZj) = 1 ^ r {(l«| 2 + |/3| 2 ) (m| - m 2 - m|J - 2m c m^ (a/?* + /?«*)} 



ti 



x A 1/2 (m-^m^m^ 



(8.9) 

with 

A (a;, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz . 

Since the splitting between the diagonal entries of gW* and their non-tilde coun- 
terparts is small, the j,b = 2 contributions to Aj and B?, in (8.6) are qualitatively approxi- 
mated by using the standard gauge couplings. Similarly, the contributions from f* in (8.8) 
are qualitatively approximated by the usual Yukawa couplings. We have checked that the 
width calculated ignoring the tilde-couplings differs from the full calculation in general by 
a few percent. 



119 



Furthermore, we find that in the mSUGRA class of models, the contributions to 
r(ti — ► cZi ) from the threshold-induced off-diagonal entries of gW* are generally two orders 
of magnitude smaller than the contribution from squark mixing (the size of which we will 
estimate in Sec. 8.3). To see this, we first write the RGEs for the gaugino couplings in terms 
of the mass-basis Yukawa matrices, using the KM matrix to convert in our "standard" 
current basis into the basis where the down quark Yukawa couplings are diagonal at m t . 
This means that the majority of the flavour structure in the down Yukawa coupling terms 
will be contained in the KM matrix. On the understanding that is now in this basis, and 
taking all three flavours of sfermions in each group to decouple at the same point, we can 
write the seed term for off-diagonal entries to the SU(2) gaugino-squark-quark coupling, 
g*2, whose RGE is found in (A. 17), as 

d(g Q h 1 



dt 16vr 2 



(f u )* fe (f u )Ii(g%- + K ik (f d y kl (f d )f m Kl n (gQ) nj , (8.10) 



where i ^ j. For the off-diagonal entries this relation holds at the threshold of the heaviest 
SUSY particle, where (gP)ij = g x 11. To estimate the size of the off-diagonal elements of g^ 
we therefore take them to be diagonal on the right-hand side of (8.10) and choose k, /, m, n 
to find the maximum contribution. Referring to the general size of the off-diagonal elements 
of the Yukawa matrices in both Fig. 7.4 and the right-hand frame of Fig. 7.5, in the region 
of the highest SUSY threshold, we find that the largest off-diagonal entry is approximately 

4g Q )23 1 



[(1 x 10~ 7 ) (1)(0.1) + (4 x 10~ 2 ) (0.2)(0.2)(1)(0.1)] , (8.11) 



dt 16tt 2 

so that this entry can at most become ~ few x 10~ 6 , coming from the down-type Yukawa 
coupling contribution. The £7(1) gaugino couplings will have much smaller entries since their 
RGEs do not contain KM dependent Yukawa coupling terms. If we now allow for additional 
flavour-violating terms arising in the SSB parameters, in a non-mSUGRA scenario, off- 
diagonal gaugino couplings are still unimportant for the width calculation. This is because, 
even if the sfermion splitting in quite large, the effect of their additional flavour-violation 
will only enter the RGEs for the gaugino couplings through the thresholds. Although with a 
completely general rotation between the sfermion mass basis and the original current basis 
there will be additional large contributions to the RGE, the effect will be limited due to the 
small size of the logarithm and will therefore tend to remain sub-dominant. For example, 
although in the case of large mixing in mfj the off-diagonal entries of the U(l) gaugino 
couplings can become proportional to (f u )33, there will be corresponding increases in the 
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mixing terms LIl,r which will simultaneously increase the terms depending on diagonal 
entries of g'* . 

Turning to the Yukawa coupling contributions in (8.8), we see that since (Wl)^- 
and {U.R)\ b are largest when j, b = 3 (i.e., when the squark current basis state is II and 
Ir respectively), the Yukawa coupling contributions are largest for (f u )32, which is ~ 1CT 7 . 
For the j,b = 2 case, the mixing is reduced, and f u oc m c /M\y, so that the Yukawa coupling 
terms are unimportant in mSUGRA. Note that in non-mSUGRA scenarios, since the squark 
mixing can become large, these terms are no longer sub-dominant, and may need to be taken 
into account. 

In mSUGRA, we are left with the two gauge coupling terms, (8.6), as the dominant 
contributions to (8.8), so that the width is dependent upon the standard gauge coupling 
terms and the size of the squark mixing given by (14l,r) 12 from (3.23). To a good approxi- 
mation we can write the mixing as 

(1 , vt _ (Mu) 23 {U L )\ 3 + (Mu) 26 (Ur){ 3 , . 

{Ul)12 ~ ml - (M fi ) 22 ' (8 ' 12a) 



_ (M fi ) 53 [U L )\ 3 + (M fi ) 56 (U R )\ 3 
mf - (Mu) 



where the (Ul,r)\ 3 can be approximately calculated using the (2 x 2) sub-matrix for the 
mixing between ti and tR. 



8.3 Single-step Estimation of the Width 

As an alternative to the full calculation of the squark mixing outlined above, 
Hikasa and Kobayashi [53] developed a scheme for estimating the width without solving 
the RGEs precisely. Working within the mSUGRA framework with squarks of around 
30 GeV, and using the approximation that Z\ ~ 7, they showed that if tree-level two-body 
decays of t\ were kinematically forbidden, t\ — ► cZ\ would be the dominant decay mode 
of t\. While their approximations and analysis are certainly valid for values that they 
considered fifteen years ago, this is not the case for top squark masses in the range of interest 
today. Since their approximate formula is widely used in the literature when calculating t\ 
decay patterns [19, 20], we estimate the decay width using their approach (including the 
straightforward extension to general mixing in the neutralino sector [54]) for comparisons 
with our full RGE calculation. 
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For models with no squark mixing at the GUT scale, in the limit that m c and 
0, cr, ul and ur do not mix with the stops so that (for the t\), (14r)\ 2 vanishes and 
we only need to find (Ul)\ 2 - We write the mass matrix for the scalars exactly in terms of 
the quantities A^ and Ar which are the II — cl and Ir — cl entries: 
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We can find an approximate value for Al and A^ at the scale Q by approximating the 
RGE running to be a straight line and evaluating the right-hand side of the RGEs at a 
single scale. If all off-diagonal entries of the Yukawa and SSB matrices are ignored, A^ and 
Ar are given by 
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(8.15) 



Since these equations do not take into account threshold effects, we do not use different 
relations in the Q < vtlr regime. If we are calculating the stop decay at a scale below mn, 
we evaluate A^ and A# at Q = run- 

To derive these relations, we see from (8.13) and (8.2a) that A^ ~ (01^)32 if we 
ignore the contribution of the Yukawa and gauge coupling terms, so we must find the weak 
scale value of (m^^. To this end, the running of rrig is written down using a single step: 



m|(Q)^m^(GUT)-ln 



Mqut 



dm 2 Q 
dt 



(8.16) 



where we take Q to be somewhere near the weak scale. We choose to use our "standard" 
current basis so that the up Yukawa matrices are approximately diagonal at the scale Q. 
We then write the down Yukawa couplings in terms of the KM matrix using (3.11) so 
that we have fd(d), which is also approximately diagonal. In this basis we ignore all terms 
proportional to off-diagonal entries of the Yukawa and SSB matrices which we assume to be 
higher order corrections. To further simplify our equation for Al, the masses of the light 
quarks are set to zero, which makes both f u (it) and fd(d) identically zero except for (f U; d)33- 
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Under these approximations, the (3,2) entry of mg(Q) reduces to the form (8.14). The 
derivation of Ar is very similar except that the important quantity to run is the left-right 
coupling, (a u ) 23 . 

Next, we approximate the mixing to be a small perturbation to the eigenstates, so 
that the physical t\, I2 and cl states are only changed a small amount by the mixing. The 
mass matrix is approximately diagonalised to obtain 



(8.17) 
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where the small perturbations are found to be 

A* L cos 6 t -A* R sin 9 t 

^ 2 2 
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, A* R cos t + A* L sin t 
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m~ — m% 

The most important relation to take from (8.17) is that of the physical t\ state 

ti(phys) =ti + ec L . 



(8.18) 
(8.19) 



(8.20) 



The mixing factor, e, is therefore the extent of mixing of the t\ with cl — which we denoted 
in (3.23) by (Wl){ 2 , and estimated in (8.12a), the form of which is directly analogous to 
(8.18). In mSUGRA scenarios, has a size of approximately 

Mqut 



|a l | 



16vr 2 



In 



K* cb K tb m 2 fZ , 



_Msvsy_ 

and the Ir — cl mixing angle has a similar magnitude, so that |e| ~ \Ai\/(few rnfy ~ 
few x 10" 4 . 

Due to the approximations made in this method, the rate calculation in (8.9) 
simplifies dramatically to become 
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8.4 The Decay Subroutine, STDEC 



We pause here for a brief description of the numerical calculation of the two- 
body decay width, for which we have developed a subroutine called STDEC. As described 
in Chapter 6, the programme RGEFLAV runs the full set of RGEs for the dimensionless and 
dimensionful couplings from the high scale, which in mSUGRA is the GUT scale, to m#, the 
scale at which the heavy Higgs bosons decouple. It does this using an iterative procedure 
in an attempt to fix the boundary conditions at three scales, Mz, Msusy and Mgut- At 
the end of the iteration the couplings are rotated into the basis where the up-type quark 
Yukawa couplings are diagonal at and then passed to STDEC. 

STDEC evolves the couplings to the scale at which the decay is to be calculated, 
which we choose to be the lightest of the eigenvalues of the left- and right-handed soft mass 
matrices for the up-type squarks, using the subroutine DECRUN. When the programme arrives 
at the required scale, the mass matrices are reconstructed using the saved eigenvalues and 
eigenvectors from the running procedure. In other words, we recreate the mass matrices 
in the basis where the squarks are diagonal using the eigenvalues of each squark at its 
decoupling scale. We then rotate from the squark mass basis to the "standard" current 
basis using the known eigenvectors in this basis. 

Once we have all the couplings at the correct scale in the quark mass basis, we 
construct the (6 x 6) mass matrix using (8.1) with (8.2) in the subroutine CUSQM — re- 
membering to use the restricted set of couplings if we are calculating the decay at a scale 
Q < mu- The mass matrix is then diagonalised in USMMA, by the EISPACK subroutine CG, 
to find the eigenvalue and eigenvector of the t\ . 

Finally, the decay rate calculation from (8.9) is carried out, unless the kinematics 
forbid this decay from occurring. To find the rate, we must know g"^, g^, g' UR , f$, and f" fl , 
along with the masses and mixings of the squarks and neutralinos. The result is compared 
to the estimate obtained using the single-step estimate in (8.21) and the partial width for 
the three-body decay [19], obtained by ISAJET. 

^If the user has chosen for the output to be in the basis where the down- type quark Yukawa couplings 
are diagonal at m t , STDEC exits without calculating the partial width. 
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Method 



Width (GeV) 



Y{tt -► bWZt) 



'single-step" estimate 



8.6 x 10- 



41 x 10- 



8 



8 



full calculation 



3.3 x 10~ 8 



Table 8.1: Partial widths for the two- and three-body decays of the t\ for the mSUGRA parameters: 
too = 250 GeV, to 1/2 = -250 GeV, A Q = -930 GeV, tan/3 = 20 and fj, < 0. The two-body decay 
width is calculated using two methods, a 'single-step' approximation, and our full RGE solution. 

8.5 Sample Numerical Results 

8.5.1 Single-step RGE Integration and the Stop Decay Rate 

As a first example, displaying the need for a complete solution of the full RGEs, 
we show results in Table 8.1 for the mSUGRA point mo = 250 GeV, m 1( / 2 = —250 GeV, 
A = -930 GeV, tan/3 = 20 and \i < 0. For this point, m~ ti ~ 181 GeV, = 102 GeV 
and W\ = 197 GeV so that tree level two-body decays of the t\ are kinematically forbidden 
but both the two-body decay that we are considering, and the three-body decay to bW Zi 
[19] are open. We see that the 'single-step' approximation of the decay width overestimates 
the actual value by a factor 13.6, the difference being entirely due to the error in calculating 
the value of e (= (jUl)\ 2 )- In the mSUGRA model with a neutral LSP, when the usual 
flavour-conserving two-body decays are kinematically forbidden, we find that the single- 
step approximation overestimates by a factor of between 10 — 25. 

If we used the single-step RGE approximation to find the branching ratios for this 
point we would conclude BR(cZ^) ~ 0.83 and BR(bW Zi) ~ 0.17 so that the three-body 
decay is subdominant. In contrast, when we use the full RGE calculation the three-body 
decay becomes dominant with BR(bW Zi) ~ 0.72, completely changing the qualitative 
picture of top squark decays. Admittedly, the change in the qualitative picture in this case 
is due to the three-body decay being close to the kinematic boundary. However, if the 
three-body decay is closed, the two-body loop decay will compete with various four-body 
decays, in which case a similar rearranging of the dominant decay modes may result from 
our RGE calculation of the two-body width. 

Note that the width obtained for the two-body decay depends on the choice of 
scale at which we calculate it. Since the t± is mainly we evaluate the decay rate using 
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the parameters at a scale equal to the lightest right-handed squark threshold. We find that 
if we instead use the scale where the heaviest right-handed squark decouples as opposed to 
the lightest, there is just a 6% change in the partial width. 

8.5.2 Model Dependence of r(*i -> cZ t ) 

To illustrate the effect of non-universal GUT scale inputs we first return to our 
sample point from the previous chapter, namely mo = 200 GeV, m 1 / 2 = —400 GeV, Ao = 
—200 GeV, tan (3 = 10 and \i > 0, and include non-zero Tq or Tu,D (while simultaneously 
setting mQ = or n?< v D y = respectively) in a variety of scenarios. Our results are 
shown in Table 8.2, beginning with Scenario (1) which is our standard mSUGRA case. In 
mSUGRA, all dependence on the ^Vl,r(u,(1) matrices is restricted to the KM combination, 
and therefore the result is independent of our range of choices in the other scenarios. We 
mention that the scenarios we examine are unrealistic both from the perspective of observing 
t\ decays — indeed, for this sample mSUGRA scenario, tree-level decays of t\ are accessible 
— as well as for that of flavour-changing neutral currents. Our purpose here is to understand 
how r(ti — > cZ\) is altered (and also which quantities it is sensitive to) as we change the 
scenarios to systematically allow increased non-universality and/or flavour- violation in the 
SSB sector. 

As expected for Scenario (2), although there is a significant splitting in the squark 
masses, the stop width is of the same order of magnitude as the mSUGRA result. We saw 
in Sec. 7.2.2 that for non-zero Tu,D entered in the "standard" current basis, there is little 
difference to the running of the usual Yukawa couplings and their tilde counterparts, and 
the same can be said for non-zero Tq in this basis. The width sees only a small effect 
because the RGEs for the gauge and Yukawa couplings only depend on the squark masses 
through the thresholds, and also because the rotation to the squark mass basis is small. t 
Any deviation from mSUGRA shows up in our diagonalisation of the squark mass matrix, 
and the difference between the first two scenarios is due to a ~ 15 — 25% change in the 
entry (Mi){ 2 . Note that although the eigenvalues of the squark mass matrices at the GUT 
scale have large splitting, the matrices themselves remain close to diagonal all the way to 
the GUT scale. 

^We emphasise that the situation would be quite different for flavour-violating decays of down-type 
squarks, since the down-squark matrix is now not aligned with the corresponding Yukawa matrix. 



126 



Scenario 



Width 



(1) 


mSUGRA — no dependence on specific Vi^fu, d) 


2.2 x 


10" 


-9 


GeV 


(2a) 


Vfl(u) = Vfl(d) = V L (u) = 11 




Tt/,D / 


3.9 x 


10" 


-9 


GeV 


(26) 






Tq/O 


1.6 x 


10" 


-9 


GeV 


(3a) 


V L («) ± 11, V fl («) = Vji(d) = 


n 


Tu,d ¥= 


3.9 x 


10" 


-9 


GeV 


(36) 






Tq/O 


2.7 x 


10" 


-5 


GeV 


(4a) 


Vfl(d) / ii, Vfl(«) = v L («) = 


n 


T[/,d + 


3.6 x 


10" 


-9 


GeV 


(46) 






Tq/O 


1.6 x 


10" 


-9 


GeV 


(5a) 


v fl («) / ii, Vfl(d) = v L («) = 


n 


T[/,d / 


5.8 x 


10" 


-3 


GeV 


(56) 






Tq/O 


1.6 x 


10" 


-9 


GeV 


(6a) 


Vfl(«) / ii, Vfl(d) / n, v L («; 




Tu,d + 


5.8 x 


10" 


-3 


GeV 


(66) 






Tq/O 


2.7 x 


10" 


-5 


GeV 



Table 8.2: A comparison of the two-body loop decay widths for six scenarios. For each scenario we 
list the basis used for our GUT scale inputs, with rotation matrices as specified by (7.2) in the text. 
In case (a) of each scenario, we take T Q = and T U>D = <fia5{10000, 40000, 90000} GeV 2 . In case 
(6), Tu. D = and T Q = diag {10000, 40000, 90000} GeV 2 . When T Q>U>D + 0, the corresponding 
TO {<2 u D}o = ^ anc ^ tne m SUGRA case is the same as in Fig. 7.4. The value of is constant in all 
scenarios to within about 5%, so its contribution to the variation of the partial width is small. 
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We can understand the effect of varying V L,R.(u,d), by considering rotating back 
from this general basis to our "standard" current basis. In effect, varying our choice of 
rotation matrices changes the GUT scale boundary conditions when viewed in our "stan- 
dard" basis. We would therefore only expect to see a change in the width for GUT scale 
inputs that correspond to a change when viewed in this basis. For example, in Scenario 
(3a) we have chosen a basis such that mfj D are unchanged when rotated back into the 
"standard" basis. Therefore entering Tjj.d / in the basis specified in Scenario (3) is the 
same as entering the same Tu,d in the "standard" current basis, i.e., Scenario (2), and we 
would expect this to be reflected in the resulting width. On the other hand, the boundary 
condition for rrig in the current basis does depend on V^(u), and we would therefore expect 
the width to be very different between Scenario (26) and Scenario (36). 

Looking, therefore, at Scenario (3) we see no change as a result of the Tu,D / 
input, but a four order of magnitude change as a result of setting Tq / 0. As already 
mentioned, since the GUT scale boundary conditions for mfj D have no dependence on 
Vi(u) or V^(d), but nig does, this difference in behaviour between (3a) and (36) was 
expected. The large jump in the width for (36) can be understood as a result of a dramatic 
increase in the level of t\ — cl mixing. We find that the Tq / input results in a change in 
both m.Q and mfj. This is because the seed terms at the GUT scale for flavour off-diagonal 
entries of mfj depend on the structure' of the matrix mg. However, although the change 
in the off-diagonal elements of mfj is in this case many orders of magnitude, the effect is 
still smaller than the increase in t\ — cl mixing coming from m^. 

Moving on to Scenario (4), the rotation matrices used in this case mean that the 
boundary condition for Scenario (4a), the Tjj,d 7^ case, is different only in m^, affecting 
the running of the rriQ v matrices to a small extent through the RGEs. As a result, we see 
that the width is different from the other scenarios by a small factor (< 10%). On the other 
hand, the boundary condition on rrig is not dependent on V^(ti), and the width in (46) is 
therefore unchanged from Scenario (26). 

In Scenario (5a), where we choose to enter the GUT scale matrices, T?u,D) in 
a basis which depends on V r(u) (thereby altering the GUT scale boundary condition on 

^The RGEs in our "standard" current basis contain KM factors next to down-type Yukawa couplings 
only. Therefore, if rtiQ is diagonal in the "standard" basis, we must use off-diagonal Yukawa coupling terms 
to obtain off-diagonal seed terms in the mfj RGE. On the other hand, if rriQ has all entries 0(1), the 
off-diagonal seed terms depend on diagonal Yukawa matrix entries. The same goes for in the RGE. 
In contrast, m 2 D always appears with KM terms in the m 2 q RGE in the "standard" basis, and therefore we 
can obtain off-diagonal terms from diagonal Yukawa matrix entries regardless of the form of m 2 D . 
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mjy), the largest contribution to the width comes from the fact that the mixing between the 
right-handed squarks has been greatly increased. In particular, the cr forms a substantial 
fraction of the t\ and couples strongly to the bino-like neutralino. This is due to the fact 
that the off-diagonal entries of mfj (expressed in our "standard" current basis) are non-zero 
at the GUT scale and remain large during the running to the weak scale as can be seen in 
Fig. 7.15. Although one might expect the contribution to (8.8) from off-diagonal entries in 
the higgsino and gaugino 'Yukawa' matrices to become significant, the enhancement that 
they receive is not sufficient to bring them to the same level as the contribution from the 
(JUr)\ 2 ^2 term in [3 which is now two orders of magnitude larger than any other entry in 
(8.8). In contrast, Scenario (56) shows no change from (26), due to the fact, once again, 
that the boundary condition is the same in this basis. 

Finally, Scenario (6a) is equal to Scenario (5a), with a small difference not shown 
in the table due to the effect of m 2 D entering the RGE for uiq. Scenarios (66) and (36) 
coincide because when Tq / the GUT scale boundary conditions only differ for different 
values of Vl(-u, d). 

In summary, we find that when we take a large amount of splitting between di- 
agonal entries of the SSB scalar mass at the GUT scale, our results are very sensitive to 
the choice of basis. Although it is possible to choose the rotation matrices carefully enough 
that the effect can be small, as in Scenarios (2a), (26), (3a), (4a), (46) and (56), when the 
rotation is completely general as in the two final scenarios in Table 8.2, there is the potential 
for large deviations from the mSUGRA result. These deviations are overwhelmingly due 
to an increase in the relevant mixings, rather than a result of the difference between the 
tilde-couplings and their counterparts. 

We move on to take a small splitting between the diagonal elements — say, 
(mo ± e GeV) 2 — so that we may consider flavour- violation effects in a more realistic sce- 
nario. For this case we take the compressed SUSY scenario [55, 56], proposed by Martin, 
where efficient neutralino annihilation to top pairs via the exchange of a light squark leads 
to the observed cold dark matter relic density. We use mSUGRA-like GUT scale inputs, 
where mo = 500 GeV, Aq = Mi, tan [3 = 10 and n > 0, but the gaugino masses are split so 
that 1.5Mi = M 2 = 3M 3 . 

Fig. 8.1 examines the effect of increasing the size, e, of the splitting between the 
diagonal entries of Tu,D, where Tjj = = diag{(500 — £) 2 , 500 2 , (500 + e) 2 } (again setting 
m WD\o = 0) f° r a fixed value of Mi = —560 GeV. We see that the effect of the general 
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M 1= -560 GeV, A =M,, m () =500 GeV, tan|3=10 
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Figure 8.1: Effect on the two-body decay width of setting m^ fl = and = To = diag {(500 — 
e) 2 ,500 2 , (500 + s) 2 } in the general basis given by (7.2) in the text. We use the "compressed SUSY" 
GUT scale inputs where m = 500 GeV, M x = -560 GeV, A = M 1: tan/3 = 10 and fi > 0, and 
\.5M\ = M 2 = 3M 3 . The width is plotted as a function of e. 

rotation, (7.2), on the size of the width rapidly increases once e is given a non-zero value, 
and that even a small splitting of e = 1 GeV increases the width by more than four orders of 
magnitude. This is not the case if we enter such a small splitting in the "standard" current 
basis, where there is no additional flavour structure introduced. Indeed, the significant 
enhancement of the width is entirely due to our choice of rotation matrices Vl,r- 

To see how the model dependence of t\ — > bW Zi can have significant phenomeno- 
logical implications, we again consider the compressed SUSY scenario of Fig. 8.1 and vary 
\Mi\ up to 640 GeV. The result is shown in Fig. 8.2. Over the whole of this region the 
two-body flavour-conserving decays of the stop, as well as its decays to sneutrinos, are kine- 
matically forbidden, but the three-body decay i\ — > bWZi is allowed until it reaches its 
kinematic limit at approx. 640 GeV. The black dashed line in the figure is the partial width 
for the three-body decay which clearly develops a phase space suppression towards the high 
| Mi | end of the graph. The partial two-body width calculated using the single-step method 
is in the range 10~ 9 GeV and competes with the three-body decay only on the extreme 
edge of phase space. It is interesting to note, however, that the result using full RGEs 
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Figure 8.2: Comparison of the partial width of the three-body decay of the stop with the two-body 
flavour-changing decay in a "compressed SUSY" scenario. The GUT scale inputs are the same as 
Fig. 8.1 except that we now vary Mi and compare the partial width obtained with Tjj.d = to 
that obtained with raf JD = and Tjj.d — diag{(498) 2 , 500 2 , (502) 2 } in the general basis given by 
(7.2) in the text. The dashed black line shows the three-body decay t\ — > bWZi while the other 
lines represent the various two-body calculations. 



(which is the lowest line at approx. 10~ 10 ) is even less important in the branching fraction. 
Although this difference between the single-step result and the full result is unimportant 
for the three-body decay, it may be important in the case that a four-body decay begins to 
compete at large \M\\. 

The highest horizontal two-body line in Fig. 8.2 represents the result of our com- 
plete calculation taking Tjj,d = diag{498 2 , 500 2 , 502 2 } and m 2 ^ D -. Q = in the general basis 
of (7.2). We see the expected enhancement from taking a small splitting between diagonal 
entries in a general current basis, as predicted in Fig. 8.1. In this case, the two-body decay is 
competitive with the decay t\ — > bW Zi over a much larger range than before, representing 
the phenomenological implications of introducing dependence in the RGEs on the specifics 
of the V l,r{u, d) matrices. We should also mention that if we take this same splitting, but 
use ~Vr(u, d) = Vx(it) = 11, we recover the result with no splitting to within around 4%. 
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Chapter 9 



Summary 

One remarkable property of supersymmetric theories with widely disparate mass 
scales is that the hierarchy between the scales is stable to radiative corrections. Thus, we 
can construct models based at the high scale and legitimately calculate their implications 
at the weak scale. The RGEs provide us with the means to use our high scale theories to 
evaluate weak scale observables, and as a result of renormalisation group evolution, simple 
models can lead to a rich structure which will soon be probed at the LHC. 

For realistic models in which SUSY is softly broken by an as yet unknown mecha- 
nism somewhere below the Planck scale, the fact that sparticle masses are different to their 
SM counterparts comes with some intriguing consequences. In a purely supersymmetric 
theory, the dimensionless couplings of gauginos and higgsinos to the fermion-sfermion sys- 
tem are, respectively, related by supersymmetry to the three gauge couplings and to the 
three Yukawa coupling matrices that specify the interactions of matter fermions with Higgs 
bosons. These relations, however, are no longer valid once SUSY is broken. 

We have reexamined the complete set of RGEs for the MSSM (excluding the 
quartic scalar couplings), keeping track of threshold corrections and flavour-mixing effects 
from quarks and squarks. To implement the decoupling of heavy particles we consider a 
series of effective theories, where the particle content is dictated by the scale, so that the 
influence of individual particles on renormalisation group evolution is removed at the scale 
of their mass. It is found that once the heaviest SUSY particle decouples the evolution 
of gaugino and higgsino couplings is no longer the same as that of their SM counterparts. 
Indeed, the gaugino couplings to fermions and sfermions develop a matrix flavour structure, 
something forbidden in unbroken supersymmetry. We have also described how to treat the 
new complications to squark decoupling due to flavour mixing among the squarks, and to 
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the electroweak symmetry breaking conditions (traditionally used to fix the magnitude of 
the superpotential fi from the observed value of Mz) from the fact that below the scale of 
SUSY breaking it is not /j, 2 but a related parameter ft, 2 that enters the scalar potential in 
the Higgs sector. 

Our complete set of MSSM RGEs can be found in Appendices A and B. We do 
not list the RGEs for quartic couplings as these are less important from a phenomenological 
perspective, though they do enter the RGEs of the SSB mass matrices and trilinear couplings 
at one-loop. In this case, we set them equal to their supersymmetric value. As mentioned 
above, the RGEs contain a number of 'Yukawa' couplings in addition to those usually 
written for the MSSM. For these couplings we use the boundary condition that they are 
equal to the corresponding SM couplings at all scales above the highest SUSY threshold. 
We write the dimensionful couplings in two regimes — above and below the scale of the 
heavy Higgs bosons, run- This is because we have a restricted set of couplings and mass 
terms below run, where the heavy Higgs bosons decouple, and some linear combinations of 
couplings do not appear in the effective theory. 

The full set of threshold RGEs have been assembled into a computer code, RGEFLAV, 
which also includes two-loop MSSM and SM RGEs. Since we expect the effect of particle 
decoupling to be of the same order as two-loop sized effects (assuming 0(1) splitting in the 
sparticle sector), we have placed these two contributions to the running on an equal footing. 
Our inputs for SSB parameters at the high scale consist of two pieces: the most general 
form for SSB parameters that does not introduce new flavour- violating structure (i.e., all 
flavour structure arises from the usual Yukawa coupling terms in the superpotential), and 
additional terms that allow the use of arbitrary values for the SSB parameters in a general 
current basis. The RGEs are then solved using an iterative procedure. 

Sample results have been presented for a range of SUSY scenarios. We have shown 
the key features of the introduction of thresholds into the Yukawa and Yukawa-like couplings 
of Higgs bosons, higgsinos and gauginos to fermions and sfermions for a simple model where 
the SUSY thresholds are clustered at two distinct scales. We were able to verify our assertion 
that threshold effects are of the same order as two-loop corrections and also see that flavour 
off-diagonal entries in the gaugino couplings can be of comparable order or even larger than 
the corresponding Yukawa coupling terms, making the inclusion of these effects essential to 
obtain two-loop accuracy of results in the study of models where the Yukawa matrices are 
the sole source of flavour-violation. 
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These conclusions apply equally well to more realistic scenarios such as the sample 
mSUGRA point we have considered, where we have shown the dependence on scale of the 
SSB trilinear and scalar mass matrices. By introducing non-universal GUT scale boundary 
conditions, we found that it is easy to obtain weak scale SSB scalar mass matrices whose 
off-diagonal terms are less than an order of magnitude below the diagonal entries. How- 
ever, since the Yukawa couplings only "feel" the effects of the SSB parameters through the 
thresholds, they are not greatly affected in the region of the scale of the SUSY masses. 

We have also considered the renormalisation of the gaugino mass parameters, and 
possible deviation (in addition to two-loop effects) from the usual MijM\ = 02/01 relation 
expected in models with gaugino mass unification. In mSUGRA, we find that the thresholds 
have an effect on this relation which is comparable to the two-loop contribution. This would 
be particularly important for any quantitative analysis at a possible e + e~ linear collider. 
Moving on to consider a more extreme model of SUSY, we find that when the particle 
spectrum is significantly separated out, as in the split SUSY scenario, there is the potential 
for violation from the MijM\ condition at the ~ 20% level. 

The analysis presented here is intended only as a sample of the range of effects 
that, to our knowledge, have not been previously considered. The introduction of flavour 
structure in the gaugino couplings could be of particular importance in future studies of 
flavour-violation in the sparticle sector. In addition, the ability to introduce non-universal 
high scale inputs, while maintaining strict control over new sources of flavour- violation could 
be useful for model building. 

As a sample of the phenomenological implications of our work, we have considered 
the two-body flavour-changing decay of the top squark, t\ — ► cZ\. We found that previous 
estimates for the partial width in mSUGRA have overestimated by between a factor of 
10 — 25 so that in regions of parameter space where the two-body decay may compete with 
three- or four-body decays, our analysis obtains a much smaller branching ratio for the 
flavour-violating mode. However, we have seen that an introduction of even small amounts 
of additional flavour structure at the GUT scale can enhance the two-body partial width by 
many orders of magnitude, resulting in an even greater competition with the kinematically 
suppressed modes. 

As anticipated by the results presented in Ch. 7, the decay rates in Sec. 8.5 are 
sensitive to the individual matrices, \~L,R,(u,d), that enter via the diagonalization of the 
Yukawa coupling matrices, and not just to the KM matrix. Indeed, this dependence of 
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physics on the separate matrices is the generic situation, while the dependence of physics 
on just the KM combination that we have become used to from studies within the SM or 
the mSUGRA frameworks, is true only in very special situations: see Table 8.2. t 

In conclusion, we find that a prerequisite to a full two-loop level treatment of the 
RGEs is the introduction of threshold effects such as those presented here. The difference 
between the usual gauge/ Yukawa couplings and their SUSY counterparts, which results 
from the decoupling of SUSY particles from the effective theory, is of importance in a 
wide range of phenomenological scenarios. In addition, threshold effects in models with a 
large misalignment between the SSB matrices and Yukawa coupling matrices can result in 
large changes to the off-diagonal entries of the Yukawa coupling matrices at the high scale. 
We have presented the full set of RGEs for the MSSM (excluding quartic couplings) and 
observed that a consistent treatment of flavour effects in the RGEs can make a significant 
difference to the widths of flavour- violating decays of supersymmetric particles. 



' Although this is well-known to many authors, we stress this here because there has been occasional 
confusion about this issue. For a different example, see p. 215 of Ref. [6]. 
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Appendix A 



RGEs for Dimensionless Couplings 



This appendix contains the RGEs with full thresholds for the dimensionless cou- 
plings of the MSSM with i?-parity conservation. Note that these RGEs are valid only in 
the basis in which the SSB sfermion mass matrices are diagonal. In any other basis they 
must be modified to account for the rotation from this basis to the sfermion mass basis so 
that the sfermions can be properly decoupled. For more information, see Sec 4.4.1. 
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Appendix B 



RGEs for Dimensionful Couplings 



This appendix contains the RGEs with full thresholds for the dimensionful cou- 
plings of the MSSM with imparity conservation. Note that we write these RGEs in the 
current basis in which the SSB sfermion mass matrices, but not the quark Yukawa matri- 
ces, are diagonal. In any other basis they must be modified to account for the rotation 
from this basis to the sfermion mass basis so that the sfermions can be properly decoupled. 
In this case the squark Oqjs become matrices ® q as discussed in Sec 4.4.1, where further 
details may be found. The RGEs for the superpotential parameter fx and the gaugino SSB 
mass parameters are, 
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The following RGEs are only valid above Q = run, where 9h = Oh = 1- We 
separate the two regimes of different Higgs boson content to simplify the resulting formulae, 
and to make explicit the parameters which remain in the theory below the heavy Higgs 
decoupling scale, Q = rein- 
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+ ^H>(g Q )I fc (g Q )H + |eg(g?)*(g?)«] (a«)« 

3 {(^ + 9^ + 4) 5,2 + 4(^ + 1 )^ 2 

+ ^( <9 Q i + ^)si}( a «)«. 
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(4tt) : 



dt 



9 + 3 













ik 



tfc 



+2^(^«)l,(^«)y 



+ 



3(fJ)«(f < i), fc + (f e t ) fc i(t)jfc+2 (9 B 



2 3 
+ -< 



' w 



-J l-d 



'J 



8 

5' 



- 3 



4 
3 
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, 6 2 

12 



®di ®Qk 



$ik +4 

,2 



ik 



2(f^-^)« J + 6(f^(f d )L 



+ 2^^(f^(fe)L(ae)fe/ 
--^-°kj + 2 



(f d ) f (fd) 



kj 



+ (M X - *M() (-^(^)*(g'°)* fc (f^)^ + 3(g ,Q )* fc (f d ) fc Kg^)^ 
-2^(f«), fc (g^)^.(^)*) 

- w w e~ h (M 2 - iJU*) (^ d )*(g Q )* fc (f^) fcJ 



32 



g (Ms - iAfS) (g?)Sfcte)H(g*% 

w Q k(f u Q )L + ^(f?)^(f?)L + A^(g /Q )i fc (g /Q )S 



'hV-u )ik\L u )kl^ v h\"-d )™V-d >kl ' jg - 

+^vK(g Q )i fc (g )H + ^Kg?): fe (g?)S 



(ad)jj 



+ 



3(fd) M (f]k + (fe)^(f e t )^ + ^% 



+ 2%A 



(ad) 



'j 



OB(s >dR )li(s ,d % + g^(g?«)£(g?«)|;- + 2^(ff )£,(§*)y 



+ &d k ( a d)ik 
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db 



V —9 a a I r,*f h 



2l 9 - 



.12 

18 



ik 



., ; S ik Sij + 6({ d )ij(fd)l k 



kl 



kl 



+ e 



Qk 



g' 2 30?'. . 

— + — ) S ik -2 
6 2 



(f u )(f u )t] J (>f^) 



-2^ (/i%"") ifc [(f«) f (fd) 



3^Vf" (2(3W*% + (g'°)* fe (ff )^ 



+ 



+ 



7 W 



kl 



3 „ ,.n,. /.n,r 8 



hi 



-3 



36 Ql 



9 d J 4 



,o 3 

5 + 4 



4 



(%+^>3 2 }K d ) u ., 
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(4 )2 djaehj 
v ' dt 



g' 2 zgl 



S tk +4 [(f e )(f e ) t 



ik 



( a e)fcj 



+ Mz t [-g' 2 SikSij + 2(f e )^(f e )L] (a e )n + 6^-^ fc (f e ) l ,(f,)L(a d ) M 



+ #e fc (a e )ifc 



V 2 «5 jfcj + 2[(f e )t(f e ) 



kj 



-29^^)10^)*) 
- Q9 w e- h (M 2 - iM' 2 ) Cg h ^(g L )U^) kj 

( a e)/j 



^(fe L )^(f e L )L + ^B^Mi'^ + 2^(g L )I fc (g L )S 



+ 



3(f ci ) H (f])a. + (fe) fe Kfe t )^+2^ 

+ 6e k ( ae ) lk [26^)1^% + 2^(f e e «)L(f e e «) / , 
" 3 { (l^ + ^ + |) 5 2 + \ (0 U + l) 52 2 } (a e)i , 



(a e )jj 
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( 47r ) n =5 #e fc %' ^ f e d J 



+ %^ fc [- 5 ' 2 ^^- + 2(f e ) l ,(f e )t fc ] (>f e ^) 



A;/ 



g' 2 3,| 

"2 2~ ' ife 



-26 s e- h ^(2(f- )u 

)ik 



)ij - {$ L )un R )k 3 ) 



-Q9- h e w ^~g h ^e)* k {i e e R )k 3 



(B.13) 
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+ 



3&)«(f!)i fc + ^ + l&w \~9 hu \ 2 } 



■ 3 {(^ + ^ + i)^ + i(^ + 1 )4( A ^)«' 
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+ Oujut [-2g'H lk + 6 [(iu) T (iuT] lk ] {m 2 u ) M 

+ 6 qM [^ fe + 6[(f u )*(f u fy K) H 
+ ^9%k {™l) kl - o Lk e il9 'H lk (mi), 



+ 60^ |/f (&)U$ d )kl + |jS| 2 {fr)U^ T - 



kl 



+ 3^ 



-29- h \6 s (Mf + Mf) 



+ 3^ (Mf + M2 ) 



3</ 2 , 9<£ 



-Jh,, 



+ 30^^ 



ma 
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d ( m 



(47T) 



h + l^l 2 ) 



-</ 2 (™k + lA| 2 ) + ^[</ 2 + <7 2 2 ] (m 2 Hd + m 2 ) 
+ 29u k 9u l g' 2 5i k {ml) kl 

+ e Q k 6 Q l [-9 ,2 sik + e[(f d y(f d f] lk ] (my, 

+ <Vl, [<7 /2 ^ + 2[(f e )*(f e f]J (mi) H 
+ gfc g| + 2 [(fe) T (fe)*]J (m|) fci 

+ 60^, |/i| 2 (**»)r*(fj")£ + e^^f^s 
+ 2^^(a e )r fe (a e )^ 



'kl 
'kl 



Jh d 



39 W 



--,h d 



-20- h \d § {M 2 + M[ 2 ) 



-Jh d 



3d w {M 2 + M' 2 ) 



--.hd 



3^ H 
2 2 



)(«4 + iai 2 ) 

+ {[6f d *fJ + 2Cf e T ] fefc + ^ 



-Jh d 



36 w e h 
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I 5 % + 2[(f u r(f u f]..|(m^ + |A| 2 ) 
1 - g % j + 2[(f d )%f d )% j }(m% d + \tf) 



~ fuJ 6ij { mu ) kk + -O^g' 8ij {ui Q ) kk 

-' 2 zgl , 8 9 l 



Qk Qi \ 



9 



5 ik 5i j (m Q ) 



+ 20^(f u )* fe (m 2 ,)^ (f„)£ + Wtf^UY* {m 2 D ) kl 
+ 2fl 4 (a d )S k (a <( )£. + 20 4 |/2| 2 (f>)* fc (f>)^ 

- 2 9 9 S {Ml + Mf) (g' Q ), fc (g' Q )l, " ™ w (Mi + <) (g% k (E Q )l 

- f % (M| + M'i) (g?)*(g?)t. 

- 46fc |^| 2 [(f«)* fc (f^)£ + (f, Q )I fc (f rf Q )£ 

- 3 (^ + ^)(^ /2 +^ + ^3 2 ) W) 



kj 



'I] 



^Mg /Q U(g ,Q )L + §Ws g )*(8°)L + |^(g?)*(g?)L 



+lo~ 9 (sf) k i(g% + 0^)1^)1 + e~ h ^)u^)l 
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4 

~3 ; 



<?% +4 [(f u ) T (f u )%j (m^ + + ^% (ml d + 



</ 2 + sl 



kl 



8 8 

+ ^U k 9' Sij (mu) kk + g^*^! 



0e fc <?% (m|) fcfc + 4^ fc ^(f u )f fc (f u )r j (m£) 



fci 



u )kj 



on Q9 

- y ^ ( M i + O (s /Ufl )L(g' UH )^ - y *s (^1 + M'i) (sHk(g R hj 

- 89- h |^| 2 - 3 + fij ) (^' 2 + i 5 2 ) (m 2 ^ 



+ [g^(g /Ufl )L(g ,Uii )^ + 3^(sDL(gr)^ 

+2^(^«)£(fj«)5-' 
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g' 2 5 tJ (m 2 Hu + |/2| 2 ) + {-^% + 4 [(U) T (U)%} (m 2 Hd + |£| 2 ) 
4 2 2 



2 

S ik Sij (m 2 D ) kl - -O^Sij (m 2 L ) 



>kk 



+ ^(a^a^ + A6 Qk |/f (f^)E(f£%- 

q on 

- -0 6 [Ml + Mf ) (g' d «)l k (g' d *) kj - -e~ g (Ml + M'i) (sH k (S dR )kj 



~ 80 h |„| 2 - 3 + ^ j ( V 2 + \g 2 ) (m 2 ,) 



1 72 , 4 2 



J, 



9 3 3' 



+2^)f fe (f^ (m 2 D ) 
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" 9'% (rn 2 Hu + |/2| 2 ) + { 5 % + 2 [(f e )*(f e ) T ] ..} (m 2 Hd + |/i| 2 ) 

+ 2^% (m 2 ^ - e Qk9 '% { m %) kk - e ik9 '% 3 ( m 2 D ) kk 
+ e ik g'% (ml) kk + e ik e Li + 3 f 2 ) Mj K ) H 

- 0e fc <?% (m|) fefc + 20 e - fc ^(f e )* fe (m|) fe/ {U% 
+ 20 e ~ fc (a e )* fe (a e )£. + 20 e ~ fc |/2| 2 (^)* fe (f e ^ 

- 20 fi (M 2 + M[ 2 ) (g' L k(g% - 6^ (Mf + M^ 2 ) (g L k(g% 



- 4^ |M| 2 (^(^S " 3 {e U + + \gl) (m|) 



+ 



+ Zfe (m 2 ) . fc [^(g /L )M(g /L )!/^ + ^(g i )«(g L )!,v + 
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=2 5 % (m 2 Hu + |/2| 2 ) + {-2 5 % + 4 [(f e ) T (f e )*] ..} (m 2 Hd + |/I| 2 ) 

" MuJ%j M) fe * + 20Q k 9% {™ 2 Q ) kk + 20 4 <7% ( m ^ 

- 20 Zfc5 % (mi) M + 26, k g'% (m|) fcfc 
+ 2^0*0%** (m|.) H + ±e Lk e Li (f e )l (mi) H 
+ 4^(a e )E(a e )^. + 4^ |/I| 2 (f^f fc (f e h %- 

- 80 5 (M 2 + M( 2 ) (g^LCi"*)** - 8^ M 2 (f|«)f fc (f e e %. 
-3(% + ^)</ 2 (m|) y 



+ 0g 



2(g /e *)l(g ,e «MB + 2(e)f fc (f e e «)^^J (mi),. 

+^ ( m |) ifc [2(g^)L(g ,e -) y ^+2(e)r i (e)r/ /l ' 
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Below the scale Q = inn, as discussed in Sees. 5.5 and 5.6.1, the trilinear couplings 
to the doublet h, and the mass parameter ml remain in the theory, with RGEs given by, 



(47T) 



~ dt 



(a u ) tfc - c U*f^) 1 l 2 ^ (c 2 - s 2 ) S kj + 2s 2 [(f u )t(f u 

V / iki o 



kj 



+ 29 h e Qk 



9 3 
^ _ 3ff| 
12 4 



<5ifc<% + 6(f !1 )ij(f u )[ 



/A: 



s(a. u ) k i - c[jl*i u 



kl 



{s 2 -c 2 )5 lk + 2s 2 [(f u )(f u )t 



(Wd) f 
S (a u ) fcj -c(/2*f tt ^) fc . 



';6 6 s (M 1 - %M[) 



32 



- - e- gS (m 3 - iM>z) (g?r ik (f u )ki(sr)h - ™ w e~ h s (m 2 - iM> 2 ) (~g h «y(g Q y ik (f: R )kj 

'3a 2 (fJ)«(£^), fc + c 2 {3(fJ)«(^) JJfc + (ft)^(f e )^}] [ s (a u ) y - c (>f£ M ) 



,hd \ +36* 



+ 0, 



^(f?)*(f?)L +^®*®L + ^(g ,Q )i*(g ,Q )S + fws )*^)*! 



+ 3*g(g?)tt(g?)£ 



- 3 



1 

36"' 



+ g^, + fh ) 9 



s{a u ) l3 -c(ffi^) 

12 i ^ 



4^ + N^ + 3 



S (a u )-c(/x%^) 
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(4*)' 



c{a d ) ij -s(jl*i^ d ) 



dt 



~- 26h9 Q k 



+ 6 d l Q k 



9 ^ + 3 f)(c 2 -s 2 )5 lk - S i (f„)(f u )t 



ik 



ik 



18 3 



Sikbij + §(U)\ k (U)ij 



+ 2r" (f d )(f d )t 

c(a d ) kj - s (/x*f 
C (ad) fc , - a (>f>) 



,12 



A-/. 



^ (c 2 - S 2 ) + 2c 2 [(f d )t(f d ) 



+ 2^^(f^(f e )[ fc 
+ [c( a<i ) ife - S (/i*f^). 

- 4^^*(fQ), fc (ft) fci (f^)^ + I^mV*" (2(f, Q ) lfc (g^)I, + (g' Q )* k (f d d R ) kj ) 



+ -^ c (M 1 -zM() 



32 



-^(^ d )*(g ,Q )* fe (f^)^ + ^g' Q M^)ki(g' dR )l, - 29^)^)10^ 
>~ g c {M 3 - iM$ (eSMUMeZ*)*, ~ ^wh c ( M 2 " iM 2) (~9 hd r(s Q )*k(^ R )kj 



+ ^%(g?)* fe (g?)^] [c(a d ) - - s (^ d ) lj 



l 'J. 



+ ^ 



w 



9 hd 



l g K \ 

c ( a d)ij - s (A*ff d Y. 



9 K 



ik- 



i 



i 



i 



36 <3i 9 d > 4 



j'2 



(% + <Ul + ^(%+^)</f 



c(a d ) - s (jS'f*-) 



J ij 
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(4tt) ; 



c(a e );j-s(/x*fe d ) . 



=29 h 9 



dt 

£ ~ 3j f) (c 2 ~ s 2 ) 6 ik + 2c 2 [(f e )(fe) t ] i J [c(a e ) fcj - s (jT&) 



kj 



+ 2(f e ); fc (f e ) ij ] [c(a e ) fe , - s (m*^) 



+ 6^^(f e ) ij (f d )J fc c( ad ) kl - S U*t d 



kl 



+ e h 9- ek [c(a e ) lfc - a (A*^) J -5' 2 (c 2 - * 2 ) <5 fcj + 2c 2 
+ 29 B c (Mi - iM[) 

- 69 w 6~ h c (M 2 - ^) (^ d )*(g L )* fc (f e efl )fc, 

+ [^(fe L )^(f e L )L + ^B(g ,L )I fc (g ,L )H + lOw^ L Mi L fk] [c(ae)y - - 
+ ^ [3s 2 (£^)fci(fJ)ifc + c 2 {3(f d ) fc ,(fi)i* + (SOwtfDifc} 



c(a e ) ij - a L*f** 



' B 



7,'hd 



+* 2 



+ 30* 



+ 0e k 

- 3 



<a e ) ij - 8 (j?f**) 

c(a e ) ife - a (>f e " d ) ] [2^(g' efl )^(g /efl )r, + 29^)1^ 



Kl 



\e u + gj + \e h ^ g' 2 + \ (9 Li + 9 h ) 5 2 } [c(a e ) - a 



ki 
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. 2 d\s 2 (m 2 Hu +\p,\ 2 )+c 2 (m* +\H\ 2 )-sc(b+b*) 
(4vr) -i ^ ' 



^ [</ 2 + <?1] (c 2 - s 2 ) 2 [, 2 (<, + |A| 2 ) + c 2 (rr& d + |/i| 2 ) - so (6 + 6*) 
+ fifc fil [-2</ 2 ( S 2 - c 2 ) % + 6s 2 [(f M ) T (f u )*] J ( m l) M 

+ ^Qk^Qi W 2 (** " + 6^ [(WQ T ] lk + 6c 2 [(W^HJ W) 

+ V* (« 2 " ^ 2 ) + 6c 2 [(f,) T (f,)1J (mi,) H 
+ ^ [V (« 2 " c 2 ) S lk + 2c 2 [(f e )*(f e ) T ] J (m 2 ) 
+ e- ek 9 h [g 12 {s 2 - c 2 ) 5 lk + 2c 2 [(fe) T (fe)*] J (ml 



kl 



'kl 



+ ze Lk e~ ei 



ik 



c(a d )i k - s (/i*f? d 



Ik 



^ e ) lk -s(jf^) ik 



c(ad)i, - s (/TfJ 
c(a e )L - a (>f e ^) 



+ 30 



2 



~9 hd 



-2B- h \0 § (M 2 + Mf) 

+3 V (M 2 + M' 2 ) 

~ \ {-M~ h 0B sc ^~9 hu ~9' hd (Mi + iM[) - YlB~ h d^sciL*~g hu t d (M 2 + iM^} 
^ {-40^/^)^)* (Mi - iM[) - ne-^scfiCg^rCg^y (M 2 - iM^)} 
S 2 (m^ + |A| 2 ) + c 2 (m^ + |A| 2 ) - S c (6 + b*) 



S/ 2 
2 2 



+ 3e w e h 



9 hd 



+c 2 \[Mi?J + 2rjT] kk + e § 6~ h 

x s 2 (m 2 Hu + |/i| 2 ) + c 2 (m 2 ^ + IAI 2 ) - sc (b 



(B.24) 



169 



With Oh = 0, RGEs for the remaining SSB scalar mass parameters take the form, 



dt 



~-0 h 



\ (c 2 - S 2 ) g>% 3 + 2. 2 [(f u )*(f u f] .. + 2c 2 [(f,)*(f d ) T ] .. 

x > (m^ + |/2| 2 ) + c 2 (m^ + |/i| 2 ) - sc (b + 6*) 
1 

g-ufc.? ~ tJ V v )kk + g 



2 1 

;6u k g' 2 5ij {™-u) kk + o e Q k 9 2 ^j ( m g) 



+ On. 0, 



Qk Qi \ 



-/2 



6^ (m 2 Q ) 



M 



+ \ e dJ% ~ lh k 9% K) fefc + \e~ej 2 k 3 (m|) 



kk 



+ 20^(f u )* fe (m&) H (f u )J + MtfzWik (m?,) H ft) 



+ 20 Jfc 0„ [c(a d )* fe - s (>f" d )*J [c(a,)^ - s (>f>) 

- ^ (M 2 + Mf ) (g'°k(g'% - 6^ (M| + M^ 2 ) (g^ fc (g% 
-f%(M 2 + M^ 2 ) (g?) tt (g?)t. 

- 48- h h 2 + (f d Q )* fc (f, Q )r; 
- 3 (% + ^) (^ /2+ ^ 2+ ^ 2 ) «), 



+ 



^Mg' Q U(g /Q )L + §Ws°)*(s g )L + |^(g?)«k(g?)L 



+ °Q k «) 



+|fls(g?)«(g?)!,-* + %(f u Q )^(f u °)S + '*g?)h«?)5] . 

(B.25) 



170 



-0 h [\(c 2 -s 2 )g'% 3 +As 2 [(f u ) T (f u )%) 

x ~s 2 (m 2 Hu + |/i| 2 ) + c 2 (m 2 Hd + |/2| 2 ) - SC (6 + 6*) 



- 12 i 2 

75 +5 3 



(m 2 

^Q fc 5% - ^<?% (m 2 D ) fefc + k fc5 % (m 2 ) 



kk 



(a u )f fe - c (>f^J Ua„)^- - c (>f u ^)* . 



on on 

- -g 9 6 ( M i + O (g' Mfl )L(g ,Ufl )^ " y e § (^1 + <) (g^)L(g^) fe , 



9 



3 5 3 ) (%) i3 



+ 0« 



ML 
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'/+-^(c 2 - S 2 ) 5 %+4c 2 [(fc) r (fc)« 



i ij 



m% u + IAI 2 ) + c 2 (m 2 Hi + - sc (b + 6*) 



" o^ fc 5% {™l) kk + -JqJ 2 ^ {m 2 Q ) kk + (m 2 ,) 



fefc 



3 o!fe «i 



;5 /2 +4<? 3 2 



tfikfy (m^) fei - -0 Lk g' 2 5ij (m|) 



fcfc 



+ ,A fc <7% (m|) fcfc +4^ fe ^(f d )f fc (f,)^. (m 2 . 



kl 



+ Ae Qk e h 



c( ad )J k - s (ins* 



c(a d )* kj -s [fi*i d 



°-e s (Ml + m( 2 ) (g^)L(g'^) fcJ - f * g (m! + M'i) (sH k (g d s «) kj 

88 h \»\ 2 (t d d «m«T kj ~ 3 fa + d dj ) Q</ 2 + ^|) (m2,) y 
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(47T) 



(47T) 



2 <, K)j 



x a 2 (m^ + |/I| 2 ) + c 2 (m^ + |/i| 2 ) - sc(b + 6*) 
+ 2^ 5 % (m 2 ;)^ - 0^% (mj) fcfc - ^<?% (m 2 D ) fefc 



+ (ml) kk + fl £fc (? £l + ^f 2 J (mi) H 

- (m|) fefe + 20 gfc e ~ ( (f e )* fe (m|) H (f e )£ 



c(a e )* fc - a 



c(a e )£ - a (>f e ^ 



kj 



20z {Ml + M?) (g' L k(g'% - 6^ (M| + M?) {z L )ik{% L ) ] 



kj 



J(g /L )^(g' L )L^ + ?(g L )^(g L )L^ + (teTik&fkA 



K) 



0' 



(g /L ) fe /(g /L )!/s + l(g L Mg%o w + (f e L )^(f e L )5^ 
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=^{-2(c 2 - S 2 )« 7 % + 4c 2 [(fe) T (fe)%} 

x a 2 (m^ + |/i| 2 ) + c 2 (m^ + |A| 2 ) - ac (6 + &*)' 
- 4^</% {ml) kk + 29Q k g'% (m 2 ,)^ + 20^% (m 2 ,), 
" 20 Zfc s% (m 2 ) kk + 2^2% (m|) fefe 

+ 29e k ee l9 '%5 ik ( m |) w + ie ik e Li {i e )l (m|) H (f e )^. 



/fcfc 



c(a e )f fc - a 



c a. 



- 86> fi (M 2 + O (g' efl )L(g' efl ) fej " 8^ |„| 2 (f e e «)f fc (f e e «)* fci 
-3(% + ^ 3 ) 5 ' 2 (m|.) 



+ 0; 



2 (S "Vifc 



e «) fc ^ + 2(f e e «)f fe (f e e «)^ 



kl h 



ml 



+ ( m IL [2(g^)L(g^)^^ + 2(e)S(f e e «)? J ^ 
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Appendix C 

DR RGEs in the SM 



The following are the RGEs for the gauge couplings, Yukawa matrices, quartic 
Higgs coupling and Higgs vacuum expectation value of the Standard Model in the DR renor- 
malisation scheme. The conventions are as in Ref. [13] such that to convert to the conven- 
tions used in the body of the text we must use X u ,d,e 

The gauge couplings run as follows: 



u,d,e' 



dgi 
dt 



-h 



Si 



16tT 2 



Eh 



9t 



(16^ 2 r (16vr 2 ) z 



^ QfcTr 



yIy, 



(C.l) 



J v / k=u,d,e 

with t = ln[Q] and = 1,2,3 corresponding to the individual gauge groups. Using 

n g = \nf (nj is the number of active flavours), the constants are defined by 

{ h\ ( \ 



and 



hj = 



( 


V o 



V ^3 / 




22 
3 



' 3 > 



\ o J 



136 
3 





\ 



102 / 



/ 19 
' 15 



n„ 



44 \ 
15 



49 
3 



w 4 



V 30 



76 
3 



50 
10 



V o 



^ \ 

io u x 

13 o 

6 u 

0/ 



/ 17 1 3 \ 
' 10 2 2 ' 



3 1 
2 2 



V 2 2 / 
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The running of the Yukawa couplings is 



dY k _ 1 a (i) , 1 «(2) 



dt 16vr 2 ^ Y fc (16tt 2 ) 

The are 



Py, + tt^Py: (C.5) 



^ = Y u {1 (YtY B - Y+Y d ) + Y 2 1 - (^ + f 52 2 + 8si) l} (C.6) 
^ = Y d {1 (YjY, - YtY u ) + y 2 l - (i ff 2 + f 5 l + 8g!) l} (C.7) 
/3« = Y e {|YtY e + y 2 l - I (g* + l} , (C.8) 



where 



y 2 = 3Tr[Yt YJ + 3Tr[Yt Y d ] + Tr[YtY e ] . (C.9) 



and the are 



/a( 2 ) _ -vr / 2"Vt"V Vt V _ vt V V^V — iv^Y Vt V 4- Hvtv v^v 
h'Y ~ 1 til2 1 u 1 ti 1 « 1 « 1 u 1 u 1 d d 4 1 d 1 ii 1 u I « t 4 I d I c2 I d I a 



+ Y 2 (fYjY, - |Yt Y u ) - X 4l + §A 2 1 - 2A (3Y£ Y u + YjY d ) 

+ ( Wrf + ^ + 20fli) Yt Y u - - i, 2 + 20, 2 ) YjY, (CIO) 

+ y 4i + [(ife + h 9 ) gf - Mai + if ahl - - 2n 9 ) 9 I 

Py\ = Y d{f Y d Y d Y rf Y d ~ Y d Y d Y i Y M ~ i Y i Y « Y d Y d + T Y i Y M Y i Y « 

+ y 2 (f Yt Y u - f Y*Y d ) - X4 1 + §A 2 1 - 2A (sY^Y, + yt Y tt ) 

+ (If 5? + W^l + 20 5 |) Y d Y d " (g* 2 " M + 20 5 2 ) Y i Y tt (C- 11 ) 
+ Y 4l + [- (fo " ^ " "Mai + ft^f - (^ - 2n 9 ) g\ 

-("¥- fn g ) gl] 1} 
= Y e {|YtY e YtY e - |y 2 YtY e - X4 1 + §A 2 1 - 6AY^Y e + (|£ 5 2 

+ W^) Y e Y e + m + [(& + 1%) + | 5l 2 5 I (C12) 

-(^-2n s ) 52 4 ] 1}, 
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with 



X 4 =1 |3TY [YtY^YtyJ +3Tr 



Y^Y.YtY, 



+ Tr 



yty Y+Y 



u YY^Y^J } 

n = (|<7? + f <?1 + 28 5 f) TV [Yt Y„ 



+ (-5)5? + f5 2 2 + 28 53 2 )Tr 



yJy, 



yty 

M- o x- c 



(C.13) 
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The Lagrangian term for the Higgs quartic coupling is 

A / , x2 



(C.15) 



where 4> is the SM Higgs field. It is not necessary to find the running of the DR coupling, 
since the difference between DR and MS A will be equivalent to a higher order effect in the 
Yukawa coupling running. However the gauge and Yukawa couplings are in the DR scheme 
for ease of computation. The RGE in the MS scheme is 



dX 



MS 



1 



dt 



id 



(i) 



+ 



1 



16tt 2 a ms (1671-2 



72 ' \- 



(2) 



(C.16) 



with the one-loop (3s given by 

C = 12X m ~ (hi + 9 ^ 2 ) Ams + i (M + bhl + ai) + 4Y 2 Ams 



(C.17) 
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and the two-loop /3s given by: 

78A ms + 18 (hi + 3 s!) *m ~ i( 2 -f ~ 10 %) 92 - Wial 



A MS 



~2% (If + f n 9 ) 9t] Ams + " 8n 9 ) ff 2 6 - | + |n 9 ) , 2 , 2 4 

_ _9_ (239 , 40„ \ 4 2 _ 27_ (59 , 40„ \ 6 
25 I 24 + 9 n 9) 9ld2 125 1 24 + 9 n 9> 9l 



+ (-^ 5l 2 + 18<7 2 2 - 128 5 |) Tr |yJY u YJY m 
+ (f g\ + 18 5 2 - 128 53 2 ) Tr 
+ (-f5 2 + 6 52 2 )Tr 
+ AMs{(is? + ¥s2 2 + H253 2 )Tr 



yty Y+Y 

e e e e 



1^2 



Y T Y 



+ 



J_n 2 , 27 2 
10 #1 -r 2 92 



YlY e 



} 



+112 5 2 ) Tr ^ 

+ §01 { ("iff 2 + 21 52 2 ) Tr [Yt Y u ] + (| ff 2 + 9 5 2 ) Tr 
+ {-fgf + llff 2 ) Tr [Y|Y e ] } - 24A 2 ^y 2 - X m H 



Y^Y, 



+ eA^Tr 



Yt Y u YjY d ] + 20Tr [3 (yJ Y u Yt Y M Yt Y tt ) 



YTY e YTY e YTY e 



+3 (Yj Y d Yt Y d Yt Yd ) + ( 

YtY u (y^ + yJy,) YlY d 



- 12Tr 



where 



H =3Tr 



Y T Y Y T Y 



+ 3Tr 



Y^Y d Yt Y , 



+ Tr 



Y T Y Y T Y 
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To find the boundary condition on the Higgs quartic coupling, we must consider the match- 
ing of MSSM Lagrangian with the SM Lagrangian at the scale where the heavy Higgs 
particles decouple. Given that 

CMssMB-l(g' 2 +g 2 )(s 2 -c 2 ) 2 (\h\ 2 y , (C.20) 

for the rotated SM-like MSSM Higgs field, h, defined in (4.28), the value of A at the scale 
run can be found by setting h = <f>. 

The Higgs vacuum expectation value (vsm) is presented in the same way, with all 
terms in the DR scheme except vsm and A which are in the MS scheme: 



din v S M 
dt 



- 7 (i) + 



l 6vr 2 '%s (167r2 y 



(2) 



MS 



(C.21) 
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with 



(C.22) 



(2) =- 5 AL,-(i<7? + f 52 2 + 28 53 2 )Tr 



1v ms 2 "MS 



u u 



-{-i- gl + 2 igl + 2%gl)Ti 



+ X4 



(C.23) 
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Appendix D 



Non- Trivial Check on the 



Coefficient of g 3 in the DR SM 

RGEs 



This is a non-trivial check on one aspect of the SM RGEs in the DR scheme. 
We will compare our result for the DR SM RGE of the bottom Yukawa coupling to the 
DR running from Ref . [42] , but in principle our procedure can also be used to check the top 
and tau Yukawa couplings. In this appendix we only check the coefficient of 53, however it 
is expected that this is the dominant contribution to the SM RGEs in this region. 

If we set 51 = 52 = A u = Ad = A e = 0, the general form of the running of the 
strong coupling is 

-//;! • s g %. (d.i) 



dt 

Similarly, the running of a general Yukawa coupling under the same conditions can be 
represented as 

Y a d t = a d + bgl. (D.2) 
We need to rearrange this to obtain the coefficients P and X in 

Aa = A°[ 5 f(t)] P [l + A^3 2 ] , (D.3) 

where A a is the coupling at some scale related to A^, the coupling at another scale. This 
form is the same as that in Ref. [42]. 
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Rearranging (D.2) we get 
d\ a 



dt 



— = [agl + bgi] j-dgs 



_7#f + 5 9l 
1 

93 



d§3 



+ 



dgs ■ 



(D.4) 



The two terms in the square brackets can be expanded as power series in to obtain (after 
multiplying through the I/53) 

d\„. ' a 
An, 1.7 \93 



{ - S ~9s + + ' b - 0(sg)) } <to. • (D.5) 



The next term in the series (cgf) would have contributed 



1 

53 



4 

c #3 



7 + ^3 



o(gl) 



so this expansion is only complete up to O(gs) since we do not keep three-loop terms. 
Going back to (D.3), we can take the natural log of either side 



In A a = In A° + P In g 2 3 (t) + In (1 + Xg 2 3 ) , 
which leads to an expression which has the same form as (D.5): 



dg. 



2P 



93 1 + X 9l Y 93 

Comparing this to (D.5) it can be seen that 



+ 2gzX(l-Xgl) 



dg3 



2P 



and 





b 


2X = = -\ 




T 


7 



(D.6) 



(D.7) 



(D.8) 



(D.9) 



These are valid for all (3, 3) entries of the Yukawa couplings and are also independent of 
the renormalisation scheme. 
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Focusing on the bottom Yukawa coupling, and taking n g = § in the RGEs (C.l) 
and (C.5), leads to 



16vr 2 a = - 8 



(I6vr 2 ) 2 6 = - 

16^ 2 7 = - 11 + 
(16tt 2 ) 2 <5 = - 102 + 



m - l l . 5 -) _ - 8 4 

V 3 3 2J 
4 5 _ _23 
3 2" ~T 
76 5__116 

IT 3~ ' 



(D.10) 
(D.ll) 

(D.12) 
(D.13) 



which, due to the choice of n g , are correct at a scale between rat and m&. Combining (D.8) 
with (D.10) and (D.12) gives 



2P = 8 



12 



23 23 
Similarly, (D.9) combined with (D.10)-(D.13) gives 



2X = 



16vr 2 



116 9 13^1 753 

I OA _v V — 

3 '529 16vr 2 ' ' 23 4vr 2 1058 



(D.14) 



(D.15) 



The equation from Ref. [42] that is equivalent to (D.3), and is valid in the region chosen 
above is 



Aa = X° a [g 2 3 (t)] 



12/23 



1 753 2 

i + — — gl 



(D.16) 



4vr 2 1058' 

Both P and X are seen to be the same as derived in (D.14) and (D.15). 

It is possible to go back to the boxed equations, (D.8) and (D.9), and calculate 
these again in the MS scheme, the RGEs of which are well known [14]. The only difference 
to (D.10)-(D.13) is that 

1019 

(D.17) 



(1^)'»=-T- 



This leads to the same P and a new X which is 

1 3731 



(D.18) 



4tt 2 3174' 

Correspondingly, the MS running of the bottom Yukawa coupling according to Ref. [42] is 



Aa = A° [gl(t)] 



12/23 



1 + 



1 3731 



4vr 2 3174 

and this supports the method used to check the DR RGEs. 



g?, 



(D.19) 
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Appendix E 

A Sample Input File for RGEFLAV 



The text that follows is a sample input file for RGEFLAV. In this example, the code 
will be run with a non-zero phase for the KM matrix (Question 2.) so that most numerical 
inputs must be in complex notation, i.e., (a,b), where a and b are real numbers, as can be 
seen: 

1. USE FULL TWO LOOP RUNNING? YES=1, N0=0 
1 

2. USE COMPLEX RUNNING WITH COMPLEX KM MATRIX? YES=1, N0=0 
1 

3. IF COMPLEX RUNNING , COMPLEX PHASE OF MU? YES=1, N0=0 


IF YES, THETA= : I.E. MU= I MU I *EXP ( iTHETA) 
3.142D0 

4. USE ISAJET'S THREHSOLDS? YES=1, N0=0 


5. USE ALL mSUGRA GUT CONDITIONS? YES=1, N0=0 
1 

IF YES GO TO 17, IF NO: 

6. UNIFIED GUT SCALE? YES=1, N0=0 
1 
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23 

24 IF NO, VALUE OF HIGH SCALE FOR SUSY INPUT: 

25 1.D19 

26 

27 7. USE mSUGRA GUT CONDITIONS FOR Ml, M2, M3? YES=1, NO=0 

28 1 

29 

30 IF NO, Ml, M2, M3 ARE: NB FOR COMPLEX RUNNING THESE ARE (M,M') 

31 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

32 

33 8. USE mSUGRA GUT CONDITIONS FOR M_{H_U}, M{H_D}? YES=1, N0=0 

34 1 

35 

36 IF NO, M_{H_U}, M_{H_D> ARE: 

37 (O.DO.O.DO) (O.DO.O.DO) 

38 

39 9. USE mSUGRA GUT CONDITIONS FOR M_Q? YES=1, N0=0 

40 1 

41 

42 IF NO, USE M~2_q = M~2_{Q0} \1 + T_Q . 

43 M_{QO} IS REAL. IT IS: 

44 O.DO 

45 

46 T_q MUST BE HERMITIAN. IT IS: 

47 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

48 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

49 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

50 

si 10. USE mSUGRA GUT CONDITIONS FOR M_U? YES=1, N0=0 

52 

53 

54 IF NO, USE 

55 M~2_U = M~2_{U0}- [ c_U \1 + R_U f"T_u f~*_u + S_U (f"T_u f"*_u)"2 ] + T_U. 

56 M_{UO} IS REAL. IT IS: 

57 . DO 

58 

59 C_U IS EITHER 1 OR 0. IT IS: 

60 
61 

62 R_U IS REAL. IT IS: 

63 O.DO 
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64 

65 S_U IS REAL. IT IS: 

66 O.DO 

67 

68 T_U MUST BE HERMITIAN. IT IS: 

69 (O.DO, O.DO) (O.DO, O.DO) (O.DO, O.DO) 

70 (O.DO, O.DO) (O.DO, O.DO) (O.DO, O.DO) 

71 (O.DO, O.DO) (O.DO, O.DO) (O.DO, O.DO) 

72 

73 11. USE mSUGRA GUT CONDITIONS FOR M_D? YES=1, N0=0 

74 

75 

76 IF NO, USE 

77 M~2_D = M~2_{D0> [ c_D \1 + R_D f~T_d f~*_d + S_D (f~T_d f~*_d)~2 ] + T_D. 

78 M_{DO} IS REAL. IT IS: 

79 O.DO 

80 

si C_D IS EITHER 1 OR 0. IT IS: 

82 

83 

84 R_D IS REAL. IT IS: 

85 O.DO 

86 

87 S_D IS REAL. IT IS: 

88 O.DO 

89 

90 T_D MUST BE HERMITIAN. 

91 (O.DO, O.DO) (O.DO, O.DO) 

92 (O.DO, O.DO) (O.DO, O.DO) 

93 (O.DO, O.DO) (O.DO, O.DO) 

94 

95 12. USE mSUGRA GUT CONDITIONS FOR M_L? YES=1, N0=0 

96 1 

97 

98 IF NO, USE M~2_L = M~2_{L0} \1 + T_L. 

99 M_{LO} IS REAL. IT IS: 
ioo O.DO 

101 

102 T_L MUST BE HERMITIAN. IT IS: 

103 (O.DO, O.DO) (O.DO, O.DO) (O.DO, O.DO) 

104 (O.DO, O.DO) (O.DO, O.DO) (O.DO, O.DO) 



IT IS: 
(O.DO, O.DO) 
(O.DO, O.DO) 
(O.DO, O.DO) 
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105 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

106 

107 13. USE mSUGRA GUT CONDITIONS FOR M_E? YES=1, N0=0 

108 1 

109 

no IF NO, USE 

in M~2_E = M~2_{E0> [ c_E \1 + R_E f~T_e f~*_e + S_E (f~T_e f~*_e)~2 ] + T_E. 

112 M_{EO> IS REAL. IT IS: 

113 O.DO 

114 

115 C_E IS EITHER 1 OR 0. IT IS: 

116 
117 

118 R_E IS REAL. IT IS: 

H9 O.DO 

120 

121 S_E IS REAL. IT IS: 

122 O.DO 

123 

124 T_E MUST BE HERMITIAN. IT IS: 

125 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

126 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

127 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

128 

129 14. USE mSUGRA GUT CONDITIONS FOR a_u? YES=1, N0=0 

130 1 
131 

132 IF NO, USE a_u = f_u [ A_{uO} \1 + W_u f~\dagger_u f_u 

133 + X_u ( f~\dagger_u f_u )~2 ] + Z_u. 

134 A_{uO} IS: 

135 (O.DO.O.DO) 

136 

137 W_u IS: 

138 (O.DO.O.DO) 

139 

140 X_u IS: 

141 (O.DO.O.DO) 

142 

143 AND Z_u IS: 

144 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 

145 (O.DO.O.DO) (O.DO.O.DO) (O.DO.O.DO) 
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146 


(O.DO.O 


.DO) (O.DO.O. DO) (O.DO.O. DO) 


147 






148 


15. USE 


mSUGRA GUT CONDITIONS FOR a_d? YES=1, N0=0 


149 


1 




150 






151 


IF NO, 


USE a_d = f_d [ A_{DO} \1 + W_d f~\dagger_d f_d 


152 




+ X_d ( f~\dagger_d f_d )~2 ] + 


153 


A_{DO} 


IS: 


154 


(O.DO.O 


.DO) 


155 






156 


W_d IS 




157 


(O.DO.O 


.DO) 


158 






159 


X_d IS 




160 


(O.DO.O 


.DO) 


161 






162 


IF NO, 


Z_d IS: 


163 


(O.DO.O 


.DO) (O.DO.O. DO) (O.DO.O. DO) 


164 


(O.DO.O 


.DO) (O.DO.O. DO) (O.DO.O. DO) 


165 


(O.DO.O 


.DO) (O.DO.O. DO) (O.DO.O. DO) 


166 






167 


16. USE 


mSUGRA GUT CONDITIONS FOR a_e? YES=1, N0=0 


168 


1 




169 






170 


IF NO, 


USE a_e = f_e [ A_{eO} \1 + W_e f~\dagger_e f_e 


171 




+ X_e ( f~\dagger_e f_e )~2 ] + 
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A_{eO> 


IS: 


173 


(O.DO.O 


.DO) 
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175 


W_e IS 




176 


(O.DO.O 


.DO) 


177 






178 


X_e IS 




179 


(O.DO.O 


.DO) 


180 






181 


Z_e IS 




182 


(O.DO.O 


.DO) (O.DO.O. DO) (O.DO.O. DO) 


183 


(O.DO.O 


.DO) (O.DO.O. DO) (O.DO.O. DO) 


184 


(O.DO.O 


.DO) (O.DO.O. DO) (O.DO.O. DO) 
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186 


17. OUTPUT IN BASIS WHERE UP QUARKS ARE DIAGONAL? 



187 



187 1 

188 

189 18. USER INPUTS ROTATION MATRICES ? YES=1, NO=0 

190 1 
191 

192 19. IF YES, GO TO 20, IF NO, V~U_L = KM ? YES=1, N0=0 

193 1 

194 

195 20. ARE V~U_R AND V~D_R IDENTITY? YES=1, N0=0 

196 1 
197 

198 21. IF USER INPUTS V'S, THEN THE ANGLES AND PHASE IN THE UNITARY 

199 MATRICES ARE: (NB IF KM IS NOT COMPLEX THE PHASE IS IGNORED) 

200 

201 V~U_L (ALPHA, BETA, GAMMA, DELTA) 

202 2.053D0 0.254D0 2 . 03D0 0.4829D0 

203 

204 V~U_R (ALPHA, BETA, GAMMA, DELTA) 

205 1.188D0 2.218D0 .763D0 0.87D0 

206 

207 V~D_R (ALPHA, BETA, GAMMA, DELTA) 

208 1.904D0 2.947D0 1 . 847D0 1 . 14D0 

209 
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Appendix F 



Numerical Instabilities Associated 
with Matrix Diagonalisation 

The problem 

We have emphasized that in order to properly implement particle decoupling into 
the RGEs, we have to be in the mass basis of the particles being decoupled. Our procedure 
for decoupling squarks, therefore, requires us to evaluate the unitary transformation from 
the given current basis to a new current basis that coincides with the squark mass basis 
(approximated, as discussed in the main text, to be the basis in which the SSB squark 
mass squared matrices are diagonal). Below the scale Q where at least one squark has 
decoupled, we not only have the rotations V l,r(u, d) (unitary matrices by construction) 
which connect the current basis with the basis in which the Yukawas are diagonal at mt, 
but also the squark rotations R, which connect the current basis to the "squark mass basis" 
that are obtained by numerically diagonalizing the SSB matrices m^. R. is of course the 
matrix of the orthogonal eigenvectors of m^. If there is a degeneracy of eigenvalues, the 
orthogonal eigenvectors are not uniquely defined. This leads to a practical problem when 
we numerically solve for the eigenvectors in the case that two eigenvalues of any SSB squark 
mass matrix with large off-diagonal components are degenerate to within ~ 1%. In this 
case, the corresponding eigenvectors, because of (system-dependent) numerical errors are 
not exactly orthogonal, and the corresponding matrix R, is not precisely unitary.''' 

^Using the g77 FORTRAN compiler with Macintosh Intel Macbook, together with the subroutine CG in the 
EISPACK collection of subroutines, we found that RlR. deviated from 11 to about one part in 10 10 compared 
to a part in 10 18 for Vl,r(u, d)' V l,r(u, d) for ~Vl,r{u, d) of the form (6.1). We obtain a similar size deviation 
from identity using the subroutine ZGEEV in the LAPACK collection of subroutines 
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The deviation from unitarity is very small, a part in 10 in our case, but is 
nonetheless orders of magnitude larger than what we can tolerate when calculating the 
smallest off-diagonal elements of m 2 . To understand why our calculation is sensitive to 
this seemingly tiny level of noise, let us imagine what would happen if we attempted to 
evolve the off-diagonal elements of mfj from Q = Mqtjt in a basis where the Yukawa 
coupling matrices all have large off-diagonal elements at the GUT scale. (This is not what 
we actually do, but we could imagine doing so since we know that we are well above all 
SUSY and Higgs field thresholds where the choice of basis should be irrelevant.) In the 
mSUGRA framework, the squark mass matrices are all given by m 2 = TOqII at Q = Mqut 
in any basis. Then, from (B.17) we see that these would develop off-diagonal components 
~ few x f 2 x 77Jq ~ f 2 x 4 x 10 4 for the case shown in Fig. 7.12, where f 2 denotes the 
size of the off-diagonal element of f Jf * . In this rough estimate we have assumed that the 
loop factor l/(167r 2 ) is compensated for by the large logarithm. In the general current basis 
where f u has comparable off-diagonal and diagonal elements, f 2 ~ 1, and the magnitude of 
the off-diagonal elements of mfj are O(10 4 ) GeV 2 . Rotating to our standard current basis 
should yield the result in Fig. 7.12. In particular, we should obtain [m 2 ;^ ~ 10~ 9 GeV 2 
because there would be large cancellations arising from the unitarity of Hjj that would 
suppress this matrix element. If instead the unitarity of R holds only to a part in 10 10 
because of numerical errors in obtaining the eigenvectors, we will find that because the 
cancellations are not perfect all off-diagonal elements of mfj will have a magnitude that is 
at least few x 100 2 x 10 -10 ~ few x 10 -6 GeV 2 , much larger than the magnitude of the 
(1,2) element in Fig. 7.12. We note here that the noise that leads to the non-unitarity of 
Vl,r(u, d) matrices at the 10 -18 level is completely irrelevant. 

The solution 

The non-unitarity of R, is only an issue when the off-diagonal entries of the squark 
mass matrix, in the basis where the Yukawas are diagonal at m t , are small compared to 
the diagonal entries. Since we, therefore, only need to consider matrices that are already 
approximately diagonal, we can associate the eigenvectors, (ei, e2, 63), with the approximate 
eigenvalues ((m 2 )n, (m 2 )22, (m 2 ,^), respectively. As an illustration, let us take a case 
where the (m 2 )23 entry is the off-diagonal entry with the largest magnitude, and (m 2 )i2 
the one with the smallest. We know the ordering quite unambiguously because above all 
squark thresholds we do not need to rotate by the matrices R, that potentially are the 
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mSUGRA: m =200 GeV, m 1/2 =-400 GeV, A Q =-200 GeV, tan|3=10, u>0 
Real variables and rotations only 




I I i I i I i I i I i I i I i ILL 

le+02 le+04 le+06 le+08 le+10 le+12 le+14 le+16 

Q (GeV) 



Figure F.l: The scale-dependence of the (1,2) and (2,1) elements of the Hermitian mfj matrix for our 
sample mSUGRA point (in the basis specified by (7.2)) calculated using two different procedures discussed in 
the text. The dashed (black) line shows the magnitude of the smallest element, i.e., | (mfr) | = | (mf/) 21 |, 
when we have used our procedure to ensure that the corresponding eigenvectors are orthogonal. The solid 
(red) line shows the same element when we do not pay attention to the orthogonality of the eigenvectors of 
raf r . The lighter dot-dashed (green) line shows the magnitude of the (1,3) element, and provides a scale for 
the size of the numerical noise discussed in the text. The noise in this curve is too small to be visible. All 
elements are zero at the GUT scale. 

origin of the noise. We need to ensure that the (m^)i2 entry does not suffer from any 
numerical noise due to the diagonalisation. This leads us to fix ei • e2 = and move any 
non-orthogonality of the eigenvectors into e2 • e3 so that the noise moves to (1x1^)23, the 
off-diagonal element with the largest magnitude. To accomplish this, we slightly modify (by 
parts in 10 10 , the limit of accuracy of the diagonalization routines) only the eigenvector e2 
from its value as given by the diagonalization routines, thereby leaving ei • e3 unaffected. 

To completely clarify what we have just described, although (as we have already 
stated) we do not need to rotate to the squark mass basis until we reach the highest squark 
threshold, we plot, in Fig. F.l, the result for |m^|i2 obtained by the two different methods 
mentioned above, in the basis where the up- type Yukawas are diagonal at mt, over the 
whole range Mz < Q < Mgut- The dashed (black) line shows the result where we have 
no rotation by Rry from Q = A/gut until the highest squark threshold, beyond which 
we implement our method for ensuring that the error from the non-orthogonality of the 
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eigenvectors of mfj only shows up in the (2,3) element. Indeed we see that this curve is 
smooth over its entire range of Q. The solid (red) line shows the result of carrying out 
the squark rotation without our fix of the eigenvectors over the entire range of Q. Note 
that there is significant noise all the way down to the low Q region where only some of 
the squarks have decoupled. This noise is largest at Q = Mgut where the eigenvalues of 
mfj are degenerate, and settles down to 1CT 5 GeV 2 , not far from our estimate above. The 
important thing is that the frozen value of this element is significantly different in the two 
cases, as a result of this noise, just before squark decoupling. It is for this range of Q (where 
the mass matrices that enter flavour-changing processes involving squarks will be evaluated) 
that we must reduce the numerical error as far as possible. The magnitude of the (1,3) 
element of mfj is shown for comparison by the dot-dashed (green) curve. It has no visible 
noise because the corresponding eigenvalues are sufficiently split, and the corresponding 
eigenvectors are orthogonal to a very high accuracy. 

The reader will be struck by the fact that the random downward fluctuations in 
the solid curve are roughly bounded by the dashed (black) curve which shows the correct 
magnitude of the matrix element. The reason for this is that the fluctuations whose typical 
magnitude is ~ 10 -5 GeV 2 need to randomly fluctuate down by four orders of magnitude to 
even reach the dashed (black) line, and even more to go below, the chance for which is very 
small. Indeed it is because we have shown results for the case where the SSB squark mass 
matrices are real that we see these fluctuations go down to even the level of the dashed 
(black) line. For the more general case the chance for both the real and the imaginary part 
of any matrix element to simultaneously fluctuate downward by this large magnitude is very 
small, so that the calculated magnitude (not shown here) is always larger than 10 -6 GeV 2 . 

After our fix of the eigenvectors, any error from the non-unitarity of Hjj is shifted 
to the largest off-diagonal element, and the only residue of the resulting noise that remains 
is in the magnitude of this element for scales close to Mqut — where the eigenvalues are 
closest — as seen in Fig. F.2. At lower scales, the eigenvalues split, and the noise level 
(whose magnitude remains the same as in Fig. F.l) becomes insignificant. 

Again, the solid (red) line shows the evolution of the magnitude of the (2,3) 
element before fixing the eigenvectors, while the dashed black line shows the same thing 
after the orthogonality fix. We see that the numerical noise has indeed moved to the (2,3) 
element which now shows fluctuations, but only close to Mgut, where the eigenvalues of 
mfj are roughly degenerate. 
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mSUGRA: m ( =200 GeV, m 1/2 =-400 GeV, A Q =-200 GeV, tan(3=10, u>0 



0.0015 




le+15 le+16 

Q (GeV) 

Figure F.2: Scale dependence of the magnitude of the (2, 3) and (3, 2) entries of mfj for the same mSUGRA 
point as Fig. F.l. We focus on the running at the extreme high scale, and compare the noise in the magnitude 
of just this element, both before and after fixing the orthogonality of the eigenvectors as described in the 
text. As in Fig. F.l, the solid (red) line shows the result before the orthogonality fix while the dashed (black) 
line shows the result after this fix when this error has been moved to the (2, 3) element which now randomly 
fluctuates close to Mgut where the eigenvalues of mfj are roughly degenerate. 
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Glossary 



CERN European Organization for Nuclear Research 

EWSB Electroweak Symmetry Breaking 

GUT Grand Unified Theory 

ISAJET a Monte Carlo particle event generator 

KM Kobayashi-Maskawa (mixing matrix) 

LHC Large Hadron Collider 

LSP Lightest Supersymmetric Particle 

MSSM Minimal Supersymmetric Standard Model 

mSUGRA Minimal Supergravity 

RGE Renormalisation Group Equation 

SM Standard Model of particle physics 

SSB Soft SUSY-Breaking 

SUSY Supersymmetry 

VEV Vacuum Expectation Value 

WMAP Wilkinson Microwave Anisotropy Probe 
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